Conjugate Gradients Squared method

C Prototype

mxArray *mlfCgs(mxArray **flag,
                mxArray **relres,
                mxArray **iter,
                mxArray **resvec,
                mxArray *A,
                mxArray *b,
                mxArray *tol,
                mxArray *maxit,
                mxArray *M1,
                mxArray *M2,
                mxArray *x0,
                ...);

C Syntax

#include "matlab.h"

mxArray *A, *b;                 /* Required input argument(s) */
mxArray *tol, *maxit;           /* Optional input argument(s) */
mxArray *M, *M1, *M2, *x0;      /* Optional input argument(s) */
mxArray *flag=NULL,*relres=NULL;/* Optional output argument(s) */
mxArray *iter=NULL,*resvec=NULL;/* Optional output argument(s) */
mxArray *x = NULL;              /* Return value */

mlfAssign(&x, mlfCgs(NULL,NULL,NULL,NULL,
                     A,b,NULL,NULL,NULL,NULL,NULL,NULL));
mlfAssign(&x, mlfCgs(NULL,NULL,NULL,NULL,
                     A,b,tol,NULL,NULL,NULL,NULL,NULL));
mlfAssign(&x, mlfCgs(NULL,NULL,NULL,NULL,
                     A,b,tol,maxit,NULL,NULL,NULL,NULL));
mlfAssign(&x, mlfCgs(NULL,NULL,NULL,NULL,
                     A,b,tol,maxit,M,NULL,NULL,NULL));
mlfAssign(&x, mlfCgs(NULL,NULL,NULL,NULL,
                     A,b,tol,maxit,M1,M2,NULL,NULL));
mlfAssign(&x, mlfCgs(NULL,NULL,NULL,NULL,
                     A,b,tol,maxit,M1,M2,x0,NULL));
mlfAssign(&x, mlfCgs(&flag,NULL,NULL,NULL,
                     A,b,tol,maxit,M1,M2,x0,NULL));
mlfAssign(&x, mlfCgs(&flag,&relres,NULL,NULL,
                     A,b,tol,maxit,M1,M2,x0,NULL));
mlfAssign(&x, mlfCgs(&flag,&relres,&iter,NULL,
                     A,b,tol,maxit,M1,M2,x0,NULL));
mlfAssign(&x, mlfCgs(&flag,&relres,&iter,&resvec,
                     A,b,tol,maxit,M1,M2,x0,NULL));

MATLAB Syntax

x = cgs(A,b)
cgs(A,b,tol)
cgs(A,b,tol,maxit)
cgs(A,b,tol,maxit,M)
cgs(A,b,tol,maxit,M1,M2)
cgs(A,b,tol,maxit,M1,M2,x0)
x = cgs(A,b,tol,maxit,M1,M2,x0)
[x,flag] = cgs(A,b,tol,maxit,M1,M2,x0)
[x,flag,relres] = cgs(A,b,tol,maxit,M1,M2,x0)
[x,flag,relres,iter] = cgs(A,b,tol,maxit,M1,M2,x0)
[x,flag,relres,iter,resvec] = cgs(A,b,tol,maxit,M1,M2,x0)

See Also

MATLAB cgs     Calling Conventions