Quasi-Minimal Residual method

C++ Prototype

mwArray qmr(const mwArray &in1,
            const mwArray &in2=mwArray::DIN,
            const mwArray &in3=mwArray::DIN,
            const mwArray &in4=mwArray::DIN,
            const mwArray &in5=mwArray::DIN,
            const mwArray &in6=mwArray::DIN,
            const mwArray &in7=mwArray::DIN,
            const mwVarargin &in8=mwVarargin::DIN,
            const mwArray &in9=mwArray::DIN,
            .
            .
            .
            const mwArray &in39=mwArray::DIN);

mwArray qmr(mwArray *out1, mwArray *out2,
            mwArray *out3, mwArray *out4,
            const mwArray &in1,
            const mwArray &in2=mwArray::DIN,
            const mwArray &in3=mwArray::DIN,
            const mwArray &in4=mwArray::DIN,
            const mwArray &in5=mwArray::DIN,
            const mwArray &in6=mwArray::DIN,
            const mwArray &in7=mwArray::DIN,
            const mwVarargin &in8=mwVarargin::DIN,
            const mwArray &in9=mwArray::DIN,
            .
            .
            .
            const mwArray &in39=mwArray::DIN);

C++ Syntax

#include "matlab.hpp"

mwArray A, b, tol, maxit, M, M1, M2, x0;// Input argument(s)
mwArray flag, relres, iter, resvec;     // Output argument(s) 
mwArray x;                              // Return value

x = qmr(A,b);
qmr(A,b,tol);
qmr(A,b,tol,maxit);
qmr(A,b,tol,maxit,M);
qmr(A,b,tol,maxit,M1,M2);
qmr(A,b,tol,maxit,M1,M2,x0);
x = qmr(A,b,tol,maxit,M1,M2,x0);
x = qmr(&flag,A,b,tol,maxit,M1,M2,x0);
x = qmr(&flag,&relres,A,b,tol,maxit,M1,M2,x0);
x = qmr(&flag,&relres,&iter,A,b,tol,maxit,M1,M2,x0);
x = qmr(&flag,&relres,&iter,&resvec,A,b,tol,maxit,M1,M2,x0);

MATLAB Syntax

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