app_hor.h

/*
 *  Definition of Lorene class App_hor
 *
 */

/*
 *   Copyright (c) 2005 Lap-Ming Lin & Jerome Novak
 *
 *
 *   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
 *
 */

#ifndef __APP_HOR_H_ 
#define __APP_HOR_H_ 

/*
 * $Id: app_hor.h,v 1.7 2022/09/30 14:41:38 j_novak Exp $
 * $Log: app_hor.h,v $
 * Revision 1.7  2022/09/30 14:41:38  j_novak
 * Corrected comments.
 *
 * Revision 1.6  2014/10/13 08:52:31  j_novak
 * Lorene classes and functions now belong to the namespace Lorene.
 *
 * Revision 1.5  2012/01/02 13:52:57  j_novak
 * New parameter 'verbose' to get less output if needed.
 *
 * Revision 1.4  2005/12/09 09:35:59  lm_lin
 *
 * Minor fix in the documentation.
 *
 * Revision 1.3  2005/12/07 11:11:30  lm_lin
 *
 * Add option to turn off screen output during iterations.
 *
 * Revision 1.2  2005/11/17 14:19:49  lm_lin
 *
 * Check the expansion function evaluated on the apparent horizon after the
 * iteration of the 2-surface converges.
 *
 * Revision 1.1  2005/10/13 08:51:14  j_novak
 * New stuff for apparent horizon finder. For the moment, there is only an
 * external function. A class should come soon...
 *
 * $Header: /cvsroot/Lorene/C++/Include/app_hor.h,v 1.7 2022/09/30 14:41:38 j_novak Exp $
 *
 */


// Headers Lorene
#include "metric.h" 

namespace Lorene {
    // Function (Apparent horizon finder) 
    //-----------------------------------------------
    
bool ah_finder(const Metric& gamma, const Sym_tensor& k_dd_in, Valeur& h, Scalar& ex_fcn,
           double a_axis, double b_axis, double c_axis, bool verbose = true, 
           bool print = false, double precis = 1.e-8, double precis_exp = 1.e-6,  
           int it_max = 200, int it_relax = 200, double relax_fac = 1.) ;
        
}
#endif