void RGRUN( ofstream& rgefile, MssmSoftsusy& m ,
double log10startscale, double log10endscale,
int number )
{
double scale ;
double step;
step = ( log10endscale - log10startscale) / number;
m.runto( pow( 10, log10startscale ) ) ;
for( int i = 0 ; i < number ; i++ ) {
scale = pow( 10, log10startscale + step*i) ;
m.runto( scale) ;
DoubleMatrix mQlmat = m.displaySoftMassSquared(mQl);
DoubleMatrix mUrmat = m.displaySoftMassSquared(mUr);
DoubleMatrix mDrmat = m.displaySoftMassSquared(mDr);
DoubleMatrix mLlmat = m.displaySoftMassSquared(mLl);
DoubleMatrix mErmat = m.displaySoftMassSquared(mEr);
double mHusqr = m.displayMh2Squared();
double mHdsqr = m.displayMh1Squared();
rgefile << scale
<< " " << m.displayGaugeCoupling(1)
<< " " << m.displayGaugeCoupling(2)
<< " " << m.displayGaugeCoupling(3)
<< " " << m.displayGaugino(1)
<< " " << m.displayGaugino(2)
<< " " << m.displayGaugino(3)
<< " " << mQlmat(3,3)
<< " " << mUrmat(3,3)
<< " " << mDrmat(3,3)
<< " " << mLlmat(3,3)
<< " " << mErmat(3,3)
<< " " << mQlmat(1,1)
<< " " << mUrmat(1,1)
<< " " << mDrmat(1,1)
<< " " << mLlmat(1,1)
<< " " << mErmat(1,1)
<< " " << mHusqr
<< " " << mHdsqr
<< endl ;
}
}
