/******************************************************************************* * This file is part of SWIFT. * Copyright (c) 2022 Bert Vandenbroucke (bert.vandenbroucke@gmail.com) * * 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_CHEMISTRY_COLIBRE_DEBUG_H #define SWIFT_CHEMISTRY_COLIBRE_DEBUG_H __attribute__((always_inline)) INLINE static void chemistry_debug_particle( const struct part* p, const struct xpart* xp) { warning("[PID%lld] chemistry_part_data:", p->id); for (int i = 0; i < chemistry_element_count; i++) { warning("[PID%lld] metal_mass_fraction[%i]=%.3e", p->id, i, p->chemistry_data.metal_mass_fraction[i]); } warning( "[PID%lld] " "metal_mass_fraction_total=%.3e, " "mass_fraction_from_SNIa=%.3e, " "metal_mass_fraction_from_SNIa=%.3e, " "mass_fraction_from_AGB=%.3e, " "metal_mass_fraction_from_AGB=%.3e, " "mass_fraction_from_SNII=%.3e, " "mass_fraction_from_NSM=%.3e, " "mass_fraction_from_CEJSN=%.3e, " "mass_fraction_from_collapsar=%.3e, " "dmass_fraction_from_NSM=%.3e, " "dmass_fraction_from_CEJSN=%.3e, " "dmass_fraction_from_collapsar=%.3e, " "metal_mass_fraction_from_SNII=%.3e, " "iron_mass_fraction_from_SNIa=%.3e, " "metal_weighted_redshift=%.3e, " "metal_diffused_redshift=%.3e, " "iron_weighted_redshift=%.3e, " "iron_diffused_redshift=%.3e, " "track_of_metal_mass_total=%.3e, " "track_of_iron_mass=%.3e", p->id, p->chemistry_data.metal_mass_fraction_total, p->chemistry_data.mass_fraction_from_SNIa, p->chemistry_data.metal_mass_fraction_from_SNIa, p->chemistry_data.mass_fraction_from_AGB, p->chemistry_data.metal_mass_fraction_from_AGB, p->chemistry_data.mass_fraction_from_SNII, p->chemistry_data.mass_fraction_from_NSM, p->chemistry_data.mass_fraction_from_CEJSN, p->chemistry_data.mass_fraction_from_collapsar, p->chemistry_data.dmass_fraction_from_NSM, p->chemistry_data.dmass_fraction_from_CEJSN, p->chemistry_data.dmass_fraction_from_collapsar, p->chemistry_data.metal_mass_fraction_from_SNII, p->chemistry_data.iron_mass_fraction_from_SNIa, p->chemistry_data.metal_weighted_redshift, p->chemistry_data.metal_diffused_redshift, p->chemistry_data.iron_weighted_redshift, p->chemistry_data.iron_diffused_redshift, p->chemistry_data.track_of_metal_mass_total, p->chemistry_data.track_of_iron_mass); for (int i = 0; i < chemistry_element_count; i++) { warning("[PID%lld] dmetal_mass_fraction[%i]=%.3e", p->id, i, p->chemistry_data.dmetal_mass_fraction[i]); } warning( "[PID%lld] " "dmetal_mass_fraction_total=%.3e, " "dmetal_mass_fraction_from_SNIa=%.3e, " "dmetal_mass_fraction_from_AGB=%.3e, " "dmetal_mass_fraction_from_SNII=%.3e, " "diron_mass_fraction_from_SNIa=%.3e, " "shear_tensor=[[%.3e,%.3e,%.3e],[%.3e,%.3e,%.3e],[%.3e,%.3e,%.3e]] " "diffusion_coefficient=%.3e", p->id, p->chemistry_data.dmetal_mass_fraction_total, p->chemistry_data.dmetal_mass_fraction_from_SNIa, p->chemistry_data.dmetal_mass_fraction_from_AGB, p->chemistry_data.dmetal_mass_fraction_from_SNII, p->chemistry_data.diron_mass_fraction_from_SNIa, p->chemistry_data.shear_tensor[0][0], p->chemistry_data.shear_tensor[0][1], p->chemistry_data.shear_tensor[0][2], p->chemistry_data.shear_tensor[1][0], p->chemistry_data.shear_tensor[1][1], p->chemistry_data.shear_tensor[1][2], p->chemistry_data.shear_tensor[2][0], p->chemistry_data.shear_tensor[2][1], p->chemistry_data.shear_tensor[2][2], p->chemistry_data.diffusion_coefficient); for (int i = 0; i < chemistry_element_count; i++) { warning("[PID%lld] diffusion_rate[%i]=%.3e", p->id, i, p->chemistry_data.diffusion_rate[i]); } } #endif /* SWIFT_CHEMISTRY_COLIBRE_DEBUG_H */