To receive notifications about scheduled maintenance, please subscribe to the mailing-list gitlab-operations@sympa.ethz.ch. You can subscribe to the mailing-list at https://sympa.ethz.ch

Commit 9acfa4f4 authored by Bastian Telgen's avatar Bastian Telgen
Browse files

remove hermitian transform

parent 3dfd9105
......@@ -232,9 +232,8 @@ int main(int argc, char **argv) {
.matrix;
}(wave_vector);
auto Tt = cpppetsc::asTransposedMatrix(&T);
const auto KT = cpppetsc::asThAT(&K, &Tt);
const auto MT = cpppetsc::asThAT(&M, &Tt);
const auto KT = cpppetsc::asThAT(&K, &T);
const auto MT = cpppetsc::asThAT(&M, &T);
const auto result = cppslepc::computeGeneralizedEigenvaluesWithSINVERT(
cpppetsc::computeElementsOfMatrix(KT),
......
......@@ -191,9 +191,8 @@ int main(int argc, char **argv) {
.matrix;
}(wave_vector);
auto Tt = cpppetsc::asTransposedMatrix(&T);
const auto KT = cpppetsc::asThAT(&K, &Tt);
const auto MT = cpppetsc::asThAT(&M, &Tt);
const auto KT = cpppetsc::asThAT(&K, &T);
const auto MT = cpppetsc::asThAT(&M, &T);
const auto result = cppslepc::computeGeneralizedEigenvaluesWithSINVERT(
cpppetsc::computeElementsOfMatrix(KT),
......
......@@ -225,9 +225,8 @@ int main(int argc, char **argv) {
}(wave_vector);
// Transform mass and stiffness matrix (eq. 13)
auto Tt = cpppetsc::asTransposedMatrix(&T);
const auto KT = cpppetsc::asThAT(&K, &Tt);
const auto MT = cpppetsc::asThAT(&M, &Tt);
const auto KT = cpppetsc::asThAT(&K, &T);
const auto MT = cpppetsc::asThAT(&M, &T);
// Solve the generalized eigenvalue problem (eq. 14)
const auto number_of_eigenvalues = Mesh::size_type(20);
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment