swap gsl_gamma_inc with gsl_gamma_inc_P for mass calculation

jamesgrimmett · jobovy · commit f2f76d51fab7 · 2025-01-09T13:10:39.000-05:00
diff --git a/galpy/potential/potential_c_ext/PowerSphericalPotentialwCutoff.c b/galpy/potential/potential_c_ext/PowerSphericalPotentialwCutoff.c
@@ -7,7 +7,7 @@
 //PowerSphericalPotentialwCutoff
 //3  arguments: amp, alpha, rc
 double mass(double r2,double alpha, double rc){
-  return 2. * M_PI * pow ( rc , 3. - alpha ) * ( gsl_sf_gamma ( 1.5 - 0.5 * alpha ) - gsl_sf_gamma_inc ( 1.5 - 0.5 * alpha , r2 / rc / rc ) );
+  return 2. * M_PI * pow ( rc , 3. - alpha ) * ( gsl_sf_gamma ( 1.5 - 0.5 * alpha ) * gsl_sf_gamma_inc_P ( 1.5 - 0.5 * alpha , r2 / rc / rc ) );
 }
 double PowerSphericalPotentialwCutoffEval(double R,double Z, double phi,
 					  double t,

Original file line number	Diff line number	Diff line change
`@@ -7,7 +7,7 @@`
`7`	`7`	`//PowerSphericalPotentialwCutoff`
`8`	`8`	`//3 arguments: amp, alpha, rc`
`9`	`9`	`double mass(double r2,double alpha, double rc){`
`10`		`- return 2. * M_PI * pow ( rc , 3. - alpha ) * ( gsl_sf_gamma ( 1.5 - 0.5 * alpha ) - gsl_sf_gamma_inc ( 1.5 - 0.5 * alpha , r2 / rc / rc ) );`
	`10`	`+ return 2. * M_PI * pow ( rc , 3. - alpha ) * ( gsl_sf_gamma ( 1.5 - 0.5 * alpha ) * gsl_sf_gamma_inc_P ( 1.5 - 0.5 * alpha , r2 / rc / rc ) );`
`11`	`11`	`}`
`12`	`12`	`double PowerSphericalPotentialwCutoffEval(double R,double Z, double phi,`
`13`	`13`	`double t,`