[paulzhang]

Not quite sure about how the matrix A is decomposed in SU2. Usually, a Gram-Schimit process is used, and an orthogonal matrix Q is first computed and then the upper triangular matrix R.

Can someone provide some details of it?

Thanks,
Paul

[hlk]

Originally Posted by paulzhang

Not quite sure about how the matrix A is decomposed in SU2. Usually, a Gram-Schimit process is used, and an orthogonal matrix Q is first computed and then the upper triangular matrix R.

Can someone provide some details of it?

Thanks,
Paul

Thanks for your question
I think this might be answered by the implementation in:
SU2/SU2_CFD/src/solver_structure.cpp. CSolver::SetSolution_Gradient_LS

[paulzhang]

Originally Posted by hlk

Thanks for your question
I think this might be answered by the implementation in:
SU2/SU2_CFD/src/solver_structure.cpp. CSolver::SetSolution_Gradient_LS

Got it. Thanks a lot for the enlightenment.
Paul