Commit 8e1ea32d authored by Matthieu Schaller's avatar Matthieu Schaller
Browse files

Added pressure to the debug information

parent 9bb68e1a
......@@ -22,18 +22,15 @@ __attribute__((always_inline)) INLINE static void hydro_debug_particle(
printf(
"x=[%.3e,%.3e,%.3e], "
"v=[%.3e,%.3e,%.3e],v_full=[%.3e,%.3e,%.3e] \n a=[%.3e,%.3e,%.3e],\n "
"h=%.3e, "
"wcount=%d, wcount_dh=%.3e, m=%.3e, dh_drho=%.3e, rho=%.3e, "
"P_over_rho2=%.3e, "
"S=%.3e, "
"dS/dt=%.3e, c=%.3e\n"
"h=%.3e, wcount=%.3f, wcount_dh=%.3e, m=%.3e, dh_drho=%.3e, rho=%.3e, "
"P=%.3e, P_over_rho2=%.3e, S=%.3e, dS/dt=%.3e, c=%.3e\n"
"divV=%.3e, rotV=[%.3e,%.3e,%.3e], balsara=%.3e \n "
"v_sig=%e dh/dt=%.3e t_begin=%d, t_end=%d\n",
p->x[0], p->x[1], p->x[2], p->v[0], p->v[1], p->v[2], xp->v_full[0],
xp->v_full[1], xp->v_full[2], p->a_hydro[0], p->a_hydro[1], p->a_hydro[2],
p->h, (int)p->density.wcount, p->density.wcount_dh, p->mass, p->rho_dh,
p->rho, p->force.P_over_rho2, p->entropy, p->force.entropy_dt,
p->force.soundspeed, p->density.div_v, p->density.rot_v[0],
p->density.rot_v[1], p->density.rot_v[2], p->force.balsara,
p->force.v_sig, p->force.h_dt, p->ti_begin, p->ti_end);
p->h, p->density.wcount, p->density.wcount_dh, p->mass, p->rho_dh, p->rho,
p->force.P_over_rho2 * p->rho * p->rho / p->rho_dh, p->force.P_over_rho2,
p->entropy, p->force.entropy_dt, p->force.soundspeed, p->density.div_v,
p->density.rot_v[0], p->density.rot_v[1], p->density.rot_v[2],
p->force.balsara, p->force.v_sig, p->force.h_dt, p->ti_begin, p->ti_end);
}
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment