A = mxCreateDoubleMatrix(M, N, mxREAL);mxDestroyArray(A);

原创 2012年03月21日 14:30:10
#include <string.h>
#include "mex.h"
void DisplayMatrix(char *Name, double *Data, int M, int N)
{ /* Display matrix data */
	int m, n;
	mexPrintf("%s = \n", Name);
	for(m = 0; m < M; m++, mexPrintf("\n"))
		for(n = 0; n < N; n++)
			mexPrintf("%8.4f ", Data[m + M*n]);
}

void CallQR(double *Data, int M, int N)
{ /* Perform QR factorization by calling the MATLAB function */
	mxArray *Q, *R, *A;
	mxArray *ppLhs[2];
	DisplayMatrix("Input", Data, M, N);
	A = mxCreateDoubleMatrix(M, N, mxREAL); /* Put input in an mxArray */
	memcpy(mxGetPr(A), Data, sizeof(double)*M*N);
	mexCallMATLAB(2, ppLhs, 1, &A, "qr"); /* Call MATLAB's qr function */
	Q = ppLhs[0];
	R = ppLhs[1];
	DisplayMatrix("Q", mxGetPr(Q), M, N);
	DisplayMatrix("R", mxGetPr(R), M, N);
	mxDestroyArray(R); /* No longer need these */
	mxDestroyArray(Q);
	mxDestroyArray(A);
}

void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[])
{
#define M_IN prhs[0]
	if(nrhs != 1 || mxGetNumberOfDimensions(M_IN) != 2 || !mxIsDouble(M_IN))
		mexErrMsgTxt("Invalid input.");
	CallQR(mxGetPr(M_IN), mxGetM(M_IN), mxGetN(M_IN));
}

A = mxCreateDoubleMatrix(M, N, mxREAL);mxDestroyArray(A);

#include #include "mex.h" void DisplayMatrix(char *Name, double *Data, int M, int N) { /* Display m...
  • autumn20080101
  • autumn20080101
  • 2012年03月21日 14:30
  • 2913

Given a matrix of m x n elements (m rows, n columns), return all elements of the matrix in spiral or

1. Given a matrix of m x n elements (m rows, n columns), return all elements of the matrix in spiral...
  • qq_18343569
  • qq_18343569
  • 2015年11月26日 20:23
  • 480

矩阵乘法---蓝桥杯

import java.util.Scanner; public class matrixMul{ public static void main(String[] args) { int N,...
  • kerrylovelj
  • kerrylovelj
  • 2014年10月27日 13:59
  • 1147

矩阵乘法(简单基础)

矩阵乘法  问题描述 给定一个N阶矩阵A,输出A的M次幂(M是非负整数) 例如:   A =   1 2   3 4   A的2次幂   7 10   15 22 输入格式 第一行是一个正整数...
  • qiqi_skystar
  • qiqi_skystar
  • 2016年01月23日 11:40
  • 1440

求C(n,m)%mod

long long pow_mod(int a,int n) { long long ans=1; while(n) { if(n&1) ans=ans*a%m...
  • acm_cxq
  • acm_cxq
  • 2016年07月11日 15:38
  • 583

N o v e m b e r 1 2 t h M o n d a y

I   h a v e   t o   s e p a r a t e   o u t   c o d e s   f r o m   t h e   w h o l e   p r o j e c ...
  • Lu_ming
  • Lu_ming
  • 2007年11月30日 21:27
  • 1845

HDU 2685-I won't tell you this is about number theory(gcd(A^m-B^m,A^n-B^n)化简)

题目地址:HDU 2685 题意:求 gcd(a^m-1,a^n-1)%k的值。 思路:根据公式gcd(A^m-B^m,A^n-B^n)= A^gcd(m,n) - B^gcd(m,n)。在此处B...
  • u013486414
  • u013486414
  • 2015年08月20日 20:53
  • 685

NOI10004 填充矩形

时间限制: 1000 ms 空间限制: 262144 KB题目描述已知矩形的大小为n×m,现用a×a的正方形填充该矩形。输入三个正整数n,m,a(n,m,a≤10^9),计算至多能填入多少正方形?(...
  • u011544909
  • u011544909
  • 2018年01月04日 14:57
  • 82

递归算法 m和A和n个B,求有多少种排列

import java.util.Scanner; public class pailie { public static void main(String[] args) { // TODO...
  • chzayi
  • chzayi
  • 2015年03月28日 14:30
  • 949

poj1426 Find The Multiple(哈夫曼思想)

Given a positive integer n, write a program to find out a nonzero multiple m of n whose decimal repr...
  • cacyth
  • cacyth
  • 2015年08月30日 23:07
  • 955
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:A = mxCreateDoubleMatrix(M, N, mxREAL);mxDestroyArray(A);
举报原因:
原因补充:

(最多只允许输入30个字)