void rho2_mix	(	double &	r,
		double	par[21],
		double &	res
	)

Definition at line 355 of file profiles.cc.

{
   //--- Sets res = Density square of the profile [Msol^2/kpc^6]
   // Input:
   //  r         Distance from the centre of the profile [kpc]
   //  par[0]    rho_1(r) normalisation: density [Msol/kpc^3] or proxy for dPdV [kpc^{-3}]
   //  par[1]    rho_1(r) scale radius [kpc]
   //  par[2]    rho_1(r) shape parameter #1
   //  par[3]    rho_1(r) shape parameter #2
   //  par[4]    rho_1(r) shape parameter #3
   //  par[5]    rho_1(r) card_profile [gENUM_PROFILE]
   //  par[6]    UNUSED
   //  par[7]    UNUSED
   //  par[8]    UNUSED
   //  par[9]    UNUSED
   //  par[10]   switch_rho - selects which combination of rho_1 and rho_2 to use
   //                        0 -> rho(r) = rho1(r)
   //                        1 -> rho(r) = rho1(r)-rho2(r)
   //                        2 -> rho(r) = rho1(r)*rho2(r)
   //  par[11]   rho_2(r) normalisation: density [Msol/kpc^3] or proxy for dPdV [kpc^{-3}]
   //  par[12]   rho_2(r) scale radius [kpc]
   //  par[13]   rho_2(r) shape parameter #1
   //  par[14]   rho_2(r) shape parameter #2
   //  par[15]   rho_2(r) shape parameter #3
   //  par[16]   rho_2(r) card_profile [gENUM_PROFILE]
   //  par[17]   UNUSED
   //  par[18]   UNUSED
   //  par[19]   UNUSED
   //  par[20]   UNUSED
   // Output:
   //  res       Value of rho^2(r)

   rho_mix(r, par, res);
   res *= res;
   return;
}