Skip to content
Snippets Groups Projects
Commit 72b1d692 authored by Carsten Gräser's avatar Carsten Gräser
Browse files

[test] Extend test of MutliTypeBlockMatrix::mv, umv, mmv, usmv

Test with non-zero matrix and print results
parent e0230a63
No related branches found
No related tags found
1 merge request!45Use hybridutilities
......@@ -71,10 +71,22 @@ int main(int argc, char** argv) try
result[_0].resize(3);
result[_1].resize(2);
multiMatrix[_0][_0] = 4200;
multiMatrix[_0][_1] = 4201;
multiMatrix[_1][_0] = 4210;
multiMatrix[_1][_1] = 4211;
multiMatrix.mv(multiVector,result);
std::cout << "mv result" << std::endl << result << std::endl;
multiMatrix.umv(multiVector,result);
std::cout << "umv result" << std::endl << result << std::endl;
multiMatrix.mmv(multiVector,result);
std::cout << "mmv result" << std::endl << result << std::endl;
multiMatrix.usmv(3.14,multiVector,result);
std::cout << "usmv result" << std::endl << result << std::endl;
///////////////////////////////////////////////////////////////////////////////////////////////////////////
// Next test: Set up a linear system with a matrix consisting of 2x2 sparse scalar matrices.
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment