virtual ::osiris::uint64 size( ) const {
::osiris::PythonState __pystate(getPythonThreadState());
if( ::osiris::PythonOverride func_size = this->get_override( "size" ) )
return func_size( );
else{
__pystate.leave();
return this->::osiris::File::size( );
}
}
/** @brief Calculates the size of the disk that contains the
* mass transfered to the bulge. */
double bulge_from_disk(double frac)
{
double x1,x2,x0,value;
/** @brief Calculates the size of the disk that contains the
* mass transfered to the bulge. The bulge is assumed
* to form with the same size. avoid doing "ln" from eq 35*/
x1=0.0;
x2=1.;
while ((func_size(x2,frac) * func_size(x1,frac))>0) {
x1=x2;
x2=x2*2;
}
x0=x1+(x2-x1)/2.;
value=func_size(x0,frac);
if (value < 0)
value = -value;
while(value>0.00001) {
if(func_size(x0,frac)*func_size(x2,frac)>0)
x2=x0;
else
x1=x0;
x0=x1+(x2-x1)/2.;
value=func_size(x0,frac);
if (value < 0)
value = -value;
}
return x0;
}
