/******************************************************************************* * This file is part of SWIFT. * Copyright (c) 2016 Matthieu Schaller (schaller@strw.leidenuniv.nl) * * This program is free software: you can redistribute it and/or modify * it under the terms of the GNU Lesser General Public License as published * by the Free Software Foundation, either version 3 of the License, or * (at your option) any later version. * * This program is distributed in the hope that it will be useful, * but WITHOUT ANY WARRANTY; without even the implied warranty of * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the * GNU General Public License for more details. * * You should have received a copy of the GNU Lesser General Public License * along with this program. If not, see . * ******************************************************************************/ #ifndef SWIFT_GADGET2_HYDRO_DEBUG_H #define SWIFT_GADGET2_HYDRO_DEBUG_H __attribute__((always_inline)) INLINE static void hydro_debug_particle( const struct part* p, const struct xpart* xp) { warning("[PID%lld] part:", p->id); warning( "[PID%lld] x=[%.3e,%.3e,%.3e], " "v=[%.3e,%.3e,%.3e]", p->id, p->x[0], p->x[1], p->x[2], p->v[0], p->v[1], p->v[2]); warning( "[PID%lld] a=[%.3e,%.3e,%.3e] 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", p->id, p->a_hydro[0], p->a_hydro[1], p->a_hydro[2], p->h, p->density.wcount, p->density.wcount_dh, p->mass, p->density.rho_dh, p->rho, hydro_get_comoving_pressure(p), p->force.P_over_rho2, p->entropy, p->entropy_dt, p->force.soundspeed); warning("[PID%lld] divV=%.3e, rotV=[%.3e,%.3e,%.3e], balsara=%.3e", p->id, p->density.div_v, p->density.rot_v[0], p->density.rot_v[1], p->density.rot_v[2], p->force.balsara); warning("[PID%lld] v_sig=%e dh/dt=%.3e time_bin=%d wakeup=%d", p->id, p->force.v_sig, p->force.h_dt, p->time_bin, p->limiter_data.wakeup); if (xp != NULL) { warning("[PID%lld] xpart:", p->id); warning("[PID%lld] v_full=[%.3e,%.3e,%.3e]", p->id, xp->v_full[0], xp->v_full[1], xp->v_full[2]); } } #endif /* SWIFT_GADGET2_HYDRO_DEBUG_H */