MATLAB: Can you use bicg with initial guess but without preconditioner. If so, how do you call it

bicg

can you use bicg with an initial guess, x0, but without specifying a preconditioner M or M1, M2?

Best Answer

Hi Michael,

I believe you can specify the initial guess without specifying the preconditioner. In place of the preconditioner, just pass '[]'. For example,

n = 100; 
on = ones(n,1); 
A = spdiags([-2*on 4*on -on],-1:1,n,n);
b = sum(A,2); 
tol = 1e-8; 
maxit = 15;
x0=zeros([n 1]);
x = bicg(A,b,tol,maxit,[],[],x0);

NOTE: The above example is just to let you know that it is syntactically possible to pass the initial guess to the function "bicg" without the preconditioner. The initial guess that I have used above is purely random and does not have a correlation with the underlying mathematics.

I hope this answers your question.

Regards,

Vidya

Related Solutions

MATLAB: 5th root with Newton-Raphson code

If n is the number whose 5th root has to be calculated:
function f = fifth_root(n)
  xold=10;
  for i = 1:100
  xnew = xold - (xold*xold*xold*xold*xold-n)/(5*xold*xold*xold*xold);
  xold=xnew;
  end
  f=xnew
  end

MATLAB: Can you vectorize this code

n=3; 
x=meshgrid(1:n,1:n);
y=x';

Best Answer

Related Solutions

MATLAB: 5th root with Newton-Raphson code

MATLAB: Can you vectorize this code

Related Question