Changeset b57b9b for test/unit_test/library/force_sinus.cpp
- Timestamp:
- Sep 20, 2012, 1:54:01 PM (13 years ago)
- Children:
- c3acbb
- Parents:
- 2d3854
- File:
-
- 1 edited
-
test/unit_test/library/force_sinus.cpp (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
test/unit_test/library/force_sinus.cpp
r2d3854 rb57b9b 41 41 #include <iostream> 42 42 43 #include "base/factory.hpp"44 43 #include "base/math.hpp" 45 44 #include "base/vector.hpp" 46 45 #include "comm/comm_mpi.hpp" 47 46 #include "comm/domain_decomposition_mpi.hpp" 47 #include "cycles/cycle_cs_periodic.hpp" 48 48 #include "level/level_operator_cs.hpp" 49 49 #include "level/level_operator.hpp" 50 50 #include "samples/discretization_poisson_fd.hpp" 51 #include "samples/techniques.hpp"52 51 #include "smoother/gsrb.hpp" 53 52 #include "solver/givens.hpp" … … 67 66 LibraryForceSinusFixture() 68 67 { 69 Factory& factory = MG::GetFactory();68 const Boundary boundary(Periodic, Periodic, Periodic); 70 69 71 Comm* comm = new CommMPI(Boundary(Periodic, Periodic, Periodic), new DomainDecompositionMPI()); 72 Interface* interface = new VMGInterfaces::InterfaceSinus(sine_factor, comm->BoundaryConditions(), 2, 6, 0.0, 1.0); 73 Discretization* discretization = new DiscretizationPoissonFD(2); 74 LevelOperator* lop = new LevelOperatorCS(Stencils::RestrictionFullWeight, Stencils::InterpolationTrilinear); 75 Smoother* smoother = new GaussSeidelRB(); 76 Solver* solver = new Givens<SolverSingular>(); 70 new CommMPI(boundary, new DomainDecompositionMPI()); 71 new VMGInterfaces::InterfaceSinus(sine_factor, boundary, 2, 6, 0.0, 1.0); 72 new DiscretizationPoissonFD(2); 73 new LevelOperatorCS(Stencils::RestrictionFullWeight, Stencils::InterpolationTrilinear); 74 new GaussSeidelRB(); 75 new Givens<SolverSingular>(); 76 new CycleCSPeriodic(2); 77 77 78 comm->Register("COMM"); 79 interface->Register("INTERFACE"); 80 discretization->Register("DISCRETIZATION"); 81 lop->Register("LEVEL_OPERATOR"); 82 smoother->Register("SMOOTHER"); 83 solver->Register("SOLVER"); 84 85 Techniques::SetCorrectionSchemePeriodic(interface->MinLevel(), interface->MaxLevel(), 2); 86 87 factory.RegisterObjectStorage("PRESMOOTHSTEPS", 3); 88 factory.RegisterObjectStorage("POSTSMOOTHSTEPS", 3); 89 factory.RegisterObjectStorage("PRECISION", 1e-10); 90 factory.RegisterObjectStorage("MAX_ITERATION", 7); 78 new ObjectStorage<vmg_int>("PRESMOOTHSTEPS", 3); 79 new ObjectStorage<vmg_int>("POSTSMOOTHSTEPS", 3); 80 new ObjectStorage<vmg_int>("MAX_ITERATION", 7); 81 new ObjectStorage<vmg_float>("PRECISION", 1e-10); 91 82 92 83 MG::PostInit();
Note:
See TracChangeset
for help on using the changeset viewer.
