Changeset b57b9b for test/unit_test/library/dirichlet_cs.cpp

Timestamp:

Sep 20, 2012, 1:54:01 PM (13 years ago)

Author:

Julian Iseringhausen <isering@…>

Children:

c3acbb

Parents:

2d3854

Message:

vmg: Simplified API.

git-svn-id: ​https://svn.version.fz-juelich.de/scafacos/trunk@2034 5161e1c8-67bf-11de-9fd5-51895aff932f

File:

• test/unit_test/library/dirichlet_cs.cpp (modified) (3 diffs)

test/unit_test/library/dirichlet_cs.cpp

@@ -34 +34 @@
 #include "base/math.hpp"
 #include "base/vector.hpp"
+#include "cycles/cycle_cs_dirichlet.hpp"
 #include "comm/comm_serial.hpp"
 #include "level/level_operator_cs.hpp"
@@ -39 +40 @@
 #include "level/stencils.hpp"
 #include "samples/discretization_poisson_fd.hpp"
-#include "samples/techniques.hpp"
 #include "smoother/gsrb.hpp"
 #ifdef HAVE_LAPACK
@@ -58 +58 @@
   LibraryDirichletCSFixture()
   {
-    Factory& factory = MG::GetFactory();
+    const Boundary boundary(Dirichlet, Dirichlet, Dirichlet);
 
-    Comm *comm = new CommSerial(Boundary(Dirichlet, Dirichlet, Dirichlet));
-    comm->Register("COMM");
+    new CommSerial(boundary);
+    new VMGInterfaces::InterfaceSinus(sine_factor, boundary, 2, 6, 0.0, 1.0);
+    new DiscretizationPoissonFD(2);
+    new LevelOperatorCS(Stencils::RestrictionFullWeight, Stencils::InterpolationTrilinear);
+    new GaussSeidelRB();
+    new Givens<SolverRegular>();
+    new CycleCSDirichlet(2);
 
-    Interface* interface = new VMGInterfaces::InterfaceSinus(sine_factor, comm->BoundaryConditions(), 2, 6, 0.0, 1.0);
-    interface->Register("INTERFACE");
-
-    Discretization* discretization = new DiscretizationPoissonFD(2);
-    discretization->Register("DISCRETIZATION");
-
-    LevelOperator* lop = new LevelOperatorCS(Stencils::RestrictionFullWeight, Stencils::InterpolationTrilinear);
-    lop->Register("LEVEL_OPERATOR");
-
-    Smoother* smoother = new GaussSeidelRB();
-    smoother->Register("SMOOTHER");
-
-#ifdef HAVE_LAPACK
-    Solver* solver = new DGESV<SolverRegular>();
-#else
-    Solver* solver = new Givens<SolverRegular>();
-#endif
-    solver->Register("SOLVER");
-
-    Techniques::SetCorrectionSchemeDirichlet(interface->MinLevel(), interface->MaxLevel(), 2);
-
-    factory.RegisterObjectStorage("PRESMOOTHSTEPS", 3);
-    factory.RegisterObjectStorage("POSTSMOOTHSTEPS", 3);
-    factory.RegisterObjectStorage("PRECISION", 1e-10);
-    factory.RegisterObjectStorage("MAX_ITERATION", 7);
+    new ObjectStorage<int>("PRESMOOTHSTEPS", 3);
+    new ObjectStorage<int>("POSTSMOOTHSTEPS", 3);
+    new ObjectStorage<int>("MAX_ITERATION", 7);
+    new ObjectStorage<vmg_float>("PRECISION", 1.0e-10);
 
     MG::PostInit();