Changeset 4571da for src/interface/interface_particles_cf.cpp

Timestamp:

Apr 27, 2012, 11:34:57 PM (14 years ago)

Author:

Julian Iseringhausen <isering@…>

Children:

1a92cf

Parents:

b2154a3

Message:

vmg: Implement fourth-order discretization of the Poisson equation.

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

File:

• src/interface/interface_particles_cf.cpp (modified) (4 diffs)

src/interface/interface_particles_cf.cpp

@@ -59 +59 @@
   int max_level, max_iterations;
   int pre_smoothing_steps, post_smoothing_steps;
-  int gamma, near_field_cells;
+  int gamma, near_field_cells, discretization_order;
   vmg_float precision;
   std::string lop_str, datafile_str;
@@ -109 +109 @@
   gamma = Helper::ToValWithDefault<int>(xml_conf.child_value("gamma"), 1);
   near_field_cells = Helper::ToValWithDefault<int>(xml_conf.child_value("near_field_cells"), 3);
+  discretization_order = Helper::ToValWithDefault<int>(xml_conf.child_value("discretization_order"), 2);
   precision = Helper::ToValWithDefault<vmg_float>(xml_conf.child_value("precision"), 1.0e-7);
   lop_str = xml_conf.child_value("level_operator");
@@ -126 +127 @@
   MG::SetInterface(interface, comm);
 
-  Discretization* discretization = new DiscretizationPoissonFD();
+  Discretization* discretization = new DiscretizationPoissonFD(discretization_order);
   discretization->Register("DISCRETIZATION");
 
@@ -174 +175 @@
     std::printf("  Precision:                    %e\n", precision);
     std::printf("  Near field cells:             %d\n", near_field_cells);
+    std::printf("  Discretization order:         %d\n", discretization_order);
   }