Changeset e85cfd for src/units/particle/interface_fcs.cpp

Timestamp:

Apr 9, 2013, 9:45:13 AM (13 years ago)

Author:

Julian Iseringhausen <isering@…>

Children:

0bd47e

Parents:

4a709e

Message:

Work.

File:

• src/units/particle/interface_fcs.cpp (modified) (4 diffs)

src/units/particle/interface_fcs.cpp

@@ -63 +63 @@
 
 #include "cycles/cycle_cs_periodic_debug.hpp"
-
+//TODO: Remove Debug
 
 using namespace VMG;
@@ -112 +112 @@
                           periodic[2] ? Periodic : Open);
 
-  const bool singular = periodic[0] * periodic[1] * periodic[2];
+  const bool singular = boundary[0] == Periodic &&
+                        boundary[1] == Periodic &&
+                        boundary[2] == Periodic;
 
   /*
@@ -121 +123 @@
     new Particle::CommMPI(boundary, new DomainDecompositionMPI(), mpi_comm);
     new DiscretizationPoissonFD(discretization_order);
-    new InterfaceParticles(boundary, 2, level, Vector(box_offset), box_size, near_field_cells, 0, 1.0);
+    new InterfaceParticles(boundary, 2, level, Vector(box_offset), box_size, near_field_cells);
     new LevelOperatorCS(Stencils::RestrictionFullWeight, Stencils::InterpolationTrilinear);
     new Givens<SolverSingular>();
-    new CycleCSPeriodicDebug(cycle_type);
+    new CycleCSPeriodic(cycle_type);
 
   }else {
@@ -130 +132 @@
     new Particle::CommMPI(boundary, new DomainDecompositionMPI(), mpi_comm);
     new DiscretizationPoissonFV(discretization_order);
-    new InterfaceParticles(boundary, 2, level, Vector(box_offset), box_size, near_field_cells, 2, 1.6);
+    new InterfaceParticles(boundary, 2, level, Vector(box_offset), box_size, near_field_cells, 9, 1.6);
     new LevelOperatorFAS(Stencils::RestrictionFullWeight, Stencils::Injection, Stencils::InterpolationTrilinear);
     new Givens<SolverRegular>();