diff --git a/.gitlab-ci.yml b/.gitlab-ci.yml
index eecf91647df09b93800b8fd9e796682a38ce8f2a..fae336c8130c604ab6e90c1ed06bf820d8b46674 100644
--- a/.gitlab-ci.yml
+++ b/.gitlab-ci.yml
@@ -119,6 +119,8 @@ gcc_9_rebuild:
msvc_16:
stage: build
+ when: manual
+
tags:
- win
- gpu
@@ -198,6 +200,7 @@ msvc_16_unit_tests:
- gpu
needs: ["msvc_16"]
+ when: manual
before_script:
- cd $CI_PROJECT_DIR/build
@@ -433,6 +436,8 @@ vf_to_phoenix:
stage: deploy
needs: ["gcc_9_python", "gcc_9_unit_tests", "gcc_9_python_bindings_test"]
+ when: manual
+
variables:
SSH_KEY: "$SSH_PRIVATE_KEY"
HOST: "phoenix.hlr.rz.tu-bs.de"
@@ -469,6 +474,8 @@ cppcheck_to_elladan:
stage: deploy
needs: ["cppcheck"]
+ when: manual
+
variables:
SSH_KEY: "$SSH_PRIVATE_KEY_ELLADAN"
HOST: "elladan.irmb.bau.tu-bs.de"
@@ -482,6 +489,8 @@ gcov_to_elladan:
stage: deploy
needs: ["gcov_gcc_9"]
+ when: manual
+
variables:
SSH_KEY: "$SSH_PRIVATE_KEY_ELLADAN"
HOST: "elladan.irmb.bau.tu-bs.de"
@@ -498,6 +507,8 @@ sonar-scanner:
refs:
- develop
+ when: manual
+
variables:
SONAR_HOST_URL: "http://finrod.irmb.bau.tu-bs.de/sonarqube"
diff --git a/CMake/3rd/boost.cmake b/CMake/3rd/boost.cmake
index 544ae2d97b3b8ef0277445cf252149986c8dfb3a..b1446128602ccba67695c7abb38c59264e395967 100644
--- a/CMake/3rd/boost.cmake
+++ b/CMake/3rd/boost.cmake
@@ -8,7 +8,9 @@ function(linkBoost components)
set(Boost_USE_STATIC_RUNTIME OFF)
else()
set(Boost_USE_STATIC_LIBS ON)
- set(Boost_USE_STATIC_RUNTIME ON)
+ if(WIN32)
+ set(Boost_USE_STATIC_RUNTIME ON)
+ endif()
endif()
set(Boost_USE_MULTITHREADED ON)
diff --git a/CMake/VirtualFluidsMacros.cmake b/CMake/VirtualFluidsMacros.cmake
index 69e56fa63d4f25adb84e539678055e7f46b62d3d..2d6e63256955f8b739b2c84e695214185471600f 100644
--- a/CMake/VirtualFluidsMacros.cmake
+++ b/CMake/VirtualFluidsMacros.cmake
@@ -112,6 +112,15 @@ function(vf_add_library)
else()
vf_get_library_name (library_name)
endif()
+
+ if(NOT DEFINED ARG_BUILDTYPE)
+ if(BUILD_SHARED_LIBS)
+ set(ARG_BUILDTYPE "shared")
+ else()
+ set(ARG_BUILDTYPE "static")
+ endif()
+ endif()
+
status("Configuring the target: ${library_name} (type=${ARG_BUILDTYPE})...")
diff --git a/CMakeLists.txt b/CMakeLists.txt
index 50073d8a0016c0345e9ae7ae3cfc5507afb6ad14..1fd7dfff17385b7079c0d285b4f326677f65b736 100644
--- a/CMakeLists.txt
+++ b/CMakeLists.txt
@@ -8,8 +8,9 @@
#################################################################################
# required cmake versions
# CMAKE 3.13: target_link_options
+# CMAKE 3.15: CMAKE_MSVC_RUNTIME_LIBRARY
#################################################################################
-cmake_minimum_required(VERSION 3.13..3.18 FATAL_ERROR)
+cmake_minimum_required(VERSION 3.15..3.18 FATAL_ERROR)
project(VirtualFluids CXX)
@@ -43,7 +44,14 @@ option(BUILD_VF_INCLUDE_WHAT_YOU_USE "Add IWYU to the targets" OFF)
option(BUILD_VF_CPPCHECK "Add cppcheck to the targets" OFF)
option(BUILD_VF_COVERAGE "Add the -coverage compiler flag." OFF)
-option(BUILD_SHARED_LIBS "" ON)
+option(BUILD_SHARED_LIBS "" OFF)
+
+# windows: use multi-threaded dynamically-linked runtime library
+if(BUILD_SHARED_LIBS)
+ set(CMAKE_MSVC_RUNTIME_LIBRARY "MultiThreaded$<$<CONFIG:Debug>:Debug>DLL")
+else()
+ set(CMAKE_MSVC_RUNTIME_LIBRARY "MultiThreaded$<$<CONFIG:Debug>:Debug>")
+endif()
option(USE_OPENMP "Include OpenMP support" ON)
diff --git a/apps/gpu/LBM/DrivenCavity/CMakeLists.txt b/apps/gpu/LBM/DrivenCavity/CMakeLists.txt
index 3c401265dc43f0900094638280c5229894138c64..8bb9ee450616604b8cc636a2e13477d76e6b43d4 100644
--- a/apps/gpu/LBM/DrivenCavity/CMakeLists.txt
+++ b/apps/gpu/LBM/DrivenCavity/CMakeLists.txt
@@ -2,8 +2,9 @@ PROJECT(DrivenCavity)
#LIST(APPEND CAB_COMPILER_ADDTIONAL_CXX_COMPILER_FLAGS "-DOMPI_SKIP_MPICXX" )
-vf_add_library(BUILDTYPE binary PRIVATE_LINK VirtualFluids_GPU basics GridGenerator FILES DrivenCavity.cpp )
+vf_add_library(BUILDTYPE binary PRIVATE_LINK basics VirtualFluids_GPU GridGenerator FILES DrivenCavity.cpp )
set_source_files_properties(DrivenCavity.cpp PROPERTIES LANGUAGE CUDA)
linkCUDA()
+linkMPI()
diff --git a/apps/gpu/LBM/DrivenCavity/DrivenCavity.cpp b/apps/gpu/LBM/DrivenCavity/DrivenCavity.cpp
index 0af28e6f6c6eed2994d72ceeb1194c399ea79d5a..da53c05db9eeed857e8b6cfed4002aa99e398994 100644
--- a/apps/gpu/LBM/DrivenCavity/DrivenCavity.cpp
+++ b/apps/gpu/LBM/DrivenCavity/DrivenCavity.cpp
@@ -9,6 +9,8 @@
#include <exception>
#include <memory>
+#include "mpi.h"
+
//////////////////////////////////////////////////////////////////////////
#include "Core/DataTypes.h"
@@ -122,7 +124,9 @@ void multipleLevel(const std::string& configPath)
Communicator* comm = Communicator::getInstanz();
SPtr<ConfigFileReader> configReader = ConfigFileReader::getNewInstance();
- SPtr<ConfigData> configData = configReader->readConfigFile(configPath);
+
+ std::cout << configPath << std::endl;
+ SPtr<ConfigData> configData = configReader->readConfigFile(configPath.c_str());
////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
@@ -350,11 +354,9 @@ int main( int argc, char* argv[])
targetPath = __FILE__;
-#ifdef _WIN32
- targetPath = targetPath.substr(0, targetPath.find_last_of('\\') + 1);
-#else
targetPath = targetPath.substr(0, targetPath.find_last_of('/') + 1);
-#endif
+
+
std::cout << targetPath << std::endl;
diff --git a/gpu.cmake b/gpu.cmake
index f4faec7286c94d1d96386c9e193c0c8acbf29f0b..d1d5415081aa131868c0bc2ee047e8a140eb9a78 100644
--- a/gpu.cmake
+++ b/gpu.cmake
@@ -61,9 +61,9 @@ IF (VF.BUILD_VF_GPU)
#add_subdirectory(targets/apps/LBM/BaselMultiGPU)
add_subdirectory(apps/gpu/LBM/DrivenCavity)
- add_subdirectory(apps/gpu/LBM/gridGeneratorTest)
- add_subdirectory(apps/gpu/LBM/TGV_3D)
- add_subdirectory(apps/gpu/LBM/TGV_3D_MultiGPU)
+ #add_subdirectory(apps/gpu/LBM/gridGeneratorTest)
+ #add_subdirectory(apps/gpu/LBM/TGV_3D)
+ #add_subdirectory(apps/gpu/LBM/TGV_3D_MultiGPU)
ELSE()
MESSAGE( STATUS "exclude Virtual Fluids GPU." )
ENDIF()
diff --git a/src/basics/CMakeLists.txt b/src/basics/CMakeLists.txt
index 0b643212ad9c5da74e711d72da085c1fe670704e..2eff85ba5c48ec475cd8690060520cc09b4a0cb6 100644
--- a/src/basics/CMakeLists.txt
+++ b/src/basics/CMakeLists.txt
@@ -1,7 +1,7 @@
include(Core/buildInfo.cmake)
-vf_add_library(BUILDTYPE static EXCLUDE buildInfo.in.cpp)
+vf_add_library(EXCLUDE buildInfo.in.cpp)
vf_get_library_name (library_name)
target_include_directories(${library_name} PUBLIC ${CMAKE_CURRENT_SOURCE_DIR}/Core)
diff --git a/src/basics/Core/Input/ConfigFileReader/ConfigFileReader.cpp b/src/basics/Core/Input/ConfigFileReader/ConfigFileReader.cpp
index dce44b5551d2aedcf0c0d70a21dbcd9e8b071c8d..8fdb9fe1bc7f5dfe10dbfb522a87a4a3187dd60b 100644
--- a/src/basics/Core/Input/ConfigFileReader/ConfigFileReader.cpp
+++ b/src/basics/Core/Input/ConfigFileReader/ConfigFileReader.cpp
@@ -6,20 +6,17 @@
#include <fstream>
#include <iostream>
-BASICS_EXPORT std::shared_ptr<ConfigFileReader> ConfigFileReader::getNewInstance()
+std::shared_ptr<ConfigFileReader> ConfigFileReader::getNewInstance()
{
return std::shared_ptr<ConfigFileReader>(new ConfigFileReader());
}
-ConfigFileReader::ConfigFileReader() = default;
-
-BASICS_EXPORT ConfigFileReader::~ConfigFileReader() = default;
-
-BASICS_EXPORT std::shared_ptr<ConfigData> ConfigFileReader::readConfigFile(const std::string &filePath) const
+std::shared_ptr<ConfigData> ConfigFileReader::readConfigFile(const char* filePath) const
{
+ std::cout << filePath << std::endl;
std::shared_ptr<ConfigDataImp> data = ConfigDataImp::getNewInstance();
std::ifstream stream;
- stream.open(filePath.c_str(), std::ios::in);
+ stream.open(filePath, std::ios::in);
if (stream.fail())
throw std::runtime_error("can not open config file!");
std::unique_ptr<input::Input> input = input::Input::makeInput(stream, "config");
diff --git a/src/basics/Core/Input/ConfigFileReader/ConfigFileReader.h b/src/basics/Core/Input/ConfigFileReader/ConfigFileReader.h
index 9d88dfce6a7337d8ffaa7c4ffe43a4fd63949914..77c93ebfa4ba8564188d8e4a5442963382cf91e3 100644
--- a/src/basics/Core/Input/ConfigFileReader/ConfigFileReader.h
+++ b/src/basics/Core/Input/ConfigFileReader/ConfigFileReader.h
@@ -1,9 +1,11 @@
#ifndef CONFIGFILEREADER_H
#define CONFIGFILEREADER_H
-#include "../Input.h"
#include <memory>
+#include <string>
+
+#include "basics_export.h"
class ConfigData;
@@ -11,11 +13,10 @@ class ConfigFileReader
{
public:
BASICS_EXPORT static std::shared_ptr<ConfigFileReader> getNewInstance();
- BASICS_EXPORT virtual ~ConfigFileReader();
- BASICS_EXPORT std::shared_ptr<ConfigData> readConfigFile(const std::string &filePath) const;
+ BASICS_EXPORT std::shared_ptr<ConfigData> readConfigFile(const char* filePath) const;
private:
- ConfigFileReader();
+ ConfigFileReader() = default;
};
#endif
diff --git a/src/basics/Core/StringUtilities/StringUtil.h b/src/basics/Core/StringUtilities/StringUtil.h
index cdf8dce290110848d5e6c50eb2fac35822dddf76..5d458f7e4f3362fe5c10dc9ab5252266762d26f5 100644
--- a/src/basics/Core/StringUtilities/StringUtil.h
+++ b/src/basics/Core/StringUtilities/StringUtil.h
@@ -36,6 +36,16 @@ public:
static BASICS_EXPORT bool endsWith(const std::string &input, const std::string &end);
+
+ template<class T>
+ static T fromString(const std::string& s)
+ {
+ std::istringstream stream (s);
+ T t;
+ stream >> t;
+ return t;
+ }
+
private:
StringUtil() = default;
;
diff --git a/src/basics/Core/buildInfo.in.cpp b/src/basics/Core/buildInfo.in.cpp
index 56f302208256ba86d7f814bfa6711f704aea3479..482f4592a41bfed0615858869ec4bb297764b0b1 100644
--- a/src/basics/Core/buildInfo.in.cpp
+++ b/src/basics/Core/buildInfo.in.cpp
@@ -1,5 +1,3 @@
-#include "buildInfo.h"
-
#include "basics_export.h"
namespace buildInfo
diff --git a/src/basics/basics/writer/WbWriterVtkXmlBinary.h b/src/basics/basics/writer/WbWriterVtkXmlBinary.h
index 393c6bb13b268805a3fc2ca0850dd75d3fdd5616..421148d90497e3628ed274439c0b2fd7636b7fd2 100644
--- a/src/basics/basics/writer/WbWriterVtkXmlBinary.h
+++ b/src/basics/basics/writer/WbWriterVtkXmlBinary.h
@@ -37,7 +37,9 @@
#include <basics/writer/WbWriter.h>
-class WbWriterVtkXmlBinary : public WbWriter
+#include "basics_export.h"
+
+class BASICS_EXPORT WbWriterVtkXmlBinary : public WbWriter
{
public:
static WbWriterVtkXmlBinary *getInstance()
diff --git a/src/basics/geometry3d/GbObject3D.h b/src/basics/geometry3d/GbObject3D.h
index f60e64c8f95dfde84a389c422bbc6b9377e68187..c562f42412c8ab75b0ebcd303a5785ab4939cdfb 100644
--- a/src/basics/geometry3d/GbObject3D.h
+++ b/src/basics/geometry3d/GbObject3D.h
@@ -49,6 +49,8 @@ class GbObject3DCreator;
#include <PointerDefinitions.h>
+#include "basics_export.h"
+
//////////////////////////////////////////////////////////////////////////
//!
//! \class GbObject3D
@@ -57,7 +59,7 @@ class GbObject3DCreator;
//!
//////////////////////////////////////////////////////////////////////////
-class GbObject3D : public ObObject
+class BASICS_EXPORT GbObject3D : public ObObject
{
public:
// abstract Methods
diff --git a/src/cpu/simulationconfig/CMakeLists.txt b/src/cpu/simulationconfig/CMakeLists.txt
index 0d6be59ca140e2ba2a1d15b7b62a1e32736bf22b..2e5e7f2eb9cba467e1ecbe3ec9cea13b3a3f322d 100644
--- a/src/cpu/simulationconfig/CMakeLists.txt
+++ b/src/cpu/simulationconfig/CMakeLists.txt
@@ -16,5 +16,5 @@ set_target_properties(simulationconfig PROPERTIES CXX_STANDARD 14 POSITION_INDEP
target_include_directories(simulationconfig PUBLIC ${VFBUILDER_INCLUDE_DIR})
target_link_libraries(simulationconfig PRIVATE VirtualFluidsCore basics muparser)
target_compile_definitions(simulationconfig PRIVATE VF_METIS VF_MPI)
-
+target_include_directories(simulationconfig PUBLIC ${CMAKE_BINARY_DIR})
diff --git a/src/gpu/GridGenerator/CMakeLists.txt b/src/gpu/GridGenerator/CMakeLists.txt
index aaabd6d5de07940194e8428f0d3c9a18a741ae33..46866937e2fe91e4d2eedd3cbbaf16a65ae8a6fc 100644
--- a/src/gpu/GridGenerator/CMakeLists.txt
+++ b/src/gpu/GridGenerator/CMakeLists.txt
@@ -1,9 +1,16 @@
-vf_add_library(BUILDTYPE shared PRIVATE_LINK basics)
+vf_add_library(PRIVATE_LINK basics)
vf_get_library_name(library_name)
linkCUDA()
linkOpenMP(${library_name})
-set_target_properties(${library_name} PROPERTIES CUDA_SEPARABLE_COMPILATION ON)
\ No newline at end of file
+
+set_target_properties(${library_name} PROPERTIES CUDA_SEPARABLE_COMPILATION ON)
+
+# according to linker error when building static libraries.
+# https://stackoverflow.com/questions/50033435/cmake-cuda-separate-compilation-static-lib-link-error-on-windows-but-not-on-ubun
+if (NOT BUILD_SHARED_LIBRARY)
+ set_target_properties(${library_name} PROPERTIES CUDA_RESOLVE_DEVICE_SYMBOLS ON)
+endif()
\ No newline at end of file
diff --git a/src/gpu/VirtualFluids_GPU/CMakeLists.txt b/src/gpu/VirtualFluids_GPU/CMakeLists.txt
index 6635323c8be77c36682bac560dc0491d08f3a54f..f533c37993a0eaf97e40b51502058dc97ae32c72 100644
--- a/src/gpu/VirtualFluids_GPU/CMakeLists.txt
+++ b/src/gpu/VirtualFluids_GPU/CMakeLists.txt
@@ -7,7 +7,7 @@ ELSE(MSVC)
ENDIF(MSVC)
-vf_add_library(BUILDTYPE shared PRIVATE_LINK ${libsToLink})
+vf_add_library(PRIVATE_LINK ${libsToLink})
linkMPI()
linkCUDA()
diff --git a/src/gpu/VirtualFluids_GPU/Calculation/Cp.cpp b/src/gpu/VirtualFluids_GPU/Calculation/Cp.cpp
index 4473a849b97ba59ccc178093c1d36084372ec98b..7ecaa362cfa4253017b021cacb549ef018683a2b 100644
--- a/src/gpu/VirtualFluids_GPU/Calculation/Cp.cpp
+++ b/src/gpu/VirtualFluids_GPU/Calculation/Cp.cpp
@@ -7,6 +7,9 @@
#include <fstream>
#include <sstream>
+#include "Core/StringUtilities/StringUtil.h"
+#include <cassert>
+
//#include <math.h>
//#include "LB.h"
diff --git a/src/gpu/VirtualFluids_GPU/Calculation/Cp.h b/src/gpu/VirtualFluids_GPU/Calculation/Cp.h
index 203566c5ae609797660f55b66aa54c1d624550d2..d11be83ce347c99f8e8a7e687a8016777c5eb3e8 100644
--- a/src/gpu/VirtualFluids_GPU/Calculation/Cp.h
+++ b/src/gpu/VirtualFluids_GPU/Calculation/Cp.h
@@ -3,7 +3,6 @@
#include "LBM/LB.h"
#include "GPU/GPU_Interface.h"
-#include "Utilities/StringUtil.hpp"
#include "Parameter/Parameter.h"
#include "GPU/CudaMemoryManager.h"
diff --git a/src/gpu/VirtualFluids_GPU/Calculation/DragLift.cpp b/src/gpu/VirtualFluids_GPU/Calculation/DragLift.cpp
index 219914cf784baa74c11397ebed20eda03266ac98..df153312b489ebbe6d1cb7057797a1a3663b583f 100644
--- a/src/gpu/VirtualFluids_GPU/Calculation/DragLift.cpp
+++ b/src/gpu/VirtualFluids_GPU/Calculation/DragLift.cpp
@@ -7,6 +7,8 @@
#include <fstream>
#include <sstream>
+#include "Core/StringUtilities/StringUtil.h"
+
//#include <math.h>
//#include "LB.h"
diff --git a/src/gpu/VirtualFluids_GPU/Calculation/DragLift.h b/src/gpu/VirtualFluids_GPU/Calculation/DragLift.h
index 31ec90fd42615dcde7a6f57931ee3d160649ba97..7c9dd586fa98434deddf91ad33d53e7a053005fd 100644
--- a/src/gpu/VirtualFluids_GPU/Calculation/DragLift.h
+++ b/src/gpu/VirtualFluids_GPU/Calculation/DragLift.h
@@ -3,7 +3,6 @@
#include "LBM/LB.h"
#include "GPU/GPU_Interface.h"
-#include "Utilities/StringUtil.hpp"
#include "Parameter/Parameter.h"
#include "GPU/CudaMemoryManager.h"
diff --git a/src/gpu/VirtualFluids_GPU/Calculation/ForceCalculations.cpp b/src/gpu/VirtualFluids_GPU/Calculation/ForceCalculations.cpp
index 77835d66658acfa029c876188264f02204134777..d6b7ca332fbd51293471a8a54e90ce44cef130d0 100644
--- a/src/gpu/VirtualFluids_GPU/Calculation/ForceCalculations.cpp
+++ b/src/gpu/VirtualFluids_GPU/Calculation/ForceCalculations.cpp
@@ -7,9 +7,9 @@
#include <stdio.h>
#include <fstream>
#include <sstream>
-#include "Utilities/StringUtil.hpp"
#include "GPU/CudaMemoryManager.h"
+#include "Core/StringUtilities/StringUtil.h"
//using namespace std;
//////////////////////////////////////////////////////////////////////////
diff --git a/src/gpu/VirtualFluids_GPU/Calculation/PlaneCalculations.cpp b/src/gpu/VirtualFluids_GPU/Calculation/PlaneCalculations.cpp
index b92ddaaa3ee3449f657d84c7b7c044a51873f00f..e8d1cee6017c43a1447b0ba671ae9c247f75b16d 100644
--- a/src/gpu/VirtualFluids_GPU/Calculation/PlaneCalculations.cpp
+++ b/src/gpu/VirtualFluids_GPU/Calculation/PlaneCalculations.cpp
@@ -9,6 +9,7 @@
#include <sstream>
//using namespace std;
//////////////////////////////////////////////////////////////////////////
+#include "Core/StringUtilities/StringUtil.h"
void setSizeOfPlane(Parameter* para, int lev, unsigned int z)
{
diff --git a/src/gpu/VirtualFluids_GPU/Calculation/PlaneCalculations.h b/src/gpu/VirtualFluids_GPU/Calculation/PlaneCalculations.h
index 67cb7e2fec293da49230abf0eceafb6588463ea1..a04e6b5c3bf403cd38282382f6db253831f30678 100644
--- a/src/gpu/VirtualFluids_GPU/Calculation/PlaneCalculations.h
+++ b/src/gpu/VirtualFluids_GPU/Calculation/PlaneCalculations.h
@@ -3,7 +3,6 @@
#include "Parameter/Parameter.h"
#include "GPU/CudaMemoryManager.h"
-#include "Utilities/StringUtil.hpp"
#include "basics/utilities/UbSystem.h"
#include <iostream>
diff --git a/src/gpu/VirtualFluids_GPU/GPU/Cumulant27.cu b/src/gpu/VirtualFluids_GPU/GPU/Cumulant27.cu
index 45259f9ac16210ecf698ee8a59471acf10892e90..bc86f1a90fb7a61c8679b261eae5176572d059c6 100644
--- a/src/gpu/VirtualFluids_GPU/GPU/Cumulant27.cu
+++ b/src/gpu/VirtualFluids_GPU/GPU/Cumulant27.cu
@@ -180,9 +180,9 @@ extern "C" __global__ void LB_Kernel_Kum_AA2016_Comp_Bulk_SP_27(real omega,
(mfbbc-mfbba)) / rho;
////////////////////////////////////////////////////////////////////////////////////
//the force be with you
- real fx = forces[0]/(pow(c2o1,level)); //zero;//0.0032653/(pow(two,level)); //0.000000005;//(two/1600000.0) / 120.0; //
- real fy = forces[1]/(pow(c2o1,level)); //zero;
- real fz = forces[2]/(pow(c2o1,level)); //zero;
+ real fx = forces[0] / (pow((double)c2o1, (double)level)); // zero;//0.0032653/(pow(two,level)); //0.000000005;//(two/1600000.0) / 120.0; //
+ real fy = forces[1] / (pow((double)c2o1, (double)level)); // zero;
+ real fz = forces[2] / (pow((double)c2o1, (double)level)); // zero;
vvx += fx;
vvy += fy;
vvz += fz;
@@ -5611,9 +5611,9 @@ extern "C" __global__ void LB_Kernel_Kum_New_Comp_SRT_SP_27(
(mfbbc - mfbba)) / rho;
////////////////////////////////////////////////////////////////////////////////////
//the force be with you
- real fx = forces[0] / (pow(c2o1, level));
- real fy = forces[1] / (pow(c2o1, level));
- real fz = forces[2] / (pow(c2o1, level));
+ real fx = forces[0] / (pow((double)c2o1, (double)level));
+ real fy = forces[1] / (pow((double)c2o1, (double)level));
+ real fz = forces[2] / (pow((double)c2o1, (double)level));
vvx += fx*c1o2;
vvy += fy*c1o2;
vvz += fz*c1o2;
diff --git a/src/gpu/VirtualFluids_GPU/GPU/Cumulant27chim.cu b/src/gpu/VirtualFluids_GPU/GPU/Cumulant27chim.cu
index 1f232b76522dc03b8641ae3066bc47a98af4a164..8bfdfb3b9421baa4290cd305542520b28aa5fd17 100644
--- a/src/gpu/VirtualFluids_GPU/GPU/Cumulant27chim.cu
+++ b/src/gpu/VirtualFluids_GPU/GPU/Cumulant27chim.cu
@@ -255,9 +255,9 @@ extern "C" __global__ void Cumulant_One_preconditioned_errorDiffusion_chim_Comp_
(mfbbc - mfbba)) / rho;
////////////////////////////////////////////////////////////////////////////////////
//the force be with you
- real fx = forces[0] / (pow(c2o1, level)); //zero;//0.0032653/(pow(two,level)); //0.000000005;//(two/1600000.0) / 120.0; //
- real fy = forces[1] / (pow(c2o1, level)); //zero;
- real fz = forces[2] / (pow(c2o1, level)); //zero;
+ real fx = forces[0] / (pow((double)c2o1, (double)level)); //zero;//0.0032653/(pow(two,level)); //0.000000005;//(two/1600000.0) / 120.0; //
+ real fy = forces[1] / (pow((double)c2o1, (double)level)); //zero;
+ real fz = forces[2] / (pow((double)c2o1, (double)level)); //zero;
vvx += fx*c1o2;
vvy += fy*c1o2;
vvz += fz*c1o2;
@@ -1170,9 +1170,9 @@ extern "C" __global__ void Cumulant_One_preconditioned_chim_Comp_SP_27(
(mfbbc - mfbba)) / rho;
////////////////////////////////////////////////////////////////////////////////////
//the force be with you
- real fx = forces[0] / (pow(c2o1, level)); //zero;//0.0032653/(pow(two,level)); //0.000000005;//(two/1600000.0) / 120.0; //
- real fy = forces[1] / (pow(c2o1, level)); //zero;
- real fz = forces[2] / (pow(c2o1, level)); //zero;
+ real fx = forces[0] / (pow((double)c2o1, (double)level)); //zero;//0.0032653/(pow(two,level)); //0.000000005;//(two/1600000.0) / 120.0; //
+ real fy = forces[1] / (pow((double)c2o1, (double)level)); //zero;
+ real fz = forces[2] / (pow((double)c2o1, (double)level)); //zero;
vvx += fx*c1o2;
vvy += fy*c1o2;
vvz += fz*c1o2;
@@ -1943,9 +1943,9 @@ extern "C" __global__ void Cumulant_One_chim_Comp_SP_27(
(mfbbc - mfbba)) * OOrho;
////////////////////////////////////////////////////////////////////////////////////
//the force be with you
- real fx = forces[0] / (pow(c2o1, level)); //zero;//0.0032653/(pow(two,level)); //0.000000005;//(two/1600000.0) / 120.0; //
- real fy = forces[1] / (pow(c2o1, level)); //zero;
- real fz = forces[2] / (pow(c2o1, level)); //zero;
+ real fx = forces[0] / (pow((double)c2o1, (double)level)); //zero;//0.0032653/(pow(two,level)); //0.000000005;//(two/1600000.0) / 120.0; //
+ real fy = forces[1] / (pow((double)c2o1, (double)level)); //zero;
+ real fz = forces[2] / (pow((double)c2o1, (double)level)); //zero;
vvx += fx*c1o2;
vvy += fy*c1o2;
vvz += fz*c1o2;
diff --git a/src/gpu/VirtualFluids_GPU/GPU/Cumulant_F3_27.cu b/src/gpu/VirtualFluids_GPU/GPU/Cumulant_F3_27.cu
index fab34f1dd10b155cf694eaa719ba82057a7b3eb1..71d839609ec400fac4addb22224f18a69a5a9035 100644
--- a/src/gpu/VirtualFluids_GPU/GPU/Cumulant_F3_27.cu
+++ b/src/gpu/VirtualFluids_GPU/GPU/Cumulant_F3_27.cu
@@ -30,7 +30,7 @@ extern "C" __global__ void LB_PostProcessor_F3_2018_Fehlberg(real omega,
{
////////////////////////////////////////////////////////////////////////////////
const unsigned x = threadIdx.x; // Globaler x-Index
- const unsigned y = blockIdx.x; // Globaler y-Index
+ const unsigned y = blockIdx.x; // Globaler y-Index
const unsigned z = blockIdx.y; // Globaler z-Index
const unsigned nx = blockDim.x;
@@ -214,9 +214,9 @@ extern "C" __global__ void LB_PostProcessor_F3_2018_Fehlberg(real omega,
(mfbbc - mfbba)) / rho;
////////////////////////////////////////////////////////////////////////////////////
//the force be with you
- real fx = forces[0] / (pow(c2o1, level)); //zero;//0.0032653/(pow(two,level)); //0.000000005;//(two/1600000.0) / 120.0; //
- real fy = forces[1] / (pow(c2o1, level)); //zero;
- real fz = forces[2] / (pow(c2o1, level)); //zero;
+ real fx = forces[0] / (pow((double)c2o1, (double)level)); //zero;//0.0032653/(pow(two,level)); //0.000000005;//(two/1600000.0) / 120.0; //
+ real fy = forces[1] / (pow((double)c2o1, (double)level)); // zero;
+ real fz = forces[2] / (pow((double)c2o1, (double)level)); // zero;
vvx += fx;
vvy += fy;
vvz += fz;
diff --git a/src/gpu/VirtualFluids_GPU/GPU/Particles.cu b/src/gpu/VirtualFluids_GPU/GPU/Particles.cu
index dfc8271ed9f79baa5cf1bee8928619c75362bed7..39c63e152c7a5ab88d8d416591a737d56a5d23b3 100644
--- a/src/gpu/VirtualFluids_GPU/GPU/Particles.cu
+++ b/src/gpu/VirtualFluids_GPU/GPU/Particles.cu
@@ -110,9 +110,9 @@ extern "C" __global__ void InitParticles( real* coordX,
}
- real localX = randArray[k] / (real)(pow(c2o1,level));
- real localY = randArray[k] / (real)(pow(c2o1,level));
- real localZ = randArray[k] / (real)(pow(c2o1,level));
+ real localX = randArray[k] / (real)(pow((double)c2o1, (double)level));
+ real localY = randArray[k] / (real)(pow((double)c2o1, (double)level));
+ real localZ = randArray[k] / (real)(pow((double)c2o1, (double)level));
real globalX = coordinateX + localX;
real globalY = coordinateY + localY;
@@ -878,9 +878,9 @@ extern "C" __global__ void MoveParticles( real* coordX,
real localY = coordParticleYlocal[kTimeStepOld];
real localZ = coordParticleZlocal[kTimeStepOld];
- x = (localX * (real)(pow(c2o1,level))) - c1o2; //-c1o4;
- y = (localY * (real)(pow(c2o1,level))) - c1o2; //-c1o4;
- z = (localZ * (real)(pow(c2o1,level))) - c1o2; //-c1o4;
+ x = (localX * (real)(pow((double)c2o1, (double)level))) - c1o2; //-c1o4;
+ y = (localY * (real)(pow((double)c2o1, (double)level))) - c1o2; //-c1o4;
+ z = (localZ * (real)(pow((double)c2o1, (double)level))) - c1o2; //-c1o4;
////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
press = d0 + x*dx + y*dy + z*dz + x*y*dxy + x*z*dxz + y*z*dyz + x*y*z*dxyz;
vx1 = (a0 + x*ax + y*ay + z*az + x*x*axx + y*y*ayy + z*z*azz + x*y*axy + x*z*axz + y*z*ayz + x*y*z*axyz);
@@ -961,9 +961,9 @@ extern "C" __global__ void MoveParticles( real* coordX,
y = tempY;
z = tempZ;
- localX = (x + c1o2) / (real)(pow(c2o1,level));
- localY = (y + c1o2) / (real)(pow(c2o1,level));
- localZ = (z + c1o2) / (real)(pow(c2o1,level));
+ localX = (x + c1o2) / (real)(pow((double)c2o1, (double)level));
+ localY = (y + c1o2) / (real)(pow((double)c2o1, (double)level));
+ localZ = (z + c1o2) / (real)(pow((double)c2o1, (double)level));
coordParticleXlocal[kTimeStep] = localX;
coordParticleYlocal[kTimeStep] = localY;
coordParticleZlocal[kTimeStep] = localZ;
@@ -1743,9 +1743,9 @@ extern "C" __global__ void MoveParticlesWithoutBCs( real* coordX,
real localY = coordParticleYlocal[kTimeStepOld];
real localZ = coordParticleZlocal[kTimeStepOld];
- x = (localX * (real)(pow(c2o1,level))) - c1o2; //-c1o4;
- y = (localY * (real)(pow(c2o1,level))) - c1o2; //-c1o4;
- z = (localZ * (real)(pow(c2o1,level))) - c1o2; //-c1o4;
+ x = (localX * (real)(pow((double)c2o1, (double)level))) - c1o2; //-c1o4;
+ y = (localY * (real)(pow((double)c2o1, (double)level))) - c1o2; //-c1o4;
+ z = (localZ * (real)(pow((double)c2o1, (double)level))) - c1o2; //-c1o4;
////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
press = d0 + x*dx + y*dy + z*dz + x*y*dxy + x*z*dxz + y*z*dyz + x*y*z*dxyz;
vx1 = (a0 + x*ax + y*ay + z*az + x*x*axx + y*y*ayy + z*z*azz + x*y*axy + x*z*axz + y*z*ayz + x*y*z*axyz);
@@ -1826,9 +1826,9 @@ extern "C" __global__ void MoveParticlesWithoutBCs( real* coordX,
y = tempY;
z = tempZ;
- localX = (x + c1o2) / (real)(pow(c2o1,level));
- localY = (y + c1o2) / (real)(pow(c2o1,level));
- localZ = (z + c1o2) / (real)(pow(c2o1,level));
+ localX = (x + c1o2) / (real)(pow((double)c2o1, (double)level));
+ localY = (y + c1o2) / (real)(pow((double)c2o1, (double)level));
+ localZ = (z + c1o2) / (real)(pow((double)c2o1, (double)level));
coordParticleXlocal[kTimeStep] = localX;
coordParticleYlocal[kTimeStep] = localY;
coordParticleZlocal[kTimeStep] = localZ;
diff --git a/src/gpu/VirtualFluids_GPU/GPU/WaleCumulant27.cu b/src/gpu/VirtualFluids_GPU/GPU/WaleCumulant27.cu
index f449e0d4923c0b8e1e472cbb2e70b562a7280d30..e4e6d9fb2b4582afb95b97826bedacd96bd80eb7 100644
--- a/src/gpu/VirtualFluids_GPU/GPU/WaleCumulant27.cu
+++ b/src/gpu/VirtualFluids_GPU/GPU/WaleCumulant27.cu
@@ -311,9 +311,9 @@ extern "C" __global__ void LB_Kernel_WaleBySoniMalav_Cum_AA2016_Comp_SP_27(
////////////////////////////////////////////////////////////////////////////////////
//the force be with you
{
- real fx = forces[0] / (pow(c2o1, level)); //zero;//0.0032653/(pow(two,level)); //0.000000005;//(two/1600000.0) / 120.0; //
- real fy = forces[1] / (pow(c2o1, level)); //zero;
- real fz = forces[2] / (pow(c2o1, level)); //zero;
+ real fx = forces[0] / (pow((double)c2o1, (double)level)); //zero;//0.0032653/(pow(two,level)); //0.000000005;//(two/1600000.0) / 120.0; //
+ real fy = forces[1] / (pow((double)c2o1, (double)level)); //zero;
+ real fz = forces[2] / (pow((double)c2o1, (double)level)); //zero;
vvx += fx*c1o2;
vvy += fy*c1o2;
vvz += fz*c1o2;
diff --git a/src/gpu/VirtualFluids_GPU/Input/ConfigFile.cpp b/src/gpu/VirtualFluids_GPU/Input/ConfigFile.cpp
index 1f6b69c6334e9a6eed9db78679008ca19c945946..613af04d7f89a646f5c8327e5f0fe697d043c6c8 100644
--- a/src/gpu/VirtualFluids_GPU/Input/ConfigFile.cpp
+++ b/src/gpu/VirtualFluids_GPU/Input/ConfigFile.cpp
@@ -1,11 +1,12 @@
#include "ConfigFile.h"
-#include "Utilities/StringUtil.hpp"
#include <errno.h>
#include <algorithm>
#include <sstream>
#include <iostream>
#include <string>
+#include "Core/StringUtilities/StringUtil.h"
+
ConfigFile::ConfigFile( const char *strConfigFile )
{
diff --git a/src/gpu/VirtualFluids_GPU/Input/VtkXmlReader.hpp b/src/gpu/VirtualFluids_GPU/Input/VtkXmlReader.hpp
index ce00ab7e8876ad73d6bfee410c92bc9b57e2266b..3c94a6a91de739da0056bf793a60f2dcede05b7e 100644
--- a/src/gpu/VirtualFluids_GPU/Input/VtkXmlReader.hpp
+++ b/src/gpu/VirtualFluids_GPU/Input/VtkXmlReader.hpp
@@ -10,7 +10,6 @@
#ifndef VTK_XML_READER_H
#define VTK_XML_READER_H
-#include "Utilities/StringUtil.hpp"
#include <boost/property_tree/ptree.hpp>
#include <boost/property_tree/xml_parser.hpp>
#include <boost/foreach.hpp>
@@ -18,6 +17,8 @@
#include <vector>
#include <iostream>
+#include "Core/StringUtilities/StringUtil.h"
+
namespace vtk_xml_reader
{
template <class T>
diff --git a/src/gpu/VirtualFluids_GPU/Kernel/Kernels/BasicKernels/Advection/Compressible/CumulantAll4/CumulantAll4CompSP27_Device.cu b/src/gpu/VirtualFluids_GPU/Kernel/Kernels/BasicKernels/Advection/Compressible/CumulantAll4/CumulantAll4CompSP27_Device.cu
index e17a5263565b2e134eb868b5b6271bfc0ff3590c..246ec79658674c8383012f97776ca60554883158 100644
--- a/src/gpu/VirtualFluids_GPU/Kernel/Kernels/BasicKernels/Advection/Compressible/CumulantAll4/CumulantAll4CompSP27_Device.cu
+++ b/src/gpu/VirtualFluids_GPU/Kernel/Kernels/BasicKernels/Advection/Compressible/CumulantAll4/CumulantAll4CompSP27_Device.cu
@@ -206,9 +206,9 @@ extern "C" __global__ void LB_Kernel_Cumulant_D3Q27All4( real omega,
(mfbbc - mfbba)) / rho;
////////////////////////////////////////////////////////////////////////////////////
//the force be with you
- real fx = forces[0] / (pow(c2o1, level)); //zero;//0.0032653/(pow(two,level)); //0.000000005;//(two/1600000.0) / 120.0; //
- real fy = forces[1] / (pow(c2o1, level)); //zero;
- real fz = forces[2] / (pow(c2o1, level)); //zero;
+ real fx = forces[0] / (pow((double)c2o1, (double)level)); //zero;//0.0032653/(pow(two,level)); //0.000000005;//(two/1600000.0) / 120.0; //
+ real fy = forces[1] / (pow((double)c2o1, (double)level)); //zero;
+ real fz = forces[2] / (pow((double)c2o1, (double)level)); //zero;
vvx += fx;
vvy += fy;
vvz += fz;
diff --git a/src/gpu/VirtualFluids_GPU/Kernel/Kernels/BasicKernels/Advection/Compressible/CumulantK15/CumulantK15Comp_Device.cu b/src/gpu/VirtualFluids_GPU/Kernel/Kernels/BasicKernels/Advection/Compressible/CumulantK15/CumulantK15Comp_Device.cu
index 183d83a6c89ca0470275302842fa267e8576cc81..4fb90b7093e9608cb5dce616e56f79d90f540067 100644
--- a/src/gpu/VirtualFluids_GPU/Kernel/Kernels/BasicKernels/Advection/Compressible/CumulantK15/CumulantK15Comp_Device.cu
+++ b/src/gpu/VirtualFluids_GPU/Kernel/Kernels/BasicKernels/Advection/Compressible/CumulantK15/CumulantK15Comp_Device.cu
@@ -202,9 +202,9 @@ extern "C" __global__ void LB_Kernel_CumulantK15Comp(real omega,
(mfbbc - mfbba)) / rho;
////////////////////////////////////////////////////////////////////////////////////
//the force be with you
- real fx = forces[0] / (pow(c2o1, level)); //zero;//0.0032653/(pow(two,level)); //0.000000005;//(two/1600000.0) / 120.0; //
- real fy = forces[1] / (pow(c2o1, level)); //zero;
- real fz = forces[2] / (pow(c2o1, level)); //zero;
+ real fx = forces[0] / (pow((double)c2o1, (double)level)); //zero;//0.0032653/(pow(two,level)); //0.000000005;//(two/1600000.0) / 120.0; //
+ real fy = forces[1] / (pow((double)c2o1, (double)level)); //zero;
+ real fz = forces[2] / (pow((double)c2o1, (double)level)); //zero;
vvx += fx*c1o2;
vvy += fy*c1o2;
vvz += fz*c1o2;
diff --git a/src/gpu/VirtualFluids_GPU/Kernel/Kernels/BasicKernels/Advection/Compressible/CumulantK15Bulk/CumulantK15BulkComp_Device.cu b/src/gpu/VirtualFluids_GPU/Kernel/Kernels/BasicKernels/Advection/Compressible/CumulantK15Bulk/CumulantK15BulkComp_Device.cu
index dfda290109dad41ed6822e2dbf95ba17264e9ae2..1a0efdf8e4e70633f4b0052763813e0639e387d3 100644
--- a/src/gpu/VirtualFluids_GPU/Kernel/Kernels/BasicKernels/Advection/Compressible/CumulantK15Bulk/CumulantK15BulkComp_Device.cu
+++ b/src/gpu/VirtualFluids_GPU/Kernel/Kernels/BasicKernels/Advection/Compressible/CumulantK15Bulk/CumulantK15BulkComp_Device.cu
@@ -202,9 +202,9 @@ extern "C" __global__ void LB_Kernel_CumulantK15BulkComp(real omega,
(mfbbc - mfbba)) / rho;
////////////////////////////////////////////////////////////////////////////////////
//the force be with you
- real fx = forces[0] / (pow(c2o1, level)); //zero;//0.0032653/(pow(two,level)); //0.000000005;//(two/1600000.0) / 120.0; //
- real fy = forces[1] / (pow(c2o1, level)); //zero;
- real fz = forces[2] / (pow(c2o1, level)); //zero;
+ real fx = forces[0] / (pow((double)c2o1, (double)level)); //zero;//0.0032653/(pow(two,level)); //0.000000005;//(two/1600000.0) / 120.0; //
+ real fy = forces[1] / (pow((double)c2o1, (double)level)); //zero;
+ real fz = forces[2] / (pow((double)c2o1, (double)level)); //zero;
vvx += fx;
vvy += fy;
vvz += fz;
diff --git a/src/gpu/VirtualFluids_GPU/Kernel/Kernels/BasicKernels/Advection/Compressible/CumulantK17/CumulantK17Comp_Device.cu b/src/gpu/VirtualFluids_GPU/Kernel/Kernels/BasicKernels/Advection/Compressible/CumulantK17/CumulantK17Comp_Device.cu
index 9d5f512bdef702973bea90cec587b7d29dfc1a2c..9ccd7dcb4ede3a493fe6c2621b03a515544bec96 100644
--- a/src/gpu/VirtualFluids_GPU/Kernel/Kernels/BasicKernels/Advection/Compressible/CumulantK17/CumulantK17Comp_Device.cu
+++ b/src/gpu/VirtualFluids_GPU/Kernel/Kernels/BasicKernels/Advection/Compressible/CumulantK17/CumulantK17Comp_Device.cu
@@ -204,9 +204,9 @@ extern "C" __global__ void LB_Kernel_CumulantK17Comp(real omega,
(mfbbc - mfbba)) / rho;
////////////////////////////////////////////////////////////////////////////////////
//the force be with you
- real fx = forces[0] / (pow(c2o1, level)); //zero;//0.0032653/(pow(two,level)); //0.000000005;//(two/1600000.0) / 120.0; //
- real fy = forces[1] / (pow(c2o1, level)); //zero;
- real fz = forces[2] / (pow(c2o1, level)); //zero;
+ real fx = forces[0] / (pow((double)c2o1, (double)level)); //zero;//0.0032653/(pow(two,level)); //0.000000005;//(two/1600000.0) / 120.0; //
+ real fy = forces[1] / (pow((double)c2o1, (double)level)); //zero;
+ real fz = forces[2] / (pow((double)c2o1, (double)level)); //zero;
vvx += fx*c1o2;
vvy += fy*c1o2;
vvz += fz*c1o2;
diff --git a/src/gpu/VirtualFluids_GPU/Kernel/Kernels/BasicKernels/Advection/Compressible/CumulantK17Bulk/CumulantK17BulkComp_Device.cu b/src/gpu/VirtualFluids_GPU/Kernel/Kernels/BasicKernels/Advection/Compressible/CumulantK17Bulk/CumulantK17BulkComp_Device.cu
index 653fd9f7bbbf98afa5e993d90ee5ee8a8ed5a71c..147b814a3c603637c8bb721e3fe96a7b6f2ce8ba 100644
--- a/src/gpu/VirtualFluids_GPU/Kernel/Kernels/BasicKernels/Advection/Compressible/CumulantK17Bulk/CumulantK17BulkComp_Device.cu
+++ b/src/gpu/VirtualFluids_GPU/Kernel/Kernels/BasicKernels/Advection/Compressible/CumulantK17Bulk/CumulantK17BulkComp_Device.cu
@@ -172,9 +172,9 @@ extern "C" __global__ void LB_Kernel_CumulantK17BulkComp(real omega,
(mfbbc - mfbba)) / rho;
////////////////////////////////////////////////////////////////////////////////////
//the force be with you
- real fx = forces[0] / (pow(c2o1, level)); //zero;//0.0032653/(pow(two,level)); //0.000000005;//(two/1600000.0) / 120.0; //
- real fy = forces[1] / (pow(c2o1, level)); //zero;
- real fz = forces[2] / (pow(c2o1, level)); //zero;
+ real fx = forces[0] / (pow((double)c2o1, (double)level)); //zero;//0.0032653/(pow(two,level)); //0.000000005;//(two/1600000.0) / 120.0; //
+ real fy = forces[1] / (pow((double)c2o1, (double)level)); //zero;
+ real fz = forces[2] / (pow((double)c2o1, (double)level)); //zero;
vvx += fx;
vvy += fy;
vvz += fz;
diff --git a/src/gpu/VirtualFluids_GPU/Kernel/Kernels/BasicKernels/Advection/Compressible/CumulantK18/CumulantK18Comp_Device.cu b/src/gpu/VirtualFluids_GPU/Kernel/Kernels/BasicKernels/Advection/Compressible/CumulantK18/CumulantK18Comp_Device.cu
index f882ecc2aa4eadb2c95b171e95c776877da28753..2249e5857cb8c5d28b53658c25bbb9c3fc3e62c0 100644
--- a/src/gpu/VirtualFluids_GPU/Kernel/Kernels/BasicKernels/Advection/Compressible/CumulantK18/CumulantK18Comp_Device.cu
+++ b/src/gpu/VirtualFluids_GPU/Kernel/Kernels/BasicKernels/Advection/Compressible/CumulantK18/CumulantK18Comp_Device.cu
@@ -204,9 +204,9 @@ extern "C" __global__ void LB_Kernel_CumulantK18Comp(
(mfbbc - mfbba)) / rho;
////////////////////////////////////////////////////////////////////////////////////
//the force be with you
- real fx = forces[0] / (pow(c2o1, level)); //zero;//0.0032653/(pow(two,level)); //0.000000005;//(two/1600000.0) / 120.0; //
- real fy = forces[1] / (pow(c2o1, level)); //zero;
- real fz = forces[2] / (pow(c2o1, level)); //zero;
+ real fx = forces[0] / (pow((double)c2o1, (double)level)); //zero;//0.0032653/(pow(two,level)); //0.000000005;//(two/1600000.0) / 120.0; //
+ real fy = forces[1] / (pow((double)c2o1, (double)level)); //zero;
+ real fz = forces[2] / (pow((double)c2o1, (double)level)); //zero;
vvx += fx;
vvy += fy;
vvz += fz;
diff --git a/src/gpu/VirtualFluids_GPU/Kernel/Kernels/BasicKernels/Advection/Compressible/CumulantK20/CumulantK20Comp_Device.cu b/src/gpu/VirtualFluids_GPU/Kernel/Kernels/BasicKernels/Advection/Compressible/CumulantK20/CumulantK20Comp_Device.cu
index 3324d52e9470d9c226c561297ee2a1347a743738..9edae07fabf8f2c18fb86055574b153e960d8dec 100644
--- a/src/gpu/VirtualFluids_GPU/Kernel/Kernels/BasicKernels/Advection/Compressible/CumulantK20/CumulantK20Comp_Device.cu
+++ b/src/gpu/VirtualFluids_GPU/Kernel/Kernels/BasicKernels/Advection/Compressible/CumulantK20/CumulantK20Comp_Device.cu
@@ -203,9 +203,9 @@ extern "C" __global__ void LB_Kernel_CumulantK20Comp(
(mfbbc - mfbba)) / rho;
////////////////////////////////////////////////////////////////////////////////////
//the force be with you
- real fx = forces[0] / (pow(c2o1, level)); //zero;//0.0032653/(pow(two,level)); //0.000000005;//(two/1600000.0) / 120.0; //
- real fy = forces[1] / (pow(c2o1, level)); //zero;
- real fz = forces[2] / (pow(c2o1, level)); //zero;
+ real fx = forces[0] / (pow((double)c2o1, (double)level)); //zero;//0.0032653/(pow(two,level)); //0.000000005;//(two/1600000.0) / 120.0; //
+ real fy = forces[1] / (pow((double)c2o1, (double)level)); //zero;
+ real fz = forces[2] / (pow((double)c2o1, (double)level)); //zero;
vvx += fx;
vvy += fy;
vvz += fz;
diff --git a/src/gpu/VirtualFluids_GPU/Kernel/Kernels/PorousMediaKernels/Advection/Compressible/CumulantOne/PMCumulantOneCompSP27_Device.cu b/src/gpu/VirtualFluids_GPU/Kernel/Kernels/PorousMediaKernels/Advection/Compressible/CumulantOne/PMCumulantOneCompSP27_Device.cu
index 2b6c98b388facf63bbbe9a9ede746269a3116fb7..945d8706415cfc66bde65390894f8729b3ef6dce 100644
--- a/src/gpu/VirtualFluids_GPU/Kernel/Kernels/PorousMediaKernels/Advection/Compressible/CumulantOne/PMCumulantOneCompSP27_Device.cu
+++ b/src/gpu/VirtualFluids_GPU/Kernel/Kernels/PorousMediaKernels/Advection/Compressible/CumulantOne/PMCumulantOneCompSP27_Device.cu
@@ -160,9 +160,9 @@ extern "C" __global__ void LB_Kernel_PM_Cum_One_Comp_SP_27(real omega,
//vvz = (two * vvz) / (two + 134.4 + 0.0068287 * sqrtf(vx2 + vy2 + vz2));
////////////////////////////////////////////////////////////////////////////////////
//the force be with you
- real fx = forces[0] / (pow(c2o1, level)); //zero;//0.0032653/(pow(two,level)); //0.000000005;//(two/1600000.0) / 120.0; //
- real fy = forces[1] / (pow(c2o1, level)); //zero;
- real fz = forces[2] / (pow(c2o1, level)); //zero;
+ real fx = forces[0] / (pow((double)c2o1, (double)level)); //zero;//0.0032653/(pow(two,level)); //0.000000005;//(two/1600000.0) / 120.0; //
+ real fy = forces[1] / (pow((double)c2o1, (double)level)); //zero;
+ real fz = forces[2] / (pow((double)c2o1, (double)level)); //zero;
vvx += fx*c1o2;
vvy += fy*c1o2;
vvz += fz*c1o2;
diff --git a/src/gpu/VirtualFluids_GPU/Kernel/Kernels/WaleKernels/Advection/Compressible/CumulantK15/WaleCumulantK15Comp_Device.cu b/src/gpu/VirtualFluids_GPU/Kernel/Kernels/WaleKernels/Advection/Compressible/CumulantK15/WaleCumulantK15Comp_Device.cu
index 3308121db745108efc31053ebbff0b707edaf8bf..a23d1c480b31befbe1ece6b83b20fa415c94a0e7 100644
--- a/src/gpu/VirtualFluids_GPU/Kernel/Kernels/WaleKernels/Advection/Compressible/CumulantK15/WaleCumulantK15Comp_Device.cu
+++ b/src/gpu/VirtualFluids_GPU/Kernel/Kernels/WaleKernels/Advection/Compressible/CumulantK15/WaleCumulantK15Comp_Device.cu
@@ -434,9 +434,9 @@ extern "C" __global__ void LB_Kernel_WaleCumulantK15Comp(
////////////////////////////////////////////////////////////////////////////////////
//the force be with you
{
- real fx = forces[0] / (pow(c2o1, level)); //zero;//0.0032653/(pow(two,level)); //0.000000005;//(two/1600000.0) / 120.0; //
- real fy = forces[1] / (pow(c2o1, level)); //zero;
- real fz = forces[2] / (pow(c2o1, level)); //zero;
+ real fx = forces[0] / (pow((double)c2o1, (double)level)); //zero;//0.0032653/(pow(two,level)); //0.000000005;//(two/1600000.0) / 120.0; //
+ real fy = forces[1] / (pow((double)c2o1, (double)level)); //zero;
+ real fz = forces[2] / (pow((double)c2o1, (double)level)); //zero;
vvx += fx*c1o2;
vvy += fy*c1o2;
vvz += fz*c1o2;
diff --git a/src/gpu/VirtualFluids_GPU/Kernel/Kernels/WaleKernels/Advection/Compressible/CumulantK15BySoniMalav/WaleBySoniMalavCumulantK15Comp_Device.cu b/src/gpu/VirtualFluids_GPU/Kernel/Kernels/WaleKernels/Advection/Compressible/CumulantK15BySoniMalav/WaleBySoniMalavCumulantK15Comp_Device.cu
index 7e685897a543a83e9fe1e5fdaee9e125140c81df..7f93133cb660c0af5a5d35ab6be37686b14f0e4c 100644
--- a/src/gpu/VirtualFluids_GPU/Kernel/Kernels/WaleKernels/Advection/Compressible/CumulantK15BySoniMalav/WaleBySoniMalavCumulantK15Comp_Device.cu
+++ b/src/gpu/VirtualFluids_GPU/Kernel/Kernels/WaleKernels/Advection/Compressible/CumulantK15BySoniMalav/WaleBySoniMalavCumulantK15Comp_Device.cu
@@ -307,9 +307,9 @@ extern "C" __global__ void LB_Kernel_WaleBySoniMalavCumulantK15Comp(
////////////////////////////////////////////////////////////////////////////////////
//the force be with you
{
- real fx = forces[0] / (pow(c2o1, level)); //zero;//0.0032653/(pow(two,level)); //0.000000005;//(two/1600000.0) / 120.0; //
- real fy = forces[1] / (pow(c2o1, level)); //zero;
- real fz = forces[2] / (pow(c2o1, level)); //zero;
+ real fx = forces[0] / (pow((double)c2o1, (double)level)); //zero;//0.0032653/(pow(two,level)); //0.000000005;//(two/1600000.0) / 120.0; //
+ real fy = forces[1] / (pow((double)c2o1, (double)level)); //zero;
+ real fz = forces[2] / (pow((double)c2o1, (double)level)); //zero;
vvx += fx*c1o2;
vvy += fy*c1o2;
vvz += fz*c1o2;
diff --git a/src/gpu/VirtualFluids_GPU/Kernel/Kernels/WaleKernels/Advection/Compressible/CumulantK17/WaleCumulantK17Comp_Device.cu b/src/gpu/VirtualFluids_GPU/Kernel/Kernels/WaleKernels/Advection/Compressible/CumulantK17/WaleCumulantK17Comp_Device.cu
index 6e9939b792e0644cb20feb207a7df49fc2c8d9d5..91e91bc0ff5e1582eab279c77e9fa0b419ce06ea 100644
--- a/src/gpu/VirtualFluids_GPU/Kernel/Kernels/WaleKernels/Advection/Compressible/CumulantK17/WaleCumulantK17Comp_Device.cu
+++ b/src/gpu/VirtualFluids_GPU/Kernel/Kernels/WaleKernels/Advection/Compressible/CumulantK17/WaleCumulantK17Comp_Device.cu
@@ -435,9 +435,9 @@ extern "C" __global__ void LB_Kernel_WaleCumulantK17Comp(
////////////////////////////////////////////////////////////////////////////////////
//the force be with you
{
- real fx = forces[0] / (pow(c2o1, level)); //zero;//0.0032653/(pow(two,level)); //0.000000005;//(two/1600000.0) / 120.0; //
- real fy = forces[1] / (pow(c2o1, level)); //zero;
- real fz = forces[2] / (pow(c2o1, level)); //zero;
+ real fx = forces[0] / (pow((double)c2o1, (double)level)); //zero;//0.0032653/(pow(two,level)); //0.000000005;//(two/1600000.0) / 120.0; //
+ real fy = forces[1] / (pow((double)c2o1, (double)level)); //zero;
+ real fz = forces[2] / (pow((double)c2o1, (double)level)); //zero;
vvx += fx*c1o2;
vvy += fy*c1o2;
vvz += fz*c1o2;
diff --git a/src/gpu/VirtualFluids_GPU/Kernel/Kernels/WaleKernels/Advection/Compressible/CumulantK17Debug/WaleCumulantK17DebugComp_Device.cu b/src/gpu/VirtualFluids_GPU/Kernel/Kernels/WaleKernels/Advection/Compressible/CumulantK17Debug/WaleCumulantK17DebugComp_Device.cu
index 1f371fea7470f156a8685972775b4e0c8ca4525f..b0e0b65064e2395f5049eeb2df3b8260125fbf4e 100644
--- a/src/gpu/VirtualFluids_GPU/Kernel/Kernels/WaleKernels/Advection/Compressible/CumulantK17Debug/WaleCumulantK17DebugComp_Device.cu
+++ b/src/gpu/VirtualFluids_GPU/Kernel/Kernels/WaleKernels/Advection/Compressible/CumulantK17Debug/WaleCumulantK17DebugComp_Device.cu
@@ -454,9 +454,9 @@ extern "C" __global__ void LB_Kernel_WaleCumulantK17DebugComp(
////////////////////////////////////////////////////////////////////////////////////
//the force be with you
{
- real fx = forces[0] / (pow(c2o1, level)); //zero;//0.0032653/(pow(two,level)); //0.000000005;//(two/1600000.0) / 120.0; //
- real fy = forces[1] / (pow(c2o1, level)); //zero;
- real fz = forces[2] / (pow(c2o1, level)); //zero;
+ real fx = forces[0] / (pow((double)c2o1, (double)level)); //zero;//0.0032653/(pow(two,level)); //0.000000005;//(two/1600000.0) / 120.0; //
+ real fy = forces[1] / (pow((double)c2o1, (double)level)); //zero;
+ real fz = forces[2] / (pow((double)c2o1, (double)level)); //zero;
vvx += fx*c1o2;
vvy += fy*c1o2;
vvz += fz*c1o2;
diff --git a/src/gpu/VirtualFluids_GPU/LBM/Simulation.cpp b/src/gpu/VirtualFluids_GPU/LBM/Simulation.cpp
index 1de408a5d5991fd8ceee59e863d5e32e6384fe67..a4dbc09a99fb4e5ffe434b6d0bfeb27a73da1052 100644
--- a/src/gpu/VirtualFluids_GPU/LBM/Simulation.cpp
+++ b/src/gpu/VirtualFluids_GPU/LBM/Simulation.cpp
@@ -19,7 +19,7 @@
#include "Output/VeloASCIIWriter.hpp"
//////////////////////////////////////////////////////////////////////////
#include "Utilities/Buffer2D.hpp"
-#include "Utilities/StringUtil.hpp"
+#include "Core/StringUtilities/StringUtil.h"
//////////////////////////////////////////////////////////////////////////
#include "Init/InitLattice.h"
#include "Init/DefineGrid.h"
diff --git a/src/gpu/VirtualFluids_GPU/Output/AnalysisData.hpp b/src/gpu/VirtualFluids_GPU/Output/AnalysisData.hpp
index 8163a612f7bcf77884d6f8c78c8e0c09a470daad..9263f5009f6bda53966197bd7076dab23093c8d9 100644
--- a/src/gpu/VirtualFluids_GPU/Output/AnalysisData.hpp
+++ b/src/gpu/VirtualFluids_GPU/Output/AnalysisData.hpp
@@ -3,7 +3,7 @@
#include "basics/utilities/UbFileOutputASCII.h"
#include "Parameter/Parameter.h"
-#include "Utilities/StringUtil.hpp"
+#include "Core/StringUtilities/StringUtil.h"
class AnalysisData
{
diff --git a/src/gpu/VirtualFluids_GPU/Output/InterfaceDebugWriter.hpp b/src/gpu/VirtualFluids_GPU/Output/InterfaceDebugWriter.hpp
index c278f5349cb4bb2d0ce7f758733e143adbaea37a..86f380705a3b89105b084fe40a7ffbae65cb10e6 100644
--- a/src/gpu/VirtualFluids_GPU/Output/InterfaceDebugWriter.hpp
+++ b/src/gpu/VirtualFluids_GPU/Output/InterfaceDebugWriter.hpp
@@ -10,7 +10,7 @@
#include "LBM/D3Q27.h"
#include "Parameter/Parameter.h"
#include "basics/utilities/UbSystem.h"
-#include "Utilities/StringUtil.hpp"
+#include "Core/StringUtilities/StringUtil.h"
#include <basics/writer/WbWriterVtkXmlBinary.h>
diff --git a/src/gpu/VirtualFluids_GPU/Output/QDebugWriter.hpp b/src/gpu/VirtualFluids_GPU/Output/QDebugWriter.hpp
index 583f36730320363dddc52bdce435cbb1f9095bf8..d661aef500b745b49d28892a01fae96e5a6d5dd3 100644
--- a/src/gpu/VirtualFluids_GPU/Output/QDebugWriter.hpp
+++ b/src/gpu/VirtualFluids_GPU/Output/QDebugWriter.hpp
@@ -10,9 +10,8 @@
#include "LBM/D3Q27.h"
#include "Parameter/Parameter.h"
#include "basics/utilities/UbSystem.h"
-#include "Utilities/StringUtil.hpp"
#include <basics/writer/WbWriterVtkXmlBinary.h>
-
+#include "Core/StringUtilities/StringUtil.h"
//using namespace std;
diff --git a/src/gpu/VirtualFluids_GPU/Particles/Particles.h b/src/gpu/VirtualFluids_GPU/Particles/Particles.h
index 4dcea815aa1ebb200e2174f32d5c85dcede53d9b..ccffb11e7afe2a5dce899644d2c65a55feed529f 100644
--- a/src/gpu/VirtualFluids_GPU/Particles/Particles.h
+++ b/src/gpu/VirtualFluids_GPU/Particles/Particles.h
@@ -4,7 +4,7 @@
#include "LBM/LB.h"
#include "GPU/GPU_Interface.h"
#include "GPU/CudaMemoryManager.h"
-#include "Utilities/StringUtil.hpp"
+#include "Core/StringUtilities/StringUtil.h"
#include "Parameter/Parameter.h"
//extern "C" void calcDragLift(Parameter* para, int lev);
diff --git a/src/gpu/VirtualFluids_GPU/Utilities/StringUtil.hpp b/src/gpu/VirtualFluids_GPU/Utilities/StringUtil.hpp
deleted file mode 100644
index 46b2fccd74f3bc0a1ed207d2d0a3d4b8d74ef834..0000000000000000000000000000000000000000
--- a/src/gpu/VirtualFluids_GPU/Utilities/StringUtil.hpp
+++ /dev/null
@@ -1,299 +0,0 @@
-#ifndef STRINGUTIL_H
-#define STRINGUTIL_H
-#include <errno.h>
-#include <algorithm>
-#include <sstream>
-#include <iostream>
-#include <string>
-#include <boost/foreach.hpp>
-#include <boost/algorithm/string.hpp>
-
-template <class T>
-bool from_string(T &t, const std::string &s, std::ios_base &(*f)(std::ios_base&))
-{
- std::istringstream iss(s);
- return !(iss>>f>>t).fail();
-}
-class StringUtil
-{
-public:
- StringUtil() {}
- ~StringUtil() {}
- // Find the given string in the source string and replace it with the
- // "replace" string, everywhere instances of that string exist.
- static void findandreplace( std::string& source, const std::string& find, const std::string& replace )
- {
- size_t j;
- for (;(j = source.find( find ))
- != std::string::npos;)
- {
- source.replace( j,
- find.length(), replace );
- }
- }
- // The following function returns a string with all-uppercase characters.
- static std::string makeUpper( const std::string& instring)
- {
- std::string temp=instring;
- transform( temp.begin(), temp.end(), temp.begin(), ::toupper );
- return temp;
- }
- // The following function returns a string with all-lowercase characters.
- static std::string makeLower( const std::string& instring)
- {
- std::string temp;
- transform( temp.begin(), temp.end(), temp.begin(), ::tolower );
- return temp;
- }
- static bool contains( const std::string& source, const char *find )
- {
- return ( 0!=strstr(source.c_str(),find) );
- }
- static std::string pad( const std::string& instring, char padchar, int length )
- {
- std::string outstring = instring;
- for ( int i=(int)outstring.length(); i<length; ++i )
- outstring += padchar;
- return outstring;
- }
- // Trim the given characters from the beginning and end of a string.
- // the default is to trim whitespace. If the string is empty or contains
- // only the trim characters, an empty string is returned.
- static std::string trim( const std::string &instring, const std::string &trimstring=std::string(" \t\n"))
- {
- if (trimstring.size()==0)
- return instring;
- std::string temp="";
- std::string::size_type begpos=instring.find_first_not_of (trimstring);
- if (begpos==std::string::npos)
- {
- return temp;
- }
- else
- {
- std::string::size_type endpos=instring.find_last_not_of (trimstring);
- temp=instring.substr(begpos, endpos-begpos+1);
- }
- return temp;
- }
- // Convert the string to an int. Note that a string exception is thrown if
- // it is invalid.
- static int toInt(const std::string & myInString)
- {
- int i=0;
- std::string inString = trim(myInString);
- if( !from_string<int>(i, inString, std::dec) )
- {
- std::string exceptionText = "StringUtils::toInt() - Not an integer: " + inString;
- throw exceptionText;
- }
- // Time to run some more checks.
- for (unsigned int j=0; j < inString.length(); j++)
- {
- if ( !isNumeric(inString[j]) )
- {
- if (j==0 && inString[j] =='-')
- {
- continue;
- }
- else
- {
- std::string exceptionText = "StringUtils::toInt() - Not an integer: " + inString;
- throw exceptionText;
- }
- }
- }
- return (i);
- }
- // Convert the string to a float. Note: A string exception is thrown if
- // it is invalid.
- static float toFloat(const std::string & myInString)
- {
- float f=0;
- std::string inString = trim(myInString);
- if( !from_string<float>(f, inString, std::dec) )
- {
- std::string exceptionText = "StringUtils::toFloat() - Not a float: " + inString;
- throw exceptionText;
- }
- // Now it runs some more checks.
- int dec_count=0;
- int e_count=0;
- for (unsigned int j=0; j < inString.length(); j++)
- {
- if ( !isNumeric(inString[j]) )
- {
- if ((j==0 || inString[j-1] == 'e' || inString[j-1] == 'E') && inString[j] =='-')
- {
- continue;
- }
- else if (inString[j]=='.')
- {
- dec_count++;
- if (dec_count > 1)
- {
- std::string exceptionText = "StringUtils::toFloat() - Not a float: " + inString;
- throw exceptionText;
- }
- continue;
- }
- else if (inString[j] == 'e' || inString[j] == 'E')
- {
- e_count++;
- if (e_count > 1)
- {
- std::string exceptionText = "StringUtils::toFloat() - Not a float: " + inString;
- throw exceptionText;
- }
- continue;
- }
- else
- {
- std::string exceptionText = "StringUtils::toFloat() - Not a float: " + inString;
- throw exceptionText;
- }
- }
- }
- return (f);
- }
- // Convert the string to a double. Note: A string exception is thrown if
- // it is invalid.
- static double toDouble(const std::string & myInString)
- {
- double d=0;
- std::string inString = trim(myInString);
- if( !from_string<double>(d, inString, std::dec) )
- {
- std::string exceptionText = "StringUtils::toDouble() - Not a double: " + inString;
- throw exceptionText;
- }
- // Now it runs some more checks.
- int dec_count=0;
- int e_count=0;
- for (unsigned int j=0; j < inString.length(); j++)
- {
- if ( !isNumeric(inString[j]) )
- {
- if ((j==0 || inString[j-1] == 'e' || inString[j-1] == 'E') && inString[j] =='-')
- {
- continue;
- }
- else if (inString[j]=='.')
- {
- dec_count++;
- if (dec_count > 1)
- {
- std::string exceptionText = "StringUtils::toDouble() - Not a double: " + inString;
- throw exceptionText;
- }
- continue;
- }
- else if (inString[j] == 'e' || inString[j] == 'E')
- {
- e_count++;
- if (e_count > 1)
- {
- std::string exceptionText = "StringUtils::toDouble() - Not a double: " + inString;
- throw exceptionText;
- }
- continue;
- }
- else
- {
- std::string exceptionText = "StringUtils::toDouble() - Not a double: " + inString;
- throw exceptionText;
- }
- }
- }
- return (d);
- }
- // Convert the string to a boolean. Note: A string exception is thrown if
- // it is invalid.
- static bool toBool(const std::string & myInString)
- {
- bool b=0;
- std::string inString = trim(myInString);
- if( !from_string<bool>(b, inString, std::boolalpha) )
- {
- std::string exceptionText = "StringUtils::toBool() - Not a bool: " + inString;
- throw exceptionText;
- }
- return (b);
- }
- // Returns true if the character is numeric.
- static bool isNumeric(char c)
- {
- return ('0' <= c && c <= '9');
- }
- // Replace environment variables in the string with their values.
- // Note: environment variables must be of the form ${ENVVAR}.
- //static std::string substituteEnvVar( const std::string &myInString )
- //{
- // std::string outString="";
- // char variable[512];
- // const char *s = myInString.c_str();
- // while(*s!=0)
- // {
- // if (*s=='$' && *(s+1)=='{')
- // {
- // // When you�ve found beginning of variable, find the end.
- // //strcpy(variable,s+2);
- // strcpy_s(variable, s+2);
- // char *end = strchr (variable,'}');
- // if (end)
- // {
- // *end='\0';
- // //char *cp = (char *)getenv(variable);
- // char *cp;
- // size_t len;
- // _dupenv_s(&cp, &len, variable );
- // if (cp)
- // //outString += (char *) getenv(variable);
- // outString += (char *) cp;
- // s = strchr(s,'}');
- // }
- // else
- // {
- // outString += *s;
- // }
- // }
- // else
- // {
- // outString += *s;
- // }
- // s++;
- // }
- // return outString;
- //}
- template<class T>
- static std::vector<T> toVector(const std::string& s)
- {
- std::vector<T> v;
- std::vector<std::string> strings;
- boost::algorithm::split(strings, s, boost::is_any_of("\t\n "));
- BOOST_FOREACH(std::string s, strings){
- if (s != "")
- {
- v.push_back(fromString<T>(s));
- }
- }
- return v;
- }
- template<class T>
- static std::string toString(const T& t)
- {
- std::ostringstream stream;
- stream << t;
- return stream.str();
- }
- template<class T>
- static T fromString(const std::string& s)
- {
- std::istringstream stream (s);
- T t;
- stream >> t;
- return t;
- }
-};
-#endif //STRINGUTIL_H
-