-
Notifications
You must be signed in to change notification settings - Fork 707
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Use more LAPACK and BLAS routines in LAPACKFullMatrix #5934
Conversation
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Thanks for attacking this issue. Please add a changelog
|
||
// Multiply rectangular mxn real matrix by real scalar CTO/CFROM | ||
void dlascl_(const char *type, | ||
const int *kl, |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Please use dealii::types::blas_int
source/lac/lapack_full_matrix.cc
Outdated
|
||
lascl(&type,&kl,&kl,&factor,&cto,&m,&n,values,&lda,&info); | ||
|
||
Assert(info >= 0, ExcInternalError()); |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I think we throw some custom error for Lapack
source/lac/lapack_full_matrix.cc
Outdated
for (size_type i=0; i<m(); ++i) | ||
for (size_type j=0; j<n(); ++j) | ||
(*this)(i,j) += a * A(i,j); | ||
// BLAS does not offer functions to add matrices |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Full stop
99f399f
to
efab33b
Compare
efab33b
to
453f4e2
Compare
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Looks good thanks for the patch
/run-tests |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
thanks for updating.
This PR is a partial remedy to #5900.
The assignment operator for
FullMatrix
is not changed due to the different memory layouts ofFullMatrix
andLAPACKFullMatrix
.