ahf2d_d.h

Go to the documentation of this file.

 
 #ifndef AHF2D_H
 #define AHF2D_H
 
 #include "cosmos.h" // Include if it needs access to Fmv0/Fmv grid data/methods
 #include <stdio.h>
 #include <math.h>
 // Include other necessary headers
 
 class Ahf2d {
 public:
     // --- Member Variables ---
     int ntheta;
     int nphi;
     double dth;
     double dphi;
     double hini;
     double var; // Corresponds to etaah in par_ahf.d
     double fac; // Corresponds to facn in par_ahf.d
     double err_p;
     double err_e;
     int ahloop;
 
     // Pointers to AHF internal grid arrays (Example for horizon shape 'h')
     double **h;
     double **ths;
     // Add other internal arrays if needed (e.g., for solver residuals, coefficients)
 
     Fmv0 *fm; // Use Fmv0* or Fmv* depending on what data is needed
 
     // --- Constructor & Destructor ---
     Ahf2d(int ntheta_in, int nphi_in, double hini_in, double var_in, double fac_in,
           double err_p_in, double err_e_in, int ahloop_in, Fmv0 *fm_in);
 
     ~Ahf2d();
 
     // --- Core Methods ---
     int find_ah(double *ah_mass = nullptr, double *ah_spin = nullptr); // Example output params
 
     void calc_expansion();
 
     int solve_elliptic(double **source, double **solution);
 
     double interpolate_to_ah(double ***f, double r, double th, double ph);
 
     double calculate_ah_mass();
 
     double calculate_ah_spin();
 
     void print_horizon_shape(const char* filename, double t);
 
     // --- Utility Functions ---
     double** matrix2d(long nrl, long nrh, long ncl, long nch);
     void free_matrix2d(double **m, long nrl, long nrh, long ncl, long nch);
 
 protected: // Or private
     // Add helper functions if needed, e.g., for calculating derivatives on the spherical grid
     double D1th(double **f, int i, int j);
     double D1ph(double **f, int i, int j);
     // Add second derivatives if needed (D2th, D2ph, Dthph)
 };
 
 
 #endif // AHF2D_H