Choose solver for the invariance equation
Contents
function [W] = solveinveq(C,b,solver)
SolveInvEq
solves C*W = R for W with given solver
switch solver case 'lsqminnorm' W = lsqminnorm(C,b); case 'linsolve' W = linsolve(C,b); case 'backslash' W = C \b; case 'pinv' W = pinv(C)*b; case 'inv' W = inv(C)*b; case 'bicg' W = bicg(C,b); case 'bicgstab' W = bicgstab(C,b); case 'cgs' W = cgs(C,b); case 'gmres' W = gmres(C,b); case 'lsqr' W = gmres(C,b); end
end