Closed perazz closed 1 month ago
Fix #817.
for Ax=b with shape(A) == [m, n]:
Ax=b
shape(A) == [m, n]
size(x,1)>=n
size(x,1)>=m
*GESV
size(x,1)==n
cc: @jvdp1 @jalvesz @loiseaujc
Thanks a lot @jvdp1. I will wait another couple of days and then merge, if there are no further comments.
Fix #817.
for
Ax=b
withshape(A) == [m, n]
:size(x,1)>=n
.size(x,1)>=m
: use as temporary storage and do not reallocate*GESV
interfacesize(x,1)==n
cc: @jvdp1 @jalvesz @loiseaujc