]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
test for matrix copy between types
authorguido <guido@0785d39b-7218-0410-832d-ea1e28bc413d>
Tue, 3 May 2005 19:25:46 +0000 (19:25 +0000)
committerguido <guido@0785d39b-7218-0410-832d-ea1e28bc413d>
Tue, 3 May 2005 19:25:46 +0000 (19:25 +0000)
git-svn-id: https://svn.dealii.org/trunk@10636 0785d39b-7218-0410-832d-ea1e28bc413d

tests/lac/Makefile
tests/lac/matrices.cc [new file with mode: 0644]
tests/results/i686-pc-linux-gnu+gcc3.2/lac/matrices.output [new file with mode: 0644]

index 42fd73394d1ebd4faab2ef4484ef1c8b2f6e7e6e..d1f09eba66f0ad6b95f09f0feb438c5e1e88478b 100644 (file)
@@ -22,7 +22,7 @@ default: run-tests
 
 tests_x = vector-vector \
        full_matrix householder tridiagonal_matrix \
-       sparsity_pattern sparse_matrices \
+       sparsity_pattern sparse_matrices matrices \
        block_vector block_vector_iterator block_matrices \
        matrix_lib matrix_out \
        solver eigen \
diff --git a/tests/lac/matrices.cc b/tests/lac/matrices.cc
new file mode 100644 (file)
index 0000000..d0b719a
--- /dev/null
@@ -0,0 +1,72 @@
+//----------------------------  full_matrix.cc,v  ---------------------------
+//    $Id$
+//    Version: $Name$ 
+//
+//    Copyright (C) 2005 by the deal.II authors
+//
+//    This file is subject to QPL and may not be  distributed
+//    without copyright and license information. Please refer
+//    to the file deal.II/doc/license.html for the  text  and
+//    further information on this license.
+//
+//----------------------------  full_matrix.cc,v  ---------------------------
+
+// Test copying between matrices (copy_from and fill)
+
+#include "../tests.h"
+
+#include <base/logstream.h>
+
+#include <lac/full_matrix.h>
+#include <lac/lapack_full_matrix.h>
+#include <lac/sparse_matrix.h>
+#include <lac/sparse_matrix_ez.h>
+
+#include <fstream>
+
+int main()
+{
+  std::ofstream logfile("matrices.output");
+  logfile.setf(std::ios::fixed);
+  logfile.precision(3);
+  deallog.attach(logfile);
+  deallog.depth_console(0);
+  deallog.threshold_double(1.e-10);
+  
+  SparseMatrixEZ<double> ez(5,4);
+  ez.set(0,0,2.);
+  ez.set(0,2,3.);
+  ez.set(0,3,4.);
+  ez.set(1,0,5.);
+  ez.set(1,1,6.);
+  ez.set(1,3,7.);
+  ez.set(2,0,8.);
+  ez.set(2,1,9.);
+  ez.set(2,2,10.);
+  ez.set(2,3,11.);
+  ez.set(4,0,12.);
+  ez.set(4,2,13.);
+  ez.set(4,3,14.);
+  
+  deallog << "FullMatrix<float>::copy_from  SparseMatrixEZ<double>"
+         << std::endl;
+  FullMatrix<float> ff;
+  ff.copy_from(ez);
+  ff.print_formatted(logfile, 0, false, 5, "~");
+  
+  deallog << "LAPACKFullMatrix<double>::copy_from  SparseMatrixEZ<double>"
+         << std::endl;
+  LAPACKFullMatrix<double> lfd;
+  lfd.copy_from(ez);
+  lfd.print_formatted(logfile, 0, false, 5, "~");
+
+  lfd.reinit(2, 3);
+  deallog << "LAPACKFullMatrix<double>::fill  SparseMatrixEZ<double>"
+         << std::endl;
+  lfd.fill(ez, 0, 0, 2, 1);
+  lfd.print_formatted(logfile, 0, false, 5, "~");
+  deallog << "LAPACKFullMatrix<double>::fill  SparseMatrixEZ<double>"
+         << std::endl;
+  lfd.fill(ez, 1, 1, 4, 2);
+  lfd.print_formatted(logfile, 0, false, 5, "~");
+}
diff --git a/tests/results/i686-pc-linux-gnu+gcc3.2/lac/matrices.output b/tests/results/i686-pc-linux-gnu+gcc3.2/lac/matrices.output
new file mode 100644 (file)
index 0000000..1301fd7
--- /dev/null
@@ -0,0 +1,19 @@
+
+DEAL::FullMatrix<float>::copy_from  SparseMatrixEZ<double>
+2.    ~     3.    4.    
+5.    6.    ~     7.    
+8.    9.    10.   11.   
+~     ~     ~     ~     
+12.   ~     13.   14.   
+DEAL::LAPACKFullMatrix<double>::copy_from  SparseMatrixEZ<double>
+2.    ~     3.    4.    
+5.    6.    ~     7.    
+8.    9.    10.   11.   
+~     ~     ~     ~     
+12.   ~     13.   14.   
+DEAL::LAPACKFullMatrix<double>::fill  SparseMatrixEZ<double>
+9.    10.   11.   
+~     ~     ~     
+DEAL::LAPACKFullMatrix<double>::fill  SparseMatrixEZ<double>
+9.    10.   11.   
+~     13.   14.   

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.