Now initilize VTXBUFF_ACCRETION to 0.0, and accretion now works as expected.

This commit is contained in:
JackHsu
2019-08-12 12:34:38 +08:00
parent 89128af44b
commit b782ead4f3
3 changed files with 5 additions and 2 deletions

View File

@@ -71,7 +71,7 @@ acmesh_create(const AcMeshInfo& mesh_info)
return mesh; return mesh;
} }
static void void
vertex_buffer_set(const VertexBufferHandle& key, const AcReal& val, AcMesh* mesh) vertex_buffer_set(const VertexBufferHandle& key, const AcReal& val, AcMesh* mesh)
{ {
const int n = acVertexBufferSize(mesh->info); const int n = acVertexBufferSize(mesh->info);

View File

@@ -48,6 +48,8 @@ extern const char* init_type_names[]; // Defined in host_memory.cc
AcMesh* acmesh_create(const AcMeshInfo& mesh_info); AcMesh* acmesh_create(const AcMeshInfo& mesh_info);
void vertex_buffer_set(const VertexBufferHandle& key, const AcReal& val, AcMesh* mesh);
void acmesh_clear(AcMesh* mesh); void acmesh_clear(AcMesh* mesh);
void acmesh_init_to(const InitType& type, AcMesh* mesh); void acmesh_init_to(const InitType& type, AcMesh* mesh);

View File

@@ -96,7 +96,6 @@ write_mesh_info(const AcMeshInfo* config)
fprintf(infotxt, "real AC_sink_pos_z %e \n", (double)config->real_params[AC_sink_pos_z]); fprintf(infotxt, "real AC_sink_pos_z %e \n", (double)config->real_params[AC_sink_pos_z]);
fprintf(infotxt, "real AC_M_sink %e \n", (double)config->real_params[AC_M_sink]); fprintf(infotxt, "real AC_M_sink %e \n", (double)config->real_params[AC_M_sink]);
fprintf(infotxt, "real AC_soft %e \n", (double)config->real_params[AC_soft]); fprintf(infotxt, "real AC_soft %e \n", (double)config->real_params[AC_soft]);
fprintf(infotxt, "real AC_G_const %e \n", (double)config->real_params[AC_G_const]); fprintf(infotxt, "real AC_G_const %e \n", (double)config->real_params[AC_G_const]);
fclose(infotxt); fclose(infotxt);
@@ -196,6 +195,8 @@ run_simulation(void)
// TODO: This need to be possible to define in astaroth.conf // TODO: This need to be possible to define in astaroth.conf
acmesh_init_to(INIT_TYPE_GAUSSIAN_RADIAL_EXPL, mesh); acmesh_init_to(INIT_TYPE_GAUSSIAN_RADIAL_EXPL, mesh);
vertex_buffer_set(VTXBUF_ACCRETION, 0.0, mesh);
acInit(mesh_info); acInit(mesh_info);
acLoad(*mesh); acLoad(*mesh);