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

Fixed some syntax errors and changed LGRAVITY to LSINK. Will test if this works later.

Browse Source
This commit is contained in:
JackHsu
2019-07-30 11:47:18 +08:00
parent e7a2734ab6
commit 0acf52d1b4
1 changed files with 10 additions and 12 deletions
Download Patch File Download Diff File Expand all files Collapse all files

22
acc/mhd_solver/stencil_process.sps
View File

@@ -65,14 +65,14 @@ continuity(in Vector uu, in Scalar lnrho) {
Vector Vector
gravity() { gravity() {
vector force_gravity; vector force_gravity;
const int3 grid_pos = (int3){threadIdx.x + blockIdx.x * blockDim.x + start.x, const Vector grid_pos = (Vector){(globalVertexIdx.x - nx_min) * dsx,
threadIdx.y + blockIdx.y * blockDim.y + start.y, (globalVertexIdx.y - ny_min) * dsy,
threadIdx.z + blockIdx.z * blockDim.z + start.z}; (globalVertexIdx.z - nz_min) * dsz};
const sink_mass = DCONST(AC_M_sink); const Scalar sink_mass = DCONST(AC_M_sink);
const 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)};
const distance = dlength_vec(grid_pos - sink_pos); const Scalar distance = length_vec(grid_pos - sink_pos);
force_gravity = (6.67e-8 * sink_mass) / (distance * distance); force_gravity = (6.67e-8 * sink_mass) / (distance * distance);
return force_gravity; return force_gravity;
}//first attempt to do a self-containing LSINK module }//first attempt to do a self-containing LSINK module
@@ -115,8 +115,8 @@ momentum(in Vector uu, in Scalar lnrho, in Scalar tt) {
(laplace_vec(uu) + Scalar(1. / 3.) * gradient_of_divergence(uu) + (laplace_vec(uu) + Scalar(1. / 3.) * gradient_of_divergence(uu) +
Scalar(2.) * mul(S, gradient(lnrho))) + zeta * gradient_of_divergence(uu); Scalar(2.) * mul(S, gradient(lnrho))) + zeta * gradient_of_divergence(uu);
#if LGRAVITY #if LSINK
mom = mom - (Vector){0, 0, -10.0}; mom = mom - (Vector){force_gravity};
#endif #endif
return mom; return mom;
@@ -136,8 +136,8 @@ momentum(in Vector uu, in Scalar lnrho) {
(laplace_vec(uu) + Scalar(1. / 3.) * gradient_of_divergence(uu) + (laplace_vec(uu) + Scalar(1. / 3.) * gradient_of_divergence(uu) +
Scalar(2.) * mul(S, gradient(lnrho))) + zeta * gradient_of_divergence(uu); Scalar(2.) * mul(S, gradient(lnrho))) + zeta * gradient_of_divergence(uu);
#if LGRAVITY #if LSINK
mom = mom - (Vector){0, 0, -10.0}; mom = mom - (Vector){force_gravity};
#endif #endif
return mom; return mom;
@@ -181,8 +181,6 @@ heat_conduction( in Scalar ss, in Scalar lnrho) {
(gamma - AcReal(1.)) * laplace(lnrho); (gamma - AcReal(1.)) * laplace(lnrho);
const Vector second_term = gamma * inv_cp_sound * gradient(ss) + const Vector second_term = gamma * inv_cp_sound * gradient(ss) +
(gamma - AcReal(1.)) * gradient(lnrho); (gamma - AcReal(1.)) * gradient(lnrho);
const Vector third_term = gamma * (inv_cp_sound * gradient(ss) +
gradient(lnrho)) + grad_ln_chi;
const Scalar chi = AC_THERMAL_CONDUCTIVITY / (exp(value(lnrho)) * cp_sound); const Scalar chi = AC_THERMAL_CONDUCTIVITY / (exp(value(lnrho)) * cp_sound);
return cp_sound * chi * (first_term + dot(second_term, third_term)); return cp_sound * chi * (first_term + dot(second_term, third_term));