Addition of two subroutines that belong in LightKrylov that, as of this writing, are only used within LightKrylov:
zero_basis: helper function to zero out a Krylov basis. This could be used in the initialisation in examples/tests as well.
sqrtm: matrix-valued function for the sqrt of a SPD matrix. In the present implementation it is included only for real matrices (rsp + rdp) and includes checks for symmetry and positive definiteness. In principle, it could be extended to complex matrices with appropriate changes to the checks, but as of now it seems superfluous.
Note: The sval function referenced in the issue is not added. It is kept in LightROM.
Addition of two subroutines that belong in
LightKrylovthat, as of this writing, are only used withinLightKrylov:zero_basis: helper function to zero out a Krylov basis. This could be used in the initialisation in examples/tests as well.sqrtm: matrix-valued function for the sqrt of a SPD matrix. In the present implementation it is included only for real matrices (rsp + rdp) and includes checks for symmetry and positive definiteness. In principle, it could be extended to complex matrices with appropriate changes to the checks, but as of now it seems superfluous.Note: The
svalfunction referenced in the issue is not added. It is kept inLightROM.Closes #84