Changeset 728149 for src/interface/interface_particles_cf.cpp
- Timestamp:
- Apr 30, 2012, 2:40:31 PM (14 years ago)
- Children:
- 759a6a
- Parents:
- d448a7
- File:
-
- 1 edited
-
src/interface/interface_particles_cf.cpp (modified) (4 diffs)
Legend:
- Unmodified
- Added
- Removed
-
src/interface/interface_particles_cf.cpp
rd448a7 r728149 59 59 int max_level, max_iterations; 60 60 int pre_smoothing_steps, post_smoothing_steps; 61 int gamma, near_field_cells, discretization_order;61 int cycle_type, near_field_cells, discretization_order; 62 62 vmg_float precision; 63 63 std::string lop_str, datafile_str; … … 107 107 pre_smoothing_steps = Helper::ToValWithDefault<int>(xml_conf.child_value("pre_smoothing_steps"), 3); 108 108 post_smoothing_steps = Helper::ToValWithDefault<int>(xml_conf.child_value("post_smoothing_steps"), 3); 109 gamma = Helper::ToValWithDefault<int>(xml_conf.child_value("gamma"), 1);109 cycle_type = Helper::ToValWithDefault<int>(xml_conf.child_value("cycle_type"), 1); 110 110 near_field_cells = Helper::ToValWithDefault<int>(xml_conf.child_value("near_field_cells"), 3); 111 111 discretization_order = Helper::ToValWithDefault<int>(xml_conf.child_value("discretization_order"), 2); … … 133 133 if (!lop_str.compare("cs")) { 134 134 lop = new LevelOperatorCS(Stencils::RestrictionFullWeight, Stencils::InterpolationTrilinear); 135 Techniques::SetCorrectionSchemePeriodic(interface->MinLevel(), interface->MaxLevel(), gamma);135 Techniques::SetCorrectionSchemePeriodic(interface->MinLevel(), interface->MaxLevel(), cycle_type); 136 136 } else if (!lop_str.compare("fas")) { 137 137 lop = new LevelOperatorFAS(Stencils::RestrictionFullWeight, Stencils::Injection, Stencils::InterpolationTrilinear); 138 Techniques::SetFullApproximationSchemePeriodic(interface->MinLevel(), interface->MaxLevel(), gamma);138 Techniques::SetFullApproximationSchemePeriodic(interface->MinLevel(), interface->MaxLevel(), cycle_type); 139 139 } else if (!lop_str.compare("cs_debug")) { 140 140 lop = new LevelOperatorCS(Stencils::RestrictionFullWeight, Stencils::InterpolationTrilinear); 141 Techniques::SetCorrectionSchemePeriodicDebug(interface->MinLevel(), interface->MaxLevel(), gamma);141 Techniques::SetCorrectionSchemePeriodicDebug(interface->MinLevel(), interface->MaxLevel(), cycle_type); 142 142 } else if (!lop_str.compare("fas_debug")) { 143 143 lop = new LevelOperatorFAS(Stencils::RestrictionFullWeight, Stencils::Injection, Stencils::InterpolationTrilinear); 144 Techniques::SetFullApproximationSchemePeriodicDebug(interface->MinLevel(), interface->MaxLevel(), gamma);144 Techniques::SetFullApproximationSchemePeriodicDebug(interface->MinLevel(), interface->MaxLevel(), cycle_type); 145 145 } else { 146 146 lop = new LevelOperatorCS(Stencils::RestrictionFullWeight, Stencils::InterpolationTrilinear); 147 Techniques::SetCorrectionSchemePeriodic(interface->MinLevel(), interface->MaxLevel(), gamma);147 Techniques::SetCorrectionSchemePeriodic(interface->MinLevel(), interface->MaxLevel(), cycle_type); 148 148 } 149 149 lop->Register("LEVEL_OPERATOR"); … … 170 170 std::printf(" Maximum level: %d\n", max_level); 171 171 std::printf(" Maximum number of iterations: %d\n", max_iterations); 172 std::printf(" Gamma: %d\n", gamma);172 std::printf(" Gamma: %d\n", cycle_type); 173 173 std::printf(" Pre smoothing steps: %d\n", pre_smoothing_steps); 174 174 std::printf(" Post smoothing steps: %d\n", post_smoothing_steps);
Note:
See TracChangeset
for help on using the changeset viewer.
