[Commits] [svn:einsteintoolkit] GRHydro_InitData/trunk/src/ (Rev. 200)
rhaas at tapir.caltech.edu
rhaas at tapir.caltech.edu
Wed Mar 27 20:48:58 CDT 2013
User: rhaas
Date: 2013/03/27 08:48 PM
Modified:
/trunk/src/
GRHydro_BondiM.c
Log:
GRHydro_Init_Data: remove rho_guess from find_bondi_solution
was used to call the newton solver twice with the very same arguments.
Not sure what the idea was.
From: Roland Haas <rhaas at tapir.caltech.edu>
File Changes:
Directory: /trunk/src/
======================
File [modified]: GRHydro_BondiM.c
Delta lines: +1 -12
===================================================================
--- trunk/src/GRHydro_BondiM.c 2013-03-28 01:48:56 UTC (rev 199)
+++ trunk/src/GRHydro_BondiM.c 2013-03-28 01:48:57 UTC (rev 200)
@@ -647,7 +647,7 @@
const int ntries = 10000;
int itry;
- CCTK_REAL rhotmp, rho_guess;
+ CCTK_REAL rhotmp;
CCTK_REAL dr,ur;
@@ -662,17 +662,12 @@
}
else {
// rhotmp = (sqrt(Qdot) - 1.) * (gamma_eos - 1.) / ( gamma_eos * K );
- // rho_guess = pow( rhotmp , (1./(gamma_eos - 1.)) );
if(r < rs) { ur = pow(r,-0.5) ; }
else { ur = 0.5*pow(r,-1.5) ; }
*rho = Mdot / (4.*M_PI * r * r * ur);
}
}
- // safe guess value for multiple tries
- rho_guess = *rho;
-
-
// set global variables needed by residual function:
r_sol = r ;
@@ -681,12 +676,6 @@
retval = gnr_bondi( rho, NEWT_DIM_B, bondi_resid);
- // first try guess if failure
- if( retval ) {
- *rho = rho_guess;
- retval = gnr_bondi( rho, NEWT_DIM_B, bondi_resid);
- }
-
// If we were unsure about the guess and solver fails, then creep from known solution to desired point:
if( retval ) {
More information about the Commits
mailing list