strot_diff_global.C

/*
 *  Methods for computing global quantities within the class Star_rot_diff
 *
 *    (see file star_rot_diff.h for documentation).
 *
 */

/*
 *   Copyright (c) 2025 Santiago Jaraba
 *
 *   This file is part of LORENE.
 *
 *   LORENE is free software; you can redistribute it and/or modify
 *   it under the terms of the GNU General Public License version 2
 *   as published by the Free Software Foundation.
 *
 *   LORENE is distributed in the hope that it will be useful,
 *   but WITHOUT ANY WARRANTY; without even the implied warranty of
 *   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 *   GNU General Public License for more details.
 *
 *   You should have received a copy of the GNU General Public License
 *   along with LORENE; if not, write to the Free Software
 *   Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
 *
 */


// C headers
#include <cmath>

// Lorene headers
#include "star_rot_diff.h"


                     //---------------------//
                     //        T/W          //
                     //---------------------//

namespace Lorene {
double Star_rot_diff::tsw() const {

    if (p_tsw == 0x0) {    // a new computation is required
    
    Scalar dens = uuu ; 

    dens.mult_rsint() ;         //  Multiplication by r sin(theta)

    if (relativistic) {
        dens = omega_field * a_car * b_car * (ener_euler + press) 
            * dens ; 
    }
    else {    // Newtonian case 
        dens = omega_field * nbar * dens ; 
    }

    dens.std_spectral_base() ; 

    double tcin = 0.5 * dens.integrale() ;
    
    if (relativistic) {
        
        Scalar dens2 = a_car * bbb * gam_euler * ener ;
        dens2.std_spectral_base() ; 
        double mass_p = dens2.integrale() ; 
        
        p_tsw = new double( tcin / ( mass_p + tcin - mass_g() ) ) ;     
       
    }
    else {      // Newtonian case 
        Scalar dens2 = 0.5 * nbar * logn ;
        dens2.std_spectral_base() ; 
        double wgrav = dens2.integrale() ; 
        p_tsw = new double( tcin / fabs(wgrav) ) ;  
    }


    }
    
    return *p_tsw ; 

}

}