]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Test lac/SchurMatrix; Indented 165/head
authordakshinai <akdaks08@gmail.com>
Fri, 19 Sep 2014 23:40:55 +0000 (18:40 -0500)
committerdakshinai <akdaks08@gmail.com>
Fri, 19 Sep 2014 23:40:55 +0000 (18:40 -0500)
tests/lac/schur_matrix_01.cc
tests/lac/schur_matrix_02.cc
tests/lac/schur_matrix_03.cc
tests/lac/schur_matrix_04.cc

index 00dfccd9e67c3497156422288b2099353886874f..06524318909d1ab862df37896d553202a99dd26b 100644 (file)
 #include <vector>
 
 template<class MA_inverse, class MB, class MDt, class MC>
-  void
-  checkVmult(BlockVector<double> &dst, const BlockVector<double> &src,
-      const MA_inverse &Ainv, const MB &B, const MDt &Dt, const MC &C)
-  {
-    deallog << "vmult" << std::endl;
+void
+checkVmult(BlockVector<double> &dst, const BlockVector<double> &src,
+           const MA_inverse &Ainv, const MB &B, const MDt &Dt, const MC &C)
+{
+  deallog << "vmult" << std::endl;
 
-    GrowingVectorMemory < BlockVector<double> > mem;
+  GrowingVectorMemory < BlockVector<double> > mem;
 
-    SchurMatrix<MA_inverse, MB, MDt, MC> S(Ainv, B, Dt, C, mem);
+  SchurMatrix<MA_inverse, MB, MDt, MC> S(Ainv, B, Dt, C, mem);
 
-    S.debug_level(1);
+  S.debug_level(1);
 
-    S.vmult(dst, src);
+  S.vmult(dst, src);
 
-    for (unsigned int i = 0; i < dst.size(); ++i)
-      deallog << dst(i) << '\t';
-    deallog << std::endl;
-  }
+  for (unsigned int i = 0; i < dst.size(); ++i)
+    deallog << dst(i) << '\t';
+  deallog << std::endl;
+}
 
 int
 main()
@@ -59,38 +59,38 @@ main()
   ivector[0] = 1;
 
   double array[] =
-    { 2 };
+  { 2 };
   std::list<double> l(&array[0], &array[1]);
 
   BlockVector<double> src(ivector, l.begin(), l.end());
   BlockVector<double> dst(ivector);
 
   const double Adata1[] =
-    { 1 };
+  { 1 };
 
   const double Adata2[] =
-    { 2 };
+  { 2 };
 
   const double Adata3[] =
-    { 4 };
+  { 4 };
 
   const double Adata4[] =
-    { 5 };
+  { 5 };
 
   const double Ddata1[] =
-    { 3 };
+  { 3 };
 
   const double Ddata2[] =
-    { 6 };
+  { 6 };
 
   const double Bdata1[] =
-    { 7 };
+  { 7 };
 
   const double Bdata2[] =
-    { 8 };
+  { 8 };
 
   const double Cdata[] =
-    { 9 };
+  { 9 };
 
   FullMatrix<double> MA1(1, 1);
   FullMatrix<double> MA2(1, 1);
@@ -128,6 +128,6 @@ main()
   Dt.enter(MD2, 0, 1);
 
   checkVmult<BlockMatrixArray<double>, BlockMatrixArray<double>,
-      BlockMatrixArray<double>, BlockMatrixArray<double> >(dst, src, Ainv, B,
-      Dt, C);
+             BlockMatrixArray<double>, BlockMatrixArray<double> >(dst, src, Ainv, B,
+                 Dt, C);
 }
index def0f2715293f869019b5eda20b70f39430d0ca2..27cbd71872dc136e8f5e4930b4b8a739638b92cd 100644 (file)
 #include <vector>
 
 template<class MA_inverse, class MB, class MDt, class MC>
-  void
-  checkResidual(BlockVector<double> &dst, const BlockVector<double> &src,
-      const BlockVector<double> &rhs, const MA_inverse &Ainv, const MB &B,
-      const MDt &Dt, const MC &C)
-  {
-    deallog << "residual" << std::endl;
+void
+checkResidual(BlockVector<double> &dst, const BlockVector<double> &src,
+              const BlockVector<double> &rhs, const MA_inverse &Ainv, const MB &B,
+              const MDt &Dt, const MC &C)
+{
+  deallog << "residual" << std::endl;
 
-    GrowingVectorMemory < BlockVector<double> > mem;
+  GrowingVectorMemory < BlockVector<double> > mem;
 
-    SchurMatrix<MA_inverse, MB, MDt, MC> S(Ainv, B, Dt, C, mem);
+  SchurMatrix<MA_inverse, MB, MDt, MC> S(Ainv, B, Dt, C, mem);
 
-    double r = S.residual(dst, src, rhs);
+  double r = S.residual(dst, src, rhs);
 
-    deallog << "Residual:" << r << std::endl;
-  }
+  deallog << "Residual:" << r << std::endl;
+}
 
 int
 main()
