]> https://gitweb.dealii.org/ - dealii.git/commitdiff
extend gla/vec_03 to complex-valued PETSc
authorDenis Davydov <davydden@gmail.com>
Tue, 29 Mar 2016 05:17:38 +0000 (07:17 +0200)
committerDenis Davydov <davydden@gmail.com>
Tue, 29 Mar 2016 18:48:46 +0000 (20:48 +0200)
tests/gla/vec_03.cc

index 2dbbe8506e92a9e3b59102151cafc8d197b36671..f33b0a0c085eb2e87f6ac01ad4acf64035b8ca54 100644 (file)
@@ -66,12 +66,12 @@ void test ()
 
   // check local values
   {
-    deallog << myid*2 << ": " << v(myid*2) << std::endl;
-    deallog << myid*2+1 << ": " << v(myid*2+1) << std::endl;
+    deallog << myid*2 << ": " << get_real_assert_zero_imag(v(myid*2)) << std::endl;
+    deallog << myid*2+1 << ": " << get_real_assert_zero_imag(v(myid*2+1)) << std::endl;
   }
 
-  Assert(v(myid*2) == myid*4.0, ExcInternalError());
-  Assert(v(myid*2+1) == myid*4.0+2.0, ExcInternalError());
+  Assert(get_real_assert_zero_imag(v(myid*2)) == myid*4.0, ExcInternalError());
+  Assert(get_real_assert_zero_imag(v(myid*2+1)) == myid*4.0+2.0, ExcInternalError());
 
   ConstraintMatrix cm;
   cm.add_line(1);
@@ -79,10 +79,10 @@ void test ()
   cm.close();
 
   if (myid==0)
-    deallog << "before: " << vb(1) << std::endl;
+    deallog << "before: " << get_real_assert_zero_imag(vb(1)) << std::endl;
   cm.distribute(vb); // this should set x(1)= 3.0 * x(2) = 12.0
   if (myid==0)
-    deallog << "after: " << vb(1) << std::endl;
+    deallog << "after: " << get_real_assert_zero_imag(vb(1)) << std::endl;
 
   // done
   if (myid==0)

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.