Remove an old commit and add a missing deallocation.

parent 9770c06a
......@@ -1211,9 +1211,7 @@ subroutine X(least_squares_vec)(nn, aa, bb, xx)
SAFE_ALLOCATE(ss(1:nn))
! MJV 2016 11 09 : TODO: this is callable with complex, but does nothing!!!
#ifdef R_TREAL
call lapack_gelss(nn, nn, 1, aa(1, 1), lead_dim(aa), xx(1), nn, ss(1), CNST(-1.0), rank, dlwork, -1, info)
SAFE_ALLOCATE(work(1:int(dlwork)))
......@@ -1229,6 +1227,8 @@ subroutine X(least_squares_vec)(nn, aa, bb, xx)
SAFE_DEALLOCATE_A(rwork)
#endif
SAFE_DEALLOCATE_A(ss)
if(info /= 0) then
write(message(1), '(5a,i5)') &
'In ', TOSTRING(X(lalg_least_squares_vec)), ', LAPACK ', TOSTRING(X(gelss)), ' returned error mess age ', info
......
Markdown is supported
0%
or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment