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