nidtec-una / krysbas-dev

A collection of Krylov Subspace-Based Adaptive Solvers (KrySBAS) written in MATLAB
https://krysbas-dev.readthedocs.io
GNU General Public License v3.0
2 stars 1 forks source link

Add possibility to acept matrix A with complex numbers #46

Open jccf19 opened 8 months ago

jccf19 commented 8 months ago

To do this. modify plane rotations:

        s = H(j + 1, j) / (sqrt(abs(H(j + 1, j))^2 + H(j, j)^2));
        c = H(j, j) / (sqrt(abs(H(j + 1, j))^2 + H(j, j)^2));