Commit 546664ac authored by Maude Le Jeune's avatar Maude Le Jeune
Browse files

No commit message

No commit message
parent a47ade7a
/*
* Project: SZ
* Date: 08/10
* Authors: M. Le Jeune
*
*/
/* Copyright (C) 2008 APC CNRS Université Paris Diderot <lejeune@apc.univ-paris7.fr>
*
* This program is free software; you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation; either version 3 of the License, or
* (at your option) any later version.
* This program 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 this program; if not, see http://www.gnu.org/licenses/gpl.html
*/
#include <cmath>
#include <iostream>
#include <string>
//#include "qromb.h"
#include <sstream>
#include "arr.h"
#include "vec3.h"
#include <gsl/gsl_integration.h>
using namespace std;
/// \param theta_min : min bound for integration (max is pi/2)
//double theta_min;
/// \param beta : beta profile parameter
//double beta = (double)2/3;
/// \param rmax : viriel radius
//double rmax;
/// \param theta_c : core radius
//double theta_c;
/// beta function to be integrated between theta_min and pi/2
/// \param theta : angular distance from cluster center
double betat (double theta, void * params);
/// beta function at center to be integrated between 0 and rmax
/// \param theta : angular distance from cluster center
double beta0 (double theta, void * params);
/// give the beta profile value.
/// \param theta_i : angular distance from cluster center
/// \param rmax : viriel radius
/// \param r_core : core radius
double beta_profile (double rmax, double r_core, double theta_i,double beta);
/// nwf function to be integrated between theta_min and pi/2
/// \param theta : angular distance from cluster center
double nwft (double theta, void * params);
/// nwf function at center to be integrated between 0 and rmax
/// \param theta : angular distance from cluster center
double nwf0 (double theta, void * params);
double xmm_profile (double rmax, double scale_radius, double theta_i, double alpha, double beta, double gamma);
/// nwf function at center to be integrated between 0 and rmax
/// \param theta : angular distance from cluster center
double nwf5 (double theta, void * params);
double y5r500 (double scale_radius, double alpha, double beta, double gamma);
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment