Changeset dfed1c for test/unit_test/periodic_fas.cpp
- Timestamp:
- Nov 22, 2011, 9:22:10 PM (14 years ago)
- Children:
- facba0
- Parents:
- 66f24d
- File:
-
- 1 edited
-
test/unit_test/periodic_fas.cpp (modified) (4 diffs)
Legend:
- Unmodified
- Added
- Removed
-
test/unit_test/periodic_fas.cpp
r66f24d rdfed1c 22 22 #endif 23 23 #include "solver/givens.hpp" 24 #include "solver/solver_ periodic.hpp"24 #include "solver/solver_singular.hpp" 25 25 #include "mg.hpp" 26 26 … … 37 37 Factory& factory = MG::GetFactory(); 38 38 39 Comm *comm = new CommSerial( Periodic);39 Comm *comm = new CommSerial(Boundary(Periodic, Periodic, Periodic)); 40 40 comm->Register("COMM"); 41 41 42 Interface* interface = new VMGInterfaces::InterfaceSinus( Periodic, 2, 6, 0.0, 1.0);42 Interface* interface = new VMGInterfaces::InterfaceSinus(comm->BoundaryConditions(), 2, 6, 0.0, 1.0); 43 43 MG::SetInterface(interface, comm); 44 44 … … 53 53 54 54 #ifdef HAVE_LAPACK 55 Solver* solver = new DSYSV<Solver Periodic>();55 Solver* solver = new DSYSV<SolverSingular>(); 56 56 #else 57 Solver* solver = new Givens<Solver Periodic>();57 Solver* solver = new Givens<SolverSingular>(); 58 58 #endif 59 59 solver->Register("SOLVER"); … … 78 78 MG::Solve(); 79 79 80 double res_init = MG::GetFactory().Get Object("INITIAL_RESIDUAL")->Cast< ObjectStorage<double> >()->Val();80 double res_init = MG::GetFactory().Get("INITIAL_RESIDUAL")->Cast< ObjectStorage<double> >()->Val(); 81 81 double res = MG::GetComm()->ComputeResidualNorm(*MG::GetSol(), *MG::GetRhs()); 82 82
Note:
See TracChangeset
for help on using the changeset viewer.
