TEST_F(TestGaussianPot, OneGaussWorksNonZero)
{
energy_test({0, 0, 0, 0}, -1);
energy_test({1, 1, 1, 1}, -0.1353352832366127);
energy_test({2, 2, 2, 2}, -0.0003354626279025118);
energy_test({1, 2, 3, 4}, -3.059023205018258e-7);
gradient_test({0, 0, 0, 0}, {0, 0, 0, 0});
gradient_test({0, 0.1, 0.2, 0.3}, {0, 0.09323938199059482, 0.1864787639811896, 0.2797181459717845});
gradient_test({1, 1, 1, 1}, {0.1353352832366127, 0.1353352832366127, 0.1353352832366127, 0.1353352832366127});
}
void star1d::dump_info() {
DEBUG_FUNCNAME;
printf("ESTER 1d model file");
printf(" (Version %s)\n", version.name.c_str());
// printf("\n1d ESTER model file (Version %d.%d rev %d",version.major,version.minor,version.rev);
// if(version.svn) printf(".svn");
// printf(")\n\n");
printf("General parameters:\n\n");
printf("\tMass = %.5f Msun (%e g)\n",M/M_SUN,M);
printf("\tRadius = %.5f Rsun (%e cm)\n",R/R_SUN,R);
printf("\tLuminosity = %.4f Lsun (%e erg/s)\n",luminosity()/L_SUN,luminosity());
printf("\tTeff = %.2f\n",Teff()(0));
printf("\tlog(geff) = %.4f\n",log10(gsup())(0));
printf("\tX0=%.4f Y0=%.4f Z0=%.4f\n",X0,Y0,Z0);
printf("\n");
if(conv==0) printf("No convective core\n\n");
else {
printf("Convective core:\n\n");
double mcc=Mcore();
printf("\tMass_core = %.5f Msun (%e g)\n",mcc/M_SUN,mcc);
double rcc=Rcore()(0);
printf("\tRadius_core (p) = %.5f Rsun (%e cm)\n",rcc/R_SUN,rcc);
printf("\tX_core/X_env = %.4f\n",Xc);
printf("\n");
}
printf("Central values:\n\n");
printf("\tTemperature = %e K\n",Tc);
printf("\tDensity = %e g/cm3\n",rhoc);
printf("\tPressure = %e dyn/cm2\n",pc);
printf("\n");
printf("Grid parameters:\n\n");
printf("\t # of domains = %d\n",ndomains);
printf("\t # of domains in convective core = %d\n",conv);
printf("\t nr = %d (",nr);
for(int n=0;n<ndomains;n++) {
printf("%d",map.gl.npts[n]);
if(n<ndomains-1) printf(",");
}
printf(")\n");
printf("\n");
printf("Additional parameters:\n\n");
printf("\tOpacity = %s\n",opa.name);
printf("\tEquation of state = %s\n",eos.name);
printf("\tNuclear reactions = %s\n",nuc.name);
printf("\tAtmosphere = %s\n",atm.name);
printf("\tsurff = %e\n",surff);
printf("\tcore_convec = %d\n",core_convec);
printf("\tenv_convec = %d\n",core_convec);
printf("\tmin_core_size = %e\n",min_core_size);
printf("\n");
printf("Tests:\n\n");
printf("\tVirial test = %e\n",virial());
printf("\tEnergy test = %e\n",energy_test());
printf("\n");
}
