Another idea: allow more than one return value from functions. Then I think functions and subroutines would be equivalent (except intent(inout)
arguments). For eigensolver:
lam, c = eig(A)
Another idea: allow more than one return value from functions. Then I think functions and subroutines would be equivalent (except intent(inout)
arguments). For eigensolver:
lam, c = eig(A)