Changeset b57b9b for test/unit_test/library/periodic_fas.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/periodic_fas.cpp (modified) (2 diffs)

test/unit_test/library/periodic_fas.cpp

@@ -34 +34 @@
 #include "base/math.hpp"
 #include "comm/comm_serial.hpp"
+#include "cycles/cycle_fas_periodic.hpp"
 #include "level/level_operator_fas.hpp"
 #include "samples/discretization_poisson_fd.hpp"
-#include "samples/techniques.hpp"
-#include "smoother/gs.hpp"
 #include "smoother/gsrb.hpp"
-#ifdef HAVE_LAPACK
-#include "solver/dsysv.hpp"
-#endif
 #include "solver/givens.hpp"
 #include "solver/solver_singular.hpp"
@@ -56 +52 @@
   LibraryPeriodicFASFixture()
   {
-    Factory& factory = MG::GetFactory();
+    const Boundary boundary(Periodic, Periodic, Periodic);
 
-    Comm *comm = new CommSerial(Boundary(Periodic, Periodic, Periodic));
-    comm->Register("COMM");
+    new CommSerial(boundary);
+    new VMGInterfaces::InterfaceSinus(sine_factor, boundary, 2, 6, 0.0, 1.0);
+    new DiscretizationPoissonFD(2);
+    new LevelOperatorFAS(Stencils::RestrictionFullWeight, Stencils::Injection, Stencils::InterpolationTrilinear);
+    new GaussSeidelRB();
+    new Givens<SolverSingular>();
+    new CycleFASPeriodic(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 LevelOperatorFAS(Stencils::RestrictionFullWeight, Stencils::Injection, Stencils::InterpolationTrilinear);
-    lop->Register("LEVEL_OPERATOR");
-
-    Smoother* smoother = new GaussSeidelRB();
-    smoother->Register("SMOOTHER");
-
-#ifdef HAVE_LAPACK
-    Solver* solver = new DSYSV<SolverSingular>();
-#else
-    Solver* solver = new Givens<SolverSingular>();
-#endif
-    solver->Register("SOLVER");
-
-    Techniques::SetFullApproximationSchemePeriodic(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();