@@ -56,43 +56,43 @@ main()
   ivector[0] = 1;
 
   double src_array[] =
-    { 2 };
+  { 2 };
   std::list<double> src_l(&src_array[0], &src_array[1]);
   BlockVector<double> src(ivector, src_l.begin(), src_l.end());
 
   double rhs_array[] =
-    { 1000 };
+  { 1000 };
   std::list<double> rhs_l(&rhs_array[0], &rhs_array[1]);
   BlockVector<double> rhs(ivector, rhs_l.begin(), rhs_l.end());
 
   BlockVector<double> dst(ivector);
 
   const double Adata1[] =
-    { 1 };
+  { 1 };
 
   const double Adata2[] =
-    { 2 };
+  { 2 };
 
   const double Adata3[] =
-    { 4 };
+  { 4 };
 
   const double Adata4[] =
-    { 5 };
+  { 5 };
 
   const double Ddata1[] =
-    { 3 };
+  { 3 };
 
   const double Ddata2[] =
-    { 6 };
+  { 6 };
 
   const double Bdata1[] =
-    { 7 };
+  { 7 };
 
   const double Bdata2[] =
-    { 8 };
+  { 8 };
 
   const double Cdata[] =
-    { 9 };
+  { 9 };
 
   FullMatrix<double> MA1(1, 1);
   FullMatrix<double> MA2(1, 1);
@@ -130,7 +130,7 @@ main()
   Dt.enter(MD2, 0, 1);
 
   checkResidual<BlockMatrixArray<double>, BlockMatrixArray<double>,
-      BlockMatrixArray<double>, BlockMatrixArray<double> >(dst, src, rhs, Ainv,
-      B, Dt, C);
+                BlockMatrixArray<double>, BlockMatrixArray<double> >(dst, src, rhs, Ainv,
+                    B, Dt, C);
 
 }
index e245493328fe05c823afc49bc7a08c61a82c4409..4056b50bee2c7daf2112d513e609e62804524413 100644 (file)
 #include <vector>
 
 template<class MA_inverse, class MB, class MDt, class MC>
-  void
-  checkPrepare_rhs(BlockVector<double> &dst, const BlockVector<double> &src,
-      const MA_inverse &Ainv, const MB &B, const MDt &Dt, const MC &C)
-  {
-    deallog << "prepare_rhs" << std::endl;
+void
+checkPrepare_rhs(BlockVector<double> &dst, const BlockVector<double> &src,
+                 const MA_inverse &Ainv, const MB &B, const MDt &Dt, const MC &C)
+{
+  deallog << "prepare_rhs" << std::endl;
 
-    GrowingVectorMemory < BlockVector<double> > mem;
+  GrowingVectorMemory < BlockVector<double> > mem;
 
-    SchurMatrix<MA_inverse, MB, MDt, MC> S(Ainv, B, Dt, C, mem);
+  SchurMatrix<MA_inverse, MB, MDt, MC> S(Ainv, B, Dt, C, mem);
 
-    S.debug_level(1);
+  S.debug_level(1);
 
-    S.prepare_rhs(dst, src);
+  S.prepare_rhs(dst, src);
 
-  }
+}
 
 int
 main()
@@ -57,7 +57,7 @@ main()
   ivector[1] = 1;
 
   double src_array[] =
-    { 9, 21 };
+  { 9, 21 };
   std::list<double> src_l(&src_array[0], &src_array[2]);
   BlockVector<double> src(ivector, src_l.begin(), src_l.end());
 
@@ -65,36 +65,36 @@ main()
   ivector2[0] = 1;
 
   double dst_array[] =
-    { 2 };
+  { 2 };
   std::list<double> dst_l(&dst_array[0], &dst_array[1]);
   BlockVector<double> dst(ivector2, dst_l.begin(), dst_l.end());
 
   const double Adata1[] =
-    { 1 };
+  { 1 };
 
   const double Adata2[] =
