diff --git a/apps/cpu/ConcreteExtrusion/ConcreteExtrusion.cpp b/apps/cpu/ConcreteExtrusion/ConcreteExtrusion.cpp
index cbbc9d8f92a2321677bd68034e332865acf4e335..1407d8c4d3f564f6bff351e7450199d983391622 100644
--- a/apps/cpu/ConcreteExtrusion/ConcreteExtrusion.cpp
+++ b/apps/cpu/ConcreteExtrusion/ConcreteExtrusion.cpp
@@ -53,7 +53,7 @@ void run(string configname)
real factorLz = config.getValue<real>("factorLz");
- SPtr<vf::mpi::Communicator> comm = vf::mpi::MPICommunicator::getInstance();
+ SPtr<vf::parallel::Communicator> comm = vf::parallel::MPICommunicator::getInstance();
int myid = comm->getProcessID();
if (myid == 0) UBLOG(logINFO, "Jet Breakup: Start!");
diff --git a/apps/cpu/JetBreakup/JetBreakup.cpp b/apps/cpu/JetBreakup/JetBreakup.cpp
index fecc1def8141e2835f1db8699e55bea8125944a5..f4a3bdc26ac6abc724aded591d691e77b2a50faa 100644
--- a/apps/cpu/JetBreakup/JetBreakup.cpp
+++ b/apps/cpu/JetBreakup/JetBreakup.cpp
@@ -52,10 +52,10 @@ void run(string configname)
real factorLy = config.getValue<real>("factorLy");
real factorLz = config.getValue<real>("factorLz");
- real dx = config.getValue<real>("dx");
+ //real dx = config.getValue<real>("dx");
- SPtr<vf::mpi::Communicator> comm = vf::mpi::MPICommunicator::getInstance();
+ SPtr<vf::parallel::Communicator> comm = vf::parallel::MPICommunicator::getInstance();
int myid = comm->getProcessID();
if (myid == 0) UBLOG(logINFO, "Jet Breakup: Start!");
@@ -141,7 +141,8 @@ void run(string configname)
real Re = rho_h * Uo * D / mu_h;
real We = rho_h * Uo * Uo * D / sigma;
- //real dx = Dg / D_LB;
+ real dx = Dg / D_LB;
+ //real D_LB = Dg / dx;
real nu_h = U_LB * D_LB / Re;
real nu_l = nu_h;
nu_h *= 0.1;
@@ -199,14 +200,31 @@ void run(string configname)
real g_minX2 = -0.5 * Ly;
real g_minX3 = -0.5 * Lz;
- // double g_maxX1 = 8.0*D;
- // double g_maxX2 = 2.5*D;
- // double g_maxX3 = 2.5*D;
-
real g_maxX1 = Lx;
real g_maxX2 = 0.5 * Ly;
real g_maxX3 = 0.5 * Lz;
+ //real g_minX1 = -0.5 * Lx;
+ //real g_minX2 = 0;
+ //real g_minX3 = -0.5 * Lz;
+
+ //real g_maxX1 = 0.5 * Lx;
+ //real g_maxX2 = Ly;
+ //real g_maxX3 = 0.5 * Lz;
+
+ //real g_minX1 = -0.5 * Lx;
+ //real g_minX2 = -0.5 * Ly;
+ //real g_minX3 = 0;
+
+ //real g_maxX1 = 0.5 * Lx;
+ //real g_maxX2 = 0.5 * Ly;
+ //real g_maxX3 = Lz;
+
+
+ // double g_maxX1 = 8.0*D;
+ // double g_maxX2 = 2.5*D;
+ // double g_maxX3 = 2.5*D;
+
SPtr<LBMUnitConverter> conv(new LBMUnitConverter());
// const int baseLevel = 0;
@@ -362,6 +380,10 @@ void run(string configname)
GbCylinder3DPtr geoInflow(new GbCylinder3D(g_minX1 - 2.0*dx, 0.0, 0.0, g_minX1, 0.0, 0.0, Dg / 2.0));
+ //GbCylinder3DPtr geoInflow(new GbCylinder3D(0.0, g_minX2 - 2.0 * dx, 0.0, 0.0, g_minX2, 0.0, Dg / 2.0));
+
+ //GbCylinder3DPtr geoInflow(new GbCylinder3D(0.0, 0.0, g_minX3 - 2.0 * dx, 0.0, 0.0, g_minX3, Dg / 2.0));
+
//GbCylinder3DPtr geoSolid(new GbCylinder3D(g_minX1 - 2.0 * dx, g_maxX2 / 2.0, g_maxX3 / 2.0, g_minX1+2.0*dx, g_maxX2 / 2.0, g_maxX3 / 2.0, 1.5*D / 2.0));
diff --git a/apps/cpu/Nozzle/nozzleSinglePhase.cpp b/apps/cpu/Nozzle/nozzleSinglePhase.cpp
index f0ffbd1f402b31ebba866ac4911702106af0739f..fe70ce94516a60ed17b559b3b8470b77341f7aaa 100644
--- a/apps/cpu/Nozzle/nozzleSinglePhase.cpp
+++ b/apps/cpu/Nozzle/nozzleSinglePhase.cpp
@@ -18,7 +18,7 @@ int main(int argc, char *argv[])
try {
- std::shared_ptr<vf::mpi::Communicator> comm = vf::mpi::MPICommunicator::getInstance();
+ std::shared_ptr<vf::parallel::Communicator> comm = vf::parallel::MPICommunicator::getInstance();
int myid = comm->getProcessID();
// bounding box
diff --git a/apps/cpu/NozzleMultiphase/nozzleMultiPhase.cpp b/apps/cpu/NozzleMultiphase/nozzleMultiPhase.cpp
index 24e52c2fcf8d6abbcb5af8b5c0e4807bf9dd62ba..ced2724ee6e8dddd3102428696a8647fee322326 100644
--- a/apps/cpu/NozzleMultiphase/nozzleMultiPhase.cpp
+++ b/apps/cpu/NozzleMultiphase/nozzleMultiPhase.cpp
@@ -18,7 +18,7 @@ int main(int argc, char *argv[])
try {
- std::shared_ptr<vf::mpi::Communicator> comm = vf::mpi::MPICommunicator::getInstance();
+ std::shared_ptr<vf::parallel::Communicator> comm = vf::parallel::MPICommunicator::getInstance();
int myid = comm->getProcessID();
// bounding box
diff --git a/apps/cpu/ShotcreteJet/jet.cpp b/apps/cpu/ShotcreteJet/jet.cpp
index f06fdefcd2194cd994f397a824478fcf8ee8d76d..18a018f81c598575313336c4c95a977d2b934b55 100644
--- a/apps/cpu/ShotcreteJet/jet.cpp
+++ b/apps/cpu/ShotcreteJet/jet.cpp
@@ -41,7 +41,7 @@ int main(int argc, char *argv[])
double restartStep = config.getValue<int>("restartStep");
int endTime = config.getValue<int>("endTime");
- std::shared_ptr<vf::mpi::Communicator> comm = vf::mpi::MPICommunicator::getInstance();
+ std::shared_ptr<vf::parallel::Communicator> comm = vf::parallel::MPICommunicator::getInstance();
int myid = comm->getProcessID();
if (myid == 0) UBLOG(logINFO, "Jet Breakup: Start!");
diff --git a/src/cpu/LiggghtsCoupling/SimulationObserver/LiggghtsCouplingSimulationObserver.cpp b/src/cpu/LiggghtsCoupling/SimulationObserver/LiggghtsCouplingSimulationObserver.cpp
index 20fb31cbf6e14c58ad092f6b8bd32e266c937787..8d4c069c7d3f93ace398acc60b07a018577a3ec9 100644
--- a/src/cpu/LiggghtsCoupling/SimulationObserver/LiggghtsCouplingSimulationObserver.cpp
+++ b/src/cpu/LiggghtsCoupling/SimulationObserver/LiggghtsCouplingSimulationObserver.cpp
@@ -1,6 +1,6 @@
#include "LiggghtsCouplingSimulationObserver.h"
#include "GbSphere3D.h"
-#include "mpi/MPICommunicator.h"
+#include "parallel/MPICommunicator.h"
#include "SimulationObserver.h"
#include "LiggghtsCoupling/3rdParty/LiggghtsCouplingWrapper.h"
#include "Grid3D.h"
diff --git a/src/cpu/VirtualFluidsCore/SimulationObservers/ShearStressSimulationObserver.cpp b/src/cpu/VirtualFluidsCore/SimulationObservers/ShearStressSimulationObserver.cpp
index 4e2e0111b70bd2b257bfe11a3f7fee64655d5c12..8a2b90698ee7f48dce5eea8f738bd69651d8b059 100644
--- a/src/cpu/VirtualFluidsCore/SimulationObservers/ShearStressSimulationObserver.cpp
+++ b/src/cpu/VirtualFluidsCore/SimulationObservers/ShearStressSimulationObserver.cpp
@@ -30,7 +30,7 @@ ShearStressSimulationObserver::ShearStressSimulationObserver(SPtr<Grid3D> grid,
for (SPtr<Block3D> block : blockVector[level]) {
UbTupleInt3 nx = grid->getBlockNX();
SPtr<ShearStressValuesArray3D> shearStressValues = SPtr<ShearStressValuesArray3D>(
- new ShearStressValuesArray3D(14, val<1>(nx) + 1, val<2>(nx) + 1, val<3>(nx) + 1, c0o1));
+ new ShearStressValuesArray3D(14, val<1>(nx) + 1, val<2>(nx) + 1, val<3>(nx) + 1, vf::basics::constant::c0o1));
block->getKernel()->getDataSet()->setShearStressValues(shearStressValues);
}
}