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

made changed to accretion_profile function, still not fully finished but close.

Browse Source
This commit is contained in:
JackHsu
2019-08-07 13:59:12 +08:00
parent 2a26b58c38
commit 3adec853df
1 changed files with 9 additions and 10 deletions
Download Patch File Download Diff File Expand all files Collapse all files

19
acc/mhd_solver/stencil_process.sps
View File

@@ -79,7 +79,6 @@ sink_gravity(int3 globalVertexIdx){
// Now this is more of a suedo-code, just write out some intuitions.
#if LSINK #if LSINK
Scalar Scalar
truelove_density(in Scalar lnrho){ truelove_density(in Scalar lnrho){
@@ -90,7 +89,7 @@ truelove_density(in Scalar lnrho){
Scalar accrection_rho = rho - TJ_rho; Scalar accrection_rho = rho - TJ_rho;
if (accrection_rho < 0){ if (accrection_rho < 0){
accrection_rho = Scalar(0); accrection_rho = Scalar(0);
} }
return accretion_rho; return accretion_rho;
@@ -98,21 +97,21 @@ truelove_density(in Scalar lnrho){
accretion_profile(int3 globalVertexIdx){ accretion_profile(int3 globalVertexIdx){
// ONE QUESTION: do I need to define grid_pos, sink_pos and distance again
// if the sink_gravity kernel will also be called once LSINK swtich is on? Seems redundant.
const Vector grid_pos = (Vector){(globalVertexIdx.x - nx_min) * dsx, const Vector grid_pos = (Vector){(globalVertexIdx.x - nx_min) * dsx,
(globalVertexIdx.y - ny_min) * dsy, (globalVertexIdx.y - ny_min) * dsy,
(globalVertexIdx.z - nz_min) * dsz}; (globalVertexIdx.z - nz_min) * dsz};
const Vector sink_pos = (Vector){DCONST_REAL(AC_sink_pos_x), const Vector sink_pos = (Vector){DCONST_REAL(AC_sink_pos_x),
DCONST_REAL(AC_sink_pos_y), DCONST_REAL(AC_sink_pos_y),
DCONST_REAL(AC_sink_pos_z)}; DCONST_REAL(AC_sink_pos_z)};
// Here I'm trying to make the accretion profile adjustable within astaroth.conf, i.e., making it possible to const Scalar profile_range = DCONST_REAL(AC_accretion_range) * dsx;
// choose how many "cells" around the sink particle we want to accrete mass from. const Scalar accretion_distance = length(grid_pos - sink_pos);
const Scalar profile = DCONST_REAL(AC_accretion_range) * dsx; if ((accretion_distance) <= profile_range){
if ((grid_pos - sink_pos) <= profile){ // calculate accretion according to chosen criterion for the grid cell.
accretion = accretion_rho * dsx * dsy * dsz;
} }
// I'm still trying to figure out what to put in if statement when the cell is within the range we set return accretion;// the returned value is effectively mass, doesn't seem right since MIIKKA said mass will be summed up in vertex buffer. But I'll keep it as it untill we discuss this.
// and what should this kernel return.
return
} }