cwpearson/astaroth
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Completed function to sum VTXBUF_ACCRETION and updates mass of sink particle. Needs further work with crashing problem.

Browse Source
This commit is contained in:
JackHsu
2019-08-12 19:02:22 +08:00
parent e7ff447510
commit 8c7176d702
3 changed files with 10 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

1
acc/mhd_solver/stencil_defines.h
View File

@@ -78,6 +78,7 @@
FUNC(AC_sink_pos_y),\ FUNC(AC_sink_pos_y),\
FUNC(AC_sink_pos_z),\ FUNC(AC_sink_pos_z),\
FUNC(AC_M_sink),\ FUNC(AC_M_sink),\
FUNC(AC_M_sink_init),\
FUNC(AC_M_sink_Msun),\ FUNC(AC_M_sink_Msun),\
FUNC(AC_soft),\ FUNC(AC_soft),\
FUNC(AC_accretion_range),\ FUNC(AC_accretion_range),\

2
src/standalone/config_loader.cc
View File

@@ -150,6 +150,8 @@ update_config(AcMeshInfo* config)
config->real_params[AC_M_sink] = config->real_params[AC_M_sink_Msun] * M_sun / config->real_params[AC_M_sink] = config->real_params[AC_M_sink_Msun] * M_sun /
config->real_params[AC_unit_mass]; config->real_params[AC_unit_mass];
config->real_params[AC_M_sink_init] = config->real_params[AC_M_sink_Msun] * M_sun /
config->real_params[AC_unit_mass];
config->real_params[AC_G_const] = G_CONST_CGS / ((config->real_params[AC_unit_velocity] * config->real_params[AC_G_const] = G_CONST_CGS / ((config->real_params[AC_unit_velocity] *
config->real_params[AC_unit_velocity]) / config->real_params[AC_unit_velocity]) /

7
src/standalone/simulation.cc
View File

@@ -240,9 +240,16 @@ run_simulation(void)
// acUpdate_sink_particle() will do the similar trick to the device. // acUpdate_sink_particle() will do the similar trick to the device.
/* Step the simulation */ /* Step the simulation */
AcReal accreted_mass = 0.0;
for (int i = 1; i < max_steps; ++i) { for (int i = 1; i < max_steps; ++i) {
const AcReal umax = acReduceVec(RTYPE_MAX, VTXBUF_UUX, VTXBUF_UUY, VTXBUF_UUZ); const AcReal umax = acReduceVec(RTYPE_MAX, VTXBUF_UUX, VTXBUF_UUY, VTXBUF_UUZ);
const AcReal dt = host_timestep(umax, mesh_info); const AcReal dt = host_timestep(umax, mesh_info);
const AcReal sum_mass = acReduceScal(RTYPE_MAX, VTXBUF_ACCRETION);
accreted_mass = accreted_mass + sum_mass;
AcReal sink_mass = AC_M_sink_init + accreted_mass;
printf("sink mass is: %e\n", sink_mass);
printf("accreted mass is: %e\n", accreted_mass);
acLoadDeviceConstant(AC_M_sink, sink_mass);
#if LFORCING #if LFORCING
const ForcingParams forcing_params = generateForcingParams(mesh_info); const ForcingParams forcing_params = generateForcingParams(mesh_info);