Skip to content
Snippets Groups Projects
Commit 33b39391 authored by Christophe Geuzaine's avatar Christophe Geuzaine
Browse files

simplify transpose

parent 26221cd5
Branches
Tags
No related merge requests found
......@@ -532,14 +532,14 @@ static gmshMatrix<double> generateLagrangeMonomialCoefficients
for (int i = 0; i < ndofs; i++) {
for (int j = 0; j < ndofs; j++) {
double dd = 1.;
for (int k = 0; k < dim; k++) dd *= pow(point(i, k), monomial(j, k));
for (int k = 0; k < dim; k++) dd *= pow(point(j, k), monomial(i, k));
Vandermonde(i, j) = dd;
}
}
gmshMatrix<double> coefficient(ndofs, ndofs);
Vandermonde.invert(coefficient);
return coefficient.transpose();
return coefficient;
}
std::map<int, gmshFunctionSpace> gmshFunctionSpaces::fs;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment