This commit is contained in:
jpekkila
2019-06-19 14:29:59 +03:00
5 changed files with 39 additions and 8 deletions

View File

@@ -545,10 +545,12 @@ normalized(const AcReal3& vec)
// Sinusoidal forcing
// https://arxiv.org/pdf/1704.04676.pdf
// NOTE: This method of forcing is depracated. However, it will remain in here
// until a corresponding scheme exists in the new code.
__constant__ AcReal3 forcing_vec;
__constant__ AcReal forcing_phi;
static __device__ __forceinline__ AcReal3
forcing(const int i, const int j, const int k)
DEPRECATED_forcing(const int i, const int j, const int k)
{
#define DOMAIN_SIZE_X (DCONST_INT(AC_nx) * DCONST_REAL(AC_dsx))
#define DOMAIN_SIZE_Y (DCONST_INT(AC_ny) * DCONST_REAL(AC_dsy))

View File

@@ -202,6 +202,7 @@ run_simulation(void)
load_config(&mesh_info);
AcMesh* mesh = acmesh_create(mesh_info);
//TODO: This need to be possible to define in astaroth.conf
acmesh_init_to(INIT_TYPE_GAUSSIAN_RADIAL_EXPL, mesh);
acInit(mesh_info);