30char valeur_equipot_out_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Valeur/valeur_equipot_out.C,v 1.5 2014/10/13 08:53:50 j_novak Exp $" ;
75#include "utilitaires.h"
79double valeur_equipot_fonc(
double,
const Param&) ;
88 int nz =
mg->get_nzone() ;
89 int nt =
mg->get_nt(0) ;
90 int np =
mg->get_np(0) ;
106 l_iso.set_etat_qcq() ;
111 parf.add_int_mod(
j, 0) ;
112 parf.add_int_mod(
k, 1) ;
123 for (
k=0;
k<np;
k++) {
125 for (
j=0;
j<nt;
j++) {
134 int nr =
mg->get_nr(
l) ;
136 for (
int i=0;
i<nr;
i++) {
137 double uux = (*this)(
l,
k,
j,
i) ;
151 "Valeur::equipot_outward: the point uu < uu0 has not been found"
153 cout <<
" for the phi index " <<
k
154 <<
" and the theta index " <<
j <<
endl ;
175 "Valeur::equipot_outward: WARNING: potentially discontinuous field !"
188 "Valeur::equipot_outward: the field has some negative value at the center !"
198 double x2 = (
mg->get_grille3d(
l2))->x[
i2] ;
199 double x3 = (
mg->get_grille3d(
l2))->x[
i3] ;
Time evolution with partial storage (*** under development ***).
Basic integer array class.
void equipot_outward(double uu0, int nz_search, double precis, int nitermax, int &niter, Itbl &l_iso, Tbl &xi_iso) const
Determines an equipotential surface of the field represented by *this (outward search).
const Mg3d * mg
Multi-grid Mgd3 on which this is defined.
Mtbl_cf * c_cf
Coefficients of the spectral expansion of the function.
void coef() const
Computes the coeffcients of *this.
int etat
Logical state (ETATNONDEF , ETATQCQ or ETATZERO ).
double zerosec(double(*f)(double, const Param &), const Param &par, double a, double b, double precis, int nitermax, int &niter, bool abort=true)
Finding the zero a function.