-    { 2 };
+  { 2 };
 
   const double Adata3[] =
-    { 4 };
+  { 4 };
 
   const double Adata4[] =
-    { 5 };
+  { 5 };
 
   const double Ddata1[] =
-    { 3 };
+  { 3 };
 
   const double Ddata2[] =
-    { 6 };
+  { 6 };
 
   const double Bdata1[] =
-    { 7 };
+  { 7 };
 
   const double Bdata2[] =
-    { 8 };
+  { 8 };
 
   const double Cdata[] =
-    { 9 };
+  { 9 };
 
   FullMatrix<double> MA1(1, 1);
   FullMatrix<double> MA2(1, 1);
@@ -132,7 +132,7 @@ main()
   Dt.enter(MD2, 0, 1);
 
   checkPrepare_rhs<BlockMatrixArray<double>, BlockMatrixArray<double>,
-      BlockMatrixArray<double>, BlockMatrixArray<double> >(dst, src, Ainv, B,
-      Dt, C);
+                   BlockMatrixArray<double>, BlockMatrixArray<double> >(dst, src, Ainv, B,
+                       Dt, C);
 
 }
index 0774c31ae247cc644e432b853e7d4fe4591148d4..010e376edec45d2f5da289d32d82722424a9797f 100644 (file)
 #include <vector>
 
 template<class MA_inverse, class MB, class MDt, class MC>
-  void
-  checkPostprocess(BlockVector<double> &dst, const BlockVector<double> &src,
-      const BlockVector<double> &rhs, const MA_inverse &Ainv, const MB &B,
-      const MDt &Dt, const MC &C)
-  {
-    deallog << "postprocess" << std::endl;
+void
+checkPostprocess(BlockVector<double> &dst, const BlockVector<double> &src,
+                 const BlockVector<double> &rhs, const MA_inverse &Ainv, const MB &B,
+                 const MDt &Dt, const MC &C)
+{
+  deallog << "postprocess" << std::endl;
 
-    GrowingVectorMemory < BlockVector<double> > mem;
+  GrowingVectorMemory < BlockVector<double> > mem;
 
-    SchurMatrix<MA_inverse, MB, MDt, MC> S(Ainv, B, Dt, C, mem);
+  SchurMatrix<MA_inverse, MB, MDt, MC> S(Ainv, B, Dt, C, mem);
 
-    S.debug_level(1);
+  S.debug_level(1);
 
-    S.postprocess(dst, src, rhs);
-  }
+  S.postprocess(dst, src, rhs);
+}
 
 int
 main()
@@ -57,7 +57,7 @@ main()
   ivector[1] = 1;
 
   double rhs_array[] =
-    { 9, 21 };
+  { 9, 21 };
   std::list<double> rhs_l(&rhs_array[0], &rhs_array[2]);
   BlockVector<double> rhs(ivector, rhs_l.begin(), rhs_l.end());
 
@@ -65,38 +65,38 @@ main()
   ivector2[0] = 1;
 
   double src_array[] =
-    { 2 };
+  { 2 };
   std::list<double> src_l(&src_array[0], &src_array[1]);
   BlockVector<double> src(ivector2, src_l.begin(), src_l.end());
 
   BlockVector<double> dst(ivector);
 
   const double Adata1[] =
-    { 1 };
+  { 1 };
 
   const double Adata2[] =
-    { 2 };
+  { 2 };
 
   const double Adata3[] =
-    { 4 };
+  { 4 };
 
   const double Adata4[] =
-    { 5 };
+  { 5 };
 
   const double Ddata1[] =
-    { 3 };
+  { 3 };
 
   const double Ddata2[] =
-    { 6 };
+  { 6 };
 
   const double Bdata1[] =
-    { 7 };
+  { 7 };
 
   const double Bdata2[] =
-    { 8 };
+  { 8 };
 
   const double Cdata[] =
-    { 9 };
+  { 9 };
 
   FullMatrix<double> MA1(1, 1);
   FullMatrix<double> MA2(1, 1);
@@ -134,7 +134,7 @@ main()
   Dt.enter(MD2, 0, 1);
 
   checkPostprocess<BlockMatrixArray<double>, BlockMatrixArray<double>,
-      BlockMatrixArray<double>, BlockMatrixArray<double> >(dst, src, rhs, Ainv,
-      B, Dt, C);
+                   BlockMatrixArray<double>, BlockMatrixArray<double> >(dst, src, rhs, Ainv,
+                       B, Dt, C);
 
 }

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.