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"); }