diff --git a/apps/benchmarks/GranularGas/MESA_PD_GranularGas.cpp b/apps/benchmarks/GranularGas/MESA_PD_GranularGas.cpp
index 01c253995a8270049349620a116ba45b0d165381..01dab71ea3bd53df776639564753e9c499818115 100644
--- a/apps/benchmarks/GranularGas/MESA_PD_GranularGas.cpp
+++ b/apps/benchmarks/GranularGas/MESA_PD_GranularGas.cpp
@@ -125,7 +125,7 @@ int main( int argc, char ** argv )
for (auto& iBlk : *forest)
{
for (auto pt : grid_generator::SCGrid(iBlk.getAABB(),
- Vector3<real_t>(params.spacing) * real_c(0.5),
+ Vector3<real_t>(params.spacing) * real_c(0.5) + params.shift,
params.spacing))
{
WALBERLA_CHECK(iBlk.getAABB().contains(pt));
diff --git a/apps/benchmarks/GranularGas/MESA_PD_KernelBenchmark.cpp b/apps/benchmarks/GranularGas/MESA_PD_KernelBenchmark.cpp
index a6cbb67471086729e009754b398dc79e299ba674..669da955b8ac17cd29216f5208ab1b65ba0cf349 100644
--- a/apps/benchmarks/GranularGas/MESA_PD_KernelBenchmark.cpp
+++ b/apps/benchmarks/GranularGas/MESA_PD_KernelBenchmark.cpp
@@ -128,7 +128,7 @@ int main( int argc, char ** argv )
for (auto& iBlk : *forest)
{
for (auto pt : grid_generator::SCGrid(iBlk.getAABB(),
- Vector3<real_t>(params.spacing) * real_c(0.5),
+ Vector3<real_t>(params.spacing) * real_c(0.5) + params.shift,
params.spacing))
{
WALBERLA_CHECK(iBlk.getAABB().contains(pt));
diff --git a/apps/benchmarks/GranularGas/MESA_PD_KernelLoadBalancing.cpp b/apps/benchmarks/GranularGas/MESA_PD_KernelLoadBalancing.cpp
index 76d90bbe9a73a2a07de60f6d61aece503db0b815..f4bedea104f0d3ab14e4786d9c3b6dc5470b337d 100644
--- a/apps/benchmarks/GranularGas/MESA_PD_KernelLoadBalancing.cpp
+++ b/apps/benchmarks/GranularGas/MESA_PD_KernelLoadBalancing.cpp
@@ -207,7 +207,7 @@ int main( int argc, char ** argv )
for (auto& iBlk : *forest)
{
for (auto pt : grid_generator::SCGrid(iBlk.getAABB(),
- Vector3<real_t>(params.spacing) * real_c(0.5),
+ Vector3<real_t>(params.spacing) * real_c(0.5) + params.shift,
params.spacing))
{
WALBERLA_CHECK(iBlk.getAABB().contains(pt));
diff --git a/apps/benchmarks/GranularGas/MESA_PD_LoadBalancing.cpp b/apps/benchmarks/GranularGas/MESA_PD_LoadBalancing.cpp
index feec74e0b8716d1c3d8a2b3ca77910e008f42ef0..7e8ea524c44c7832bc167211675786551c264b21 100644
--- a/apps/benchmarks/GranularGas/MESA_PD_LoadBalancing.cpp
+++ b/apps/benchmarks/GranularGas/MESA_PD_LoadBalancing.cpp
@@ -203,7 +203,7 @@ int main( int argc, char ** argv )
for (auto& iBlk : *forest)
{
for (auto pt : grid_generator::SCGrid(iBlk.getAABB(),
- Vector3<real_t>(params.spacing) * real_c(0.5),
+ Vector3<real_t>(params.spacing) * real_c(0.5) + params.shift,
params.spacing))
{
WALBERLA_CHECK(iBlk.getAABB().contains(pt));
diff --git a/apps/benchmarks/GranularGas/PE_GranularGas.cpp b/apps/benchmarks/GranularGas/PE_GranularGas.cpp
index ece49a7452f0f086629f2688e7004ae9017c2825..be3fd973c4e5dcfacbee18e2ea6ae1ca893b793d 100644
--- a/apps/benchmarks/GranularGas/PE_GranularGas.cpp
+++ b/apps/benchmarks/GranularGas/PE_GranularGas.cpp
@@ -190,7 +190,7 @@ int main( int argc, char ** argv )
for (auto& currentBlock : *forest)
{
for (auto it = grid_generator::SCIterator(currentBlock.getAABB().getIntersection(generationDomain),
- Vector3<real_t>(params.spacing) * real_c(0.5),
+ Vector3<real_t>(params.spacing) * real_c(0.5) + params.shift,
params.spacing);
it != grid_generator::SCIterator();
++it)
diff --git a/apps/benchmarks/GranularGas/PE_LoadBalancing.cpp b/apps/benchmarks/GranularGas/PE_LoadBalancing.cpp
index 42bcac4ece8efb4df0bc09b8fd2267a4259bd92c..84e2e8f41f7fc4eb723c3780b1e9b0d24681bf5b 100644
--- a/apps/benchmarks/GranularGas/PE_LoadBalancing.cpp
+++ b/apps/benchmarks/GranularGas/PE_LoadBalancing.cpp
@@ -268,7 +268,7 @@ int main( int argc, char ** argv )
for (auto& currentBlock : *forest)
{
for (auto it = grid_generator::SCIterator(currentBlock.getAABB().getIntersection(generationDomain),
- Vector3<real_t>(params.spacing) * real_c(0.5),
+ Vector3<real_t>(params.spacing) * real_c(0.5) + params.shift,
params.spacing);
it != grid_generator::SCIterator();
++it)
diff --git a/apps/benchmarks/GranularGas/Parameters.cpp b/apps/benchmarks/GranularGas/Parameters.cpp
index 83a517db8af38a5effe1efa9ed900a7232b8d328..9b9e904ca5744ee472ee91fed5ea4b8f3a94de2a 100644
--- a/apps/benchmarks/GranularGas/Parameters.cpp
+++ b/apps/benchmarks/GranularGas/Parameters.cpp
@@ -42,6 +42,9 @@ void loadFromConfig(Parameters& params, const Config::BlockHandle& cfg)
params.spacing = cfg.getParameter<real_t>("spacing", real_t(1.0) );
WALBERLA_LOG_INFO_ON_ROOT("spacing: " << params.spacing);
+ params.shift = cfg.getParameter<Vec3>("shift", Vec3(real_t(0.1), real_t(0.1), real_t(0.1)) );
+ WALBERLA_LOG_INFO_ON_ROOT("shift: " << params.shift);
+
params.radius = cfg.getParameter<real_t>("radius", real_t(0.5) );
WALBERLA_LOG_INFO_ON_ROOT("radius: " << params.radius);
@@ -135,6 +138,7 @@ void saveToSQL(const Parameters& params,
realProperties["spacing"] = double_c(params.spacing);
+
realProperties["radius"] = double_c(params.radius);
diff --git a/apps/benchmarks/GranularGas/Parameters.h b/apps/benchmarks/GranularGas/Parameters.h
index de6152baadc959e794a01e323987620acbc8601e..263b1ff5a8732d8babc110ceb7268b613260edae 100644
--- a/apps/benchmarks/GranularGas/Parameters.h
+++ b/apps/benchmarks/GranularGas/Parameters.h
@@ -39,6 +39,7 @@ struct Parameters
std::string sorting = "none";
Vec3 normal = Vec3(real_t(1.0), real_t(1.0), real_t(1.0));
real_t spacing = real_t(1.0);
+ Vec3 shift = Vec3(real_t(0.1), real_t(0.1), real_t(0.1));
real_t radius = real_t(0.5);
bool bBarrier = false;
bool storeNodeTimings = false;
diff --git a/apps/benchmarks/GranularGas/generateConfig.py b/apps/benchmarks/GranularGas/generateConfig.py
index 5a1bd8f753ea7824816ab93a84d436bded2fdd61..195534574c85a9fee6efe988b8dab8b0f77a7fb1 100755
--- a/apps/benchmarks/GranularGas/generateConfig.py
+++ b/apps/benchmarks/GranularGas/generateConfig.py
@@ -7,6 +7,7 @@ cfg = Config()
cfg.addParameter("sorting", "std::string", '"none"')
cfg.addParameter("normal", "Vec3", "Vec3(real_t(1.0), real_t(1.0), real_t(1.0))")
cfg.addParameter("spacing", "real_t", "real_t(1.0)")
+cfg.addParameter("shift", "Vec3", "Vec3(real_t(0.1), real_t(0.1), real_t(0.1))")
cfg.addParameter("radius", "real_t", "real_t(0.5)")
cfg.addParameter("bBarrier", "bool", "false")
cfg.addParameter("storeNodeTimings", "bool", "false")