search

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
4 stars 1 forks source link

Encapsulate plane rotation #30

Closed jhabriel closed 10 months ago

jhabriel commented 12 months ago

QR rotation block must be encapsulated, i.e., lines 299 - 316.

Consider first whether the built-in function can be used


[R, g] = plane_rotations(H, beta)

% H: Upper Hessenberg matrix (this is an m+1-by-m matrix)
% beta: Norm of the residual (scalar)
% R: Upper triangular matrix (m-by-m)
% g: *decide on name* (m-by-1)
gusespinola commented 11 months ago

There is a recommendation for this encapsulation in the wiki for GMRES.

https://en.wikipedia.org/wiki/Generalized_minimal_residual_method