Create a sparse identifty matrix when it makes sense.
parent
fd76ce53af
commit
24fbc0d923
|
@ -123,7 +123,11 @@ for its = 1:maxit
|
||||||
fjac2=fjac'*fjac;
|
fjac2=fjac'*fjac;
|
||||||
temp=max(sum(abs(fjac2)));
|
temp=max(sum(abs(fjac2)));
|
||||||
if temp>0
|
if temp>0
|
||||||
p=-(fjac2+sqrt(nn*eps)*temp*eye(nn))\(fjac'*fvec);
|
if issparse(fjac)
|
||||||
|
p=-(fjac2+sqrt(nn*eps)*temp*speye(nn))\(fjac'*fvec);
|
||||||
|
else
|
||||||
|
p=-(fjac2+sqrt(nn*eps)*temp*eye(nn))\(fjac'*fvec);
|
||||||
|
end
|
||||||
else
|
else
|
||||||
errorflag = true;
|
errorflag = true;
|
||||||
errorcode = 5;
|
errorcode = 5;
|
||||||
|
|
Loading…
Reference in New Issue