#include <Metric.h>

Data Fields
int	opt

double	lab

double	L [1024]

double	S [1024]

Friends
double	LengthInterpole (const Metric &Ma, const Metric &Mb, R2 AB)

double	abscisseInterpole (const Metric &Ma, const Metric &Mb, R2, double s, int optim)

Detailed Description

Definition at line 85 of file Metric.h.

Friends And Related Function Documentation

◆ LengthInterpole

double LengthInterpole	(	const Metric &	Ma,
		const Metric &	Mb,
		R2	AB
	)

friend

Definition at line 158 of file Metric.cpp.

                                                                      {/*{{{*/
       /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Metric.cpp/LengthInterpole)*/
  
       double k=1./2.;
       int level=0;
       static int kkk=0;
       static  Metric Ms1[32],Ms2[32];
       static double lMs1[32],lMs2[32];
       static double K[32];
       double l=0,sss=0;
       Ms1[level]=Ma;
       Ms2[level]=Mb;
       double sa =  Ma.Length(AB.x,AB.y);
       double sb =  Mb.Length(AB.x,AB.y);
       lMs1[level]=sa;
       lMs2[level]=sb;
       K[level]=k;
       level++;
       int i=0;
       double * L= LastMetricInterpole.L, *S = LastMetricInterpole.S;
       double  sstop = 0.1; // Max(0.6,(sa+sb)/5000);
       while (level) {
          level--;
          Metric M1=Ms1[level];
          Metric M2=Ms2[level];
          k=K[level];
          double s1=  lMs1[level];
          double s2=  lMs2[level];
  
          double s= (s1+s2)*k;
          if( s > sstop   && level < 30 && i < 500-level ) {
             Metric Mi(0.5,M1,0.5,M2);
             double si = Mi.Length(AB.x,AB.y);
             if( Abs((s1+s2)-(si+si)) > s1*0.001) 
               {
                k=k/2;
                // we begin by the end to walk in the correct direction from a to b
                // due to the stack 
                Ms1[level]=Mi;
                Ms2[level]=M2;
                lMs1[level]=si;
                lMs2[level]=s2;
                K[level]=k;
                level++;
                Ms1[level]=M1;
                Ms2[level]=Mi;
                lMs1[level]=s1;
                lMs2[level]=si;
                K[level]=k;
                level++;
               }
             else
              L[i]= l += s,S[i]=sss+=k,i++;
          }
          else 
           L[i]= l += s,S[i]=sss+=k,i++;
       }
       // warning for optimisation S is in [0:0.5] not in [0:1]
       if (i>=512){
          _error_("i>=512");
       }
       LastMetricInterpole.lab=l;
       LastMetricInterpole.opt=i;
       if (i>200 && kkk++<10) _printf_("WARNING: LengthInterpole: ( i=" << i << " l=" << l << " sss=" << sss << " ) " << sstop << "\n"); 
       return l;
    }

◆ abscisseInterpole

double abscisseInterpole	(	const Metric &	Ma,
		const Metric &	Mb,
		R2	,
		double	s,
		int	optim
	)

friend

Definition at line 327 of file Metric.cpp.

                                                                                           { /*{{{*/
       /*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Metric.cpp/abscisseInterpole)*/
  
       if(!optim)  LengthInterpole(Ma,Mb,AB);
       double l  = s* LastMetricInterpole.lab,r;
       int j=LastMetricInterpole.opt-1;
  
       double * L= LastMetricInterpole.L, *S = LastMetricInterpole.S;
       // warning for optimisation S is the abcisse in [0:0.5]
       // and L is le lenght 
       if(l<=L[0]){
          r=2*S[0]*l/L[0];
       }
       else if (l>=L[j]){
          r=1;
       }
       else{
          int i=0;
          while (j-i>1){
             int k;
             k= (i+j)/2;
             if(l<=L[k]){
                j=k;// l<=L[j] 
             }
             else{
                i=k; //  L[i]<l
             }
          };
          if (i==j){
             r = 2*S[i];
          }
          else{
             r =  2*(S[i]*(L[j]-l)+ S[j]*(l-L[i]))/(L[j]-L[i]);
          }
       }
       if (r>1 || r<0){
          _error_("r>1 || r<0");
       }
       return r ;
    }

Field Documentation

◆ opt

int bamg::SaveMetricInterpole::opt

Definition at line 89 of file Metric.h.

◆ lab

double bamg::SaveMetricInterpole::lab

Definition at line 90 of file Metric.h.

◆ L

double bamg::SaveMetricInterpole::L[1024]

Definition at line 91 of file Metric.h.

◆ S

double bamg::SaveMetricInterpole::S[1024]