//----------------------------------------------------
// project: floater_delta     header generated by YAO version v.9
//----------------------------------------------------

// ¤ ¤ ¤ ¤ ¤ ¤ ¤ ¤ LES FONCTIONS PREDEFINIES
/*M1_TBTOGT  */ 
void YtabTOgrad_r_float(int frompdt, int topdt, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_r_float; ++Yws){
				 YG_r_float(Yws, Yw1, Ywt) = tab[Y3windice]; ++Y3windice;}}}
}/*M1_GEQPDT*/ 
void YgradEQPdelta_r_float(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_r_float; ++Yws){
				 YG_r_float(Yws, Yw1, Ywt) = YD_r_float(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_DEQPST*/ 
void YdeltaEQPstate_r_float(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_r_float; ++Yws){
				 YD_r_float(Yws, Yw1, Ywt) = YS_r_float(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_DEQPGT*/ 
void YdeltaEQPgrad_r_float(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_r_float; ++Yws){
				 YD_r_float(Yws, Yw1, Ywt) = YG_r_float(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_SEQPDT*/ 
void YstateEQPdelta_r_float(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_r_float; ++Yws){
				 YS_r_float(Yws, Yw1, Ywt) = YD_r_float(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_SEQPODT*/ 
void YstateEQPOdelta_r_float(int frompdt, int topdt, double pfact, char *codop)
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_r_float; ++Yws){
				 YS_r_float(Yws, Yw1, Ywt) = YD_r_float(Yws, Yw1, Ywt) * (pfact*(1+Ygetval(codop)));}}}
}/*M1_SEQAPTDT*/ 
void YstateEQAPTdelta_r_float(int frompdt, int topdt, double pfact, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_r_float; ++Yws){
				 YS_r_float(Yws, Yw1, Ywt) = YD_r_float(Yws, Yw1, Ywt) + pfact * tab[Y3windice++];}}}
}/*M1_GEQPST*/ 
void YgradEQPstate_r_float(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_r_float; ++Yws){
				 YG_r_float(Yws, Yw1, Ywt) = YS_r_float(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_GEQPOST*/ 
void YgradEQPOstate_r_float(int frompdt, int topdt, double pfact, char *codop)
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_r_float; ++Yws){
				 YG_r_float(Yws, Yw1, Ywt) = YS_r_float(Yws, Yw1, Ywt) * (pfact*(1+Ygetval(codop)));}}}
}/*M1_STOTBT  */ 
void YstateTOtab_r_float(int frompdt, int topdt, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_r_float; ++Yws){
				 tab[Y3windice] = YS_r_float(Yws, Yw1, Ywt); ++Y3windice;}}}
}/*M1_GTOTBT  */ 
void YgradTOtab_r_float(int frompdt, int topdt, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_r_float; ++Yws){
				 tab[Y3windice] = YG_r_float(Yws, Yw1, Ywt); ++Y3windice;}}}
}/*M1_RGRADT*/ 
void Yrazgrad_r_float()
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=0; Ywt<YNBALLTIME_T_float; ++Ywt){
			 for (int Yws=0; Yws<YNBS_r_float; ++Yws){
				 Yr_float(Yw1)->Ygrad[Ywt][Yws] = 0.0;}}}
}/*M1_SETAT*/ 
void Ysetstate_r_float(double val)
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=0; Ywt<YNBALLTIME_T_float;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_r_float; ++Yws){
				 Yr_float(Yw1)->Ystate[Ywt][Yws] = val;}}}
}/*M1_YIOUT*/ 
void Yio_r_float(int Yws, int Yw1, int Yw2, int Yw3, int Ywt, double val)	
{  Yiorwf ("r_float", Yws, Yw1, Yw2, Yw3, Ywt, &YS_r_float(Yws, Yw1, Ywt), val);
}/*M00_GTEST */
int Ytestdf_r_float (int modop, int KeKo, float pdx, float ptol, int ynbi, int ynbs, char *nmmod, YREAL yinput[], YREAL ystate[], r_float *Modul)
{	 int wk, wj, nbKo=0; YREAL wSj[YMAX_NBS],wdxk,wJjk,wyj,wDjk,wdifjk;
	 YREAL svx; int resKo; char okko[2][3]={"ok", "KO"}; 
	 if (!Ydftestijkt(Yid_r_float)) return(0); 
	 memset (Yjac, 0, YMAX_JAC_NBI*YMAX_JAC_NBS*sizeof(YREAL));
/*M01_GTESTT*/
	 Modul->forward( Yting[0], Yting[1], Yting[2], Yting[3] );
	 for (wj=0; wj<ynbs; ++wj) wSj[wj] = ystate[YTemps*ynbs + wj];
	 Modul->backward( Yting[0], Yting[1], Yting[2], Yting[3] );
	 for (wk=0; wk<ynbi; ++wk)
	 {	 svx=yinput[wk];
			 wdxk = svx*pdx;
			 if (modop||svx==0.0) wdxk = pdx; 
			 yinput[wk] += wdxk;
			 Modul->forward( Yting[0], Yting[1], Yting[2], Yting[3] );
			 for (wj=0; wj<ynbs; ++wj)
			 { wJjk = Yjac[wj][wk];
				 wyj = ystate[YTemps*ynbs + wj];
				 wDjk = (wyj - wSj[wj]) / wdxk;
				 wdifjk = fabs(wDjk - wJjk);
				 resKo=(wdifjk/fabs(wJjk))>ptol;  nbKo+=resKo;
				 if(resKo || !KeKo) 
				 printf(" %7s | J[%2i,%2i] |% 12.6e |% 12.6e |% 12.6e :  %s |% 12.6e | (%i %i %i %i)\n",
								 nmmod, wj+1, wk+1, wJjk, wDjk, wdifjk, okko[resKo], svx, Yi+1, Yj+1, Yk+1, Yt+1);
				} 
			  yinput[wk] = svx;
	 }
	 return(nbKo);
}/*M1_CREU*/ 
void Ycreate_r_float()  
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 if ((Yr_float(Yw1) = new r_float ())==NULL) Ycreate_errnew();}
}/*M1_TBTOGT  */ 
void YtabTOgrad_ur(int frompdt, int topdt, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_ur; ++Yws){
				 YG_ur(Yws, Yw1, Ywt) = tab[Y3windice]; ++Y3windice;}}}
}/*M1_GEQPDT*/ 
void YgradEQPdelta_ur(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_ur; ++Yws){
				 YG_ur(Yws, Yw1, Ywt) = YD_ur(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_DEQPST*/ 
void YdeltaEQPstate_ur(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_ur; ++Yws){
				 YD_ur(Yws, Yw1, Ywt) = YS_ur(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_DEQPGT*/ 
void YdeltaEQPgrad_ur(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_ur; ++Yws){
				 YD_ur(Yws, Yw1, Ywt) = YG_ur(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_SEQPDT*/ 
void YstateEQPdelta_ur(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_ur; ++Yws){
				 YS_ur(Yws, Yw1, Ywt) = YD_ur(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_SEQPODT*/ 
void YstateEQPOdelta_ur(int frompdt, int topdt, double pfact, char *codop)
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_ur; ++Yws){
				 YS_ur(Yws, Yw1, Ywt) = YD_ur(Yws, Yw1, Ywt) * (pfact*(1+Ygetval(codop)));}}}
}/*M1_SEQAPTDT*/ 
void YstateEQAPTdelta_ur(int frompdt, int topdt, double pfact, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_ur; ++Yws){
				 YS_ur(Yws, Yw1, Ywt) = YD_ur(Yws, Yw1, Ywt) + pfact * tab[Y3windice++];}}}
}/*M1_GEQPST*/ 
void YgradEQPstate_ur(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_ur; ++Yws){
				 YG_ur(Yws, Yw1, Ywt) = YS_ur(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_GEQPOST*/ 
void YgradEQPOstate_ur(int frompdt, int topdt, double pfact, char *codop)
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_ur; ++Yws){
				 YG_ur(Yws, Yw1, Ywt) = YS_ur(Yws, Yw1, Ywt) * (pfact*(1+Ygetval(codop)));}}}
}/*M1_STOTBT  */ 
void YstateTOtab_ur(int frompdt, int topdt, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_ur; ++Yws){
				 tab[Y3windice] = YS_ur(Yws, Yw1, Ywt); ++Y3windice;}}}
}/*M1_GTOTBT  */ 
void YgradTOtab_ur(int frompdt, int topdt, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_ur; ++Yws){
				 tab[Y3windice] = YG_ur(Yws, Yw1, Ywt); ++Y3windice;}}}
}/*M1_RGRADT*/ 
void Yrazgrad_ur()
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=0; Ywt<YNBALLTIME_T_float; ++Ywt){
			 for (int Yws=0; Yws<YNBS_ur; ++Yws){
				 Yur(Yw1)->Ygrad[Ywt][Yws] = 0.0;}}}
}/*M1_SETAT*/ 
void Ysetstate_ur(double val)
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=0; Ywt<YNBALLTIME_T_float;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_ur; ++Yws){
				 Yur(Yw1)->Ystate[Ywt][Yws] = val;}}}
}/*M1_YIOUT*/ 
void Yio_ur(int Yws, int Yw1, int Yw2, int Yw3, int Ywt, double val)	
{  Yiorwf ("ur", Yws, Yw1, Yw2, Yw3, Ywt, &YS_ur(Yws, Yw1, Ywt), val);
}/*M00_GTEST */
int Ytestdf_ur (int modop, int KeKo, float pdx, float ptol, int ynbi, int ynbs, char *nmmod, YREAL yinput[], YREAL ystate[], ur *Modul)
{	 int wk, wj, nbKo=0; YREAL wSj[YMAX_NBS],wdxk,wJjk,wyj,wDjk,wdifjk;
	 YREAL svx; int resKo; char okko[2][3]={"ok", "KO"}; 
	 if (!Ydftestijkt(Yid_ur)) return(0); 
	 memset (Yjac, 0, YMAX_JAC_NBI*YMAX_JAC_NBS*sizeof(YREAL));
/*M01_GTESTT*/
	 Modul->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4], Yting[5], Yting[6], Yting[7], Yting[8], Yting[9], Yting[10], Yting[11], Yting[12], Yting[13], Yting[14], Yting[15], Yting[16], Yting[17] );
	 for (wj=0; wj<ynbs; ++wj) wSj[wj] = ystate[YTemps*ynbs + wj];
	 Modul->backward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4], Yting[5], Yting[6], Yting[7], Yting[8], Yting[9], Yting[10], Yting[11], Yting[12], Yting[13], Yting[14], Yting[15], Yting[16], Yting[17] );
	 for (wk=0; wk<ynbi; ++wk)
	 {	 svx=yinput[wk];
			 wdxk = svx*pdx;
			 if (modop||svx==0.0) wdxk = pdx; 
			 yinput[wk] += wdxk;
			 Modul->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4], Yting[5], Yting[6], Yting[7], Yting[8], Yting[9], Yting[10], Yting[11], Yting[12], Yting[13], Yting[14], Yting[15], Yting[16], Yting[17] );
			 for (wj=0; wj<ynbs; ++wj)
			 { wJjk = Yjac[wj][wk];
				 wyj = ystate[YTemps*ynbs + wj];
				 wDjk = (wyj - wSj[wj]) / wdxk;
				 wdifjk = fabs(wDjk - wJjk);
				 resKo=(wdifjk/fabs(wJjk))>ptol;  nbKo+=resKo;
				 if(resKo || !KeKo) 
				 printf(" %7s | J[%2i,%2i] |% 12.6e |% 12.6e |% 12.6e :  %s |% 12.6e | (%i %i %i %i)\n",
								 nmmod, wj+1, wk+1, wJjk, wDjk, wdifjk, okko[resKo], svx, Yi+1, Yj+1, Yk+1, Yt+1);
				} 
			  yinput[wk] = svx;
	 }
	 return(nbKo);
}/*M1_CREU*/ 
void Ycreate_ur()  
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 if ((Yur(Yw1) = new ur ())==NULL) Ycreate_errnew();}
}/*M1_TBTOGT  */ 
void YtabTOgrad_uinter(int frompdt, int topdt, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_uinter; ++Yws){
				 YG_uinter(Yws, Yw1, Ywt) = tab[Y3windice]; ++Y3windice;}}}
}/*M1_GEQPDT*/ 
void YgradEQPdelta_uinter(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_uinter; ++Yws){
				 YG_uinter(Yws, Yw1, Ywt) = YD_uinter(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_DEQPST*/ 
void YdeltaEQPstate_uinter(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_uinter; ++Yws){
				 YD_uinter(Yws, Yw1, Ywt) = YS_uinter(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_DEQPGT*/ 
void YdeltaEQPgrad_uinter(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_uinter; ++Yws){
				 YD_uinter(Yws, Yw1, Ywt) = YG_uinter(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_SEQPDT*/ 
void YstateEQPdelta_uinter(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_uinter; ++Yws){
				 YS_uinter(Yws, Yw1, Ywt) = YD_uinter(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_SEQPODT*/ 
void YstateEQPOdelta_uinter(int frompdt, int topdt, double pfact, char *codop)
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_uinter; ++Yws){
				 YS_uinter(Yws, Yw1, Ywt) = YD_uinter(Yws, Yw1, Ywt) * (pfact*(1+Ygetval(codop)));}}}
}/*M1_SEQAPTDT*/ 
void YstateEQAPTdelta_uinter(int frompdt, int topdt, double pfact, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_uinter; ++Yws){
				 YS_uinter(Yws, Yw1, Ywt) = YD_uinter(Yws, Yw1, Ywt) + pfact * tab[Y3windice++];}}}
}/*M1_GEQPST*/ 
void YgradEQPstate_uinter(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_uinter; ++Yws){
				 YG_uinter(Yws, Yw1, Ywt) = YS_uinter(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_GEQPOST*/ 
void YgradEQPOstate_uinter(int frompdt, int topdt, double pfact, char *codop)
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_uinter; ++Yws){
				 YG_uinter(Yws, Yw1, Ywt) = YS_uinter(Yws, Yw1, Ywt) * (pfact*(1+Ygetval(codop)));}}}
}/*M1_STOTBT  */ 
void YstateTOtab_uinter(int frompdt, int topdt, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_uinter; ++Yws){
				 tab[Y3windice] = YS_uinter(Yws, Yw1, Ywt); ++Y3windice;}}}
}/*M1_GTOTBT  */ 
void YgradTOtab_uinter(int frompdt, int topdt, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_uinter; ++Yws){
				 tab[Y3windice] = YG_uinter(Yws, Yw1, Ywt); ++Y3windice;}}}
}/*M1_RGRADT*/ 
void Yrazgrad_uinter()
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=0; Ywt<YNBALLTIME_T_float; ++Ywt){
			 for (int Yws=0; Yws<YNBS_uinter; ++Yws){
				 Yuinter(Yw1)->Ygrad[Ywt][Yws] = 0.0;}}}
}/*M1_SETAT*/ 
void Ysetstate_uinter(double val)
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=0; Ywt<YNBALLTIME_T_float;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_uinter; ++Yws){
				 Yuinter(Yw1)->Ystate[Ywt][Yws] = val;}}}
}/*M1_YIOUT*/ 
void Yio_uinter(int Yws, int Yw1, int Yw2, int Yw3, int Ywt, double val)	
{  Yiorwf ("uinter", Yws, Yw1, Yw2, Yw3, Ywt, &YS_uinter(Yws, Yw1, Ywt), val);
}/*M1_CREY*/ 
void Ycreate_uinter() 	
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 if ((Yuinter(Yw1) = new Yaouinter ())==NULL) Ycreate_errnew();}
}/*M1_TBTOGT  */ 
void YtabTOgrad_xinter(int frompdt, int topdt, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_xinter; ++Yws){
				 YG_xinter(Yws, Yw1, Ywt) = tab[Y3windice]; ++Y3windice;}}}
}/*M1_GEQPDT*/ 
void YgradEQPdelta_xinter(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_xinter; ++Yws){
				 YG_xinter(Yws, Yw1, Ywt) = YD_xinter(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_DEQPST*/ 
void YdeltaEQPstate_xinter(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_xinter; ++Yws){
				 YD_xinter(Yws, Yw1, Ywt) = YS_xinter(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_DEQPGT*/ 
void YdeltaEQPgrad_xinter(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_xinter; ++Yws){
				 YD_xinter(Yws, Yw1, Ywt) = YG_xinter(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_SEQPDT*/ 
void YstateEQPdelta_xinter(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_xinter; ++Yws){
				 YS_xinter(Yws, Yw1, Ywt) = YD_xinter(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_SEQPODT*/ 
void YstateEQPOdelta_xinter(int frompdt, int topdt, double pfact, char *codop)
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_xinter; ++Yws){
				 YS_xinter(Yws, Yw1, Ywt) = YD_xinter(Yws, Yw1, Ywt) * (pfact*(1+Ygetval(codop)));}}}
}/*M1_SEQAPTDT*/ 
void YstateEQAPTdelta_xinter(int frompdt, int topdt, double pfact, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_xinter; ++Yws){
				 YS_xinter(Yws, Yw1, Ywt) = YD_xinter(Yws, Yw1, Ywt) + pfact * tab[Y3windice++];}}}
}/*M1_GEQPST*/ 
void YgradEQPstate_xinter(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_xinter; ++Yws){
				 YG_xinter(Yws, Yw1, Ywt) = YS_xinter(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_GEQPOST*/ 
void YgradEQPOstate_xinter(int frompdt, int topdt, double pfact, char *codop)
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_xinter; ++Yws){
				 YG_xinter(Yws, Yw1, Ywt) = YS_xinter(Yws, Yw1, Ywt) * (pfact*(1+Ygetval(codop)));}}}
}/*M1_STOTBT  */ 
void YstateTOtab_xinter(int frompdt, int topdt, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_xinter; ++Yws){
				 tab[Y3windice] = YS_xinter(Yws, Yw1, Ywt); ++Y3windice;}}}
}/*M1_GTOTBT  */ 
void YgradTOtab_xinter(int frompdt, int topdt, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_xinter; ++Yws){
				 tab[Y3windice] = YG_xinter(Yws, Yw1, Ywt); ++Y3windice;}}}
}/*M1_RGRADT*/ 
void Yrazgrad_xinter()
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=0; Ywt<YNBALLTIME_T_float; ++Ywt){
			 for (int Yws=0; Yws<YNBS_xinter; ++Yws){
				 Yxinter(Yw1)->Ygrad[Ywt][Yws] = 0.0;}}}
}/*M1_SETAT*/ 
void Ysetstate_xinter(double val)
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 for (int Ywt=0; Ywt<YNBALLTIME_T_float;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_xinter; ++Yws){
				 Yxinter(Yw1)->Ystate[Ywt][Yws] = val;}}}
}/*M1_YIOUT*/ 
void Yio_xinter(int Yws, int Yw1, int Yw2, int Yw3, int Ywt, double val)	
{  Yiorwf ("xinter", Yws, Yw1, Yw2, Yw3, Ywt, &YS_xinter(Yws, Yw1, Ywt), val);
}/*M1_CREY*/ 
void Ycreate_xinter() 	
{
	 for (int Yw1=0; Yw1<YA1_S_float; ++Yw1){
		 if ((Yxinter(Yw1) = new Yaoxinter ())==NULL) Ycreate_errnew();}
}/*M2_TBTOG   */ 
void YtabTOgrad_u(YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_euler; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler;  ++Yw2){
			 for (int Yws=0; Yws<YNBS_u; ++Yws){
				 YG_u(Yws, Yw1, Yw2) = tab[Y3windice]; ++Y3windice;}}}
}/*M2_GEQPD */ 
void YgradEQPdelta_u(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_euler; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler;  ++Yw2){
			 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_u; ++Yws){
				 YG_u(Yws, Yw1, Yw2) = YD_u(Yws, Yw1, Yw2, Ywt) * pfact;}}}}
}/*M2_DEQPS */ 
void YdeltaEQPstate_u(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_euler; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler;  ++Yw2){
			 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_u; ++Yws){
				 YD_u(Yws, Yw1, Yw2, Ywt) = YS_u(Yws, Yw1, Yw2) * pfact;}}}}
}/*M2_DEQPG */ 
void YdeltaEQPgrad_u(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_euler; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler;  ++Yw2){
			 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_u; ++Yws){
				 YD_u(Yws, Yw1, Yw2, Ywt) = YG_u(Yws, Yw1, Yw2) * pfact;}}}}
}/*M2_SEQPD */ 
void YstateEQPdelta_u(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_euler; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler;  ++Yw2){
			 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_u; ++Yws){
				 YS_u(Yws, Yw1, Yw2) = YD_u(Yws, Yw1, Yw2, Ywt) * pfact;}}}}
}/*M2_SEQPOD */ 
void YstateEQPOdelta_u(int frompdt, int topdt, double pfact, char *codop)
{
	 for (int Yw1=0; Yw1<YA1_S_euler; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler;  ++Yw2){
			 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_u; ++Yws){
				 YS_u(Yws, Yw1, Yw2) = YD_u(Yws, Yw1, Yw2, Ywt) * (pfact*(1+Ygetval(codop)));}}}}
}/*M2_SEQAPTD */ 
void YstateEQAPTdelta_u(int frompdt, int topdt, double pfact, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_euler; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler;  ++Yw2){
			 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_u; ++Yws){
				 YS_u(Yws, Yw1, Yw2) = YD_u(Yws, Yw1, Yw2, Ywt) + pfact * tab[Y3windice++];}}}}
}/*M2_GEQPS */ 
void YgradEQPstate_u(double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_euler; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler;  ++Yw2){
			 for (int Yws=0; Yws<YNBS_u; ++Yws){
				 YG_u(Yws, Yw1, Yw2) = YS_u(Yws, Yw1, Yw2) * pfact;}}}
}/*M2_GEQPOS */ 
void YgradEQPOstate_u(double pfact, char *codop)
{
	 for (int Yw1=0; Yw1<YA1_S_euler; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler;  ++Yw2){
			 for (int Yws=0; Yws<YNBS_u; ++Yws){
				 YG_u(Yws, Yw1, Yw2) = YS_u(Yws, Yw1, Yw2) * (pfact*(1+Ygetval(codop)));}}}
}/*M2_STOTB   */ 
void YstateTOtab_u(YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_euler; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler;  ++Yw2){
			 for (int Yws=0; Yws<YNBS_u; ++Yws){
				 tab[Y3windice] = YS_u(Yws, Yw1, Yw2); ++Y3windice;}}}
}/*M2_GTOTB   */ 
void YgradTOtab_u(YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_euler; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler;  ++Yw2){
			 for (int Yws=0; Yws<YNBS_u; ++Yws){
				 tab[Y3windice] = YG_u(Yws, Yw1, Yw2); ++Y3windice;}}}
}/*M2_RGRAD*/ 
void Yrazgrad_u()
{
	 for (int Yw1=0; Yw1<YA1_S_euler; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler; ++Yw2){
			 for (int Yws=0; Yws<YNBS_u; ++Yws){
				 Yu(Yw1, Yw2)->Ygrad[Yws] = 0.0;}}}
}/*M2_SETA*/ 
void Ysetstate_u(double val)
{
	 for (int Yw1=0; Yw1<YA1_S_euler; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler;  ++Yw2){
			 for (int Yws=0; Yws<YNBS_u; ++Yws){
				 Yu(Yw1, Yw2)->Ystate[Yws] = val;}}}
}/*M2_YIOTU*/ 
void Yio_u(int Yws, int Yw1, int Yw2, int Yw3, int Ywt, double val)	
{  Yiorwf ("u", Yws, Yw1, Yw2, Yw3, Ywt, &YS_u(Yws, Yw1, Yw2), val);
}/*M2_CREY*/ 
void Ycreate_u()  
{
	 for (int Yw1=0; Yw1<YA1_S_euler; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler; ++Yw2){
			 if ((Yu(Yw1, Yw2)= new Yaou ())==NULL) Ycreate_errnew();}}
}/*M2_TBTOG   */ 
void YtabTOgrad_v(YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_euler; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler;  ++Yw2){
			 for (int Yws=0; Yws<YNBS_v; ++Yws){
				 YG_v(Yws, Yw1, Yw2) = tab[Y3windice]; ++Y3windice;}}}
}/*M2_GEQPD */ 
void YgradEQPdelta_v(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_euler; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler;  ++Yw2){
			 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_v; ++Yws){
				 YG_v(Yws, Yw1, Yw2) = YD_v(Yws, Yw1, Yw2, Ywt) * pfact;}}}}
}/*M2_DEQPS */ 
void YdeltaEQPstate_v(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_euler; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler;  ++Yw2){
			 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_v; ++Yws){
				 YD_v(Yws, Yw1, Yw2, Ywt) = YS_v(Yws, Yw1, Yw2) * pfact;}}}}
}/*M2_DEQPG */ 
void YdeltaEQPgrad_v(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_euler; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler;  ++Yw2){
			 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_v; ++Yws){
				 YD_v(Yws, Yw1, Yw2, Ywt) = YG_v(Yws, Yw1, Yw2) * pfact;}}}}
}/*M2_SEQPD */ 
void YstateEQPdelta_v(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_euler; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler;  ++Yw2){
			 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_v; ++Yws){
				 YS_v(Yws, Yw1, Yw2) = YD_v(Yws, Yw1, Yw2, Ywt) * pfact;}}}}
}/*M2_SEQPOD */ 
void YstateEQPOdelta_v(int frompdt, int topdt, double pfact, char *codop)
{
	 for (int Yw1=0; Yw1<YA1_S_euler; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler;  ++Yw2){
			 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_v; ++Yws){
				 YS_v(Yws, Yw1, Yw2) = YD_v(Yws, Yw1, Yw2, Ywt) * (pfact*(1+Ygetval(codop)));}}}}
}/*M2_SEQAPTD */ 
void YstateEQAPTdelta_v(int frompdt, int topdt, double pfact, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_euler; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler;  ++Yw2){
			 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_v; ++Yws){
				 YS_v(Yws, Yw1, Yw2) = YD_v(Yws, Yw1, Yw2, Ywt) + pfact * tab[Y3windice++];}}}}
}/*M2_GEQPS */ 
void YgradEQPstate_v(double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_euler; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler;  ++Yw2){
			 for (int Yws=0; Yws<YNBS_v; ++Yws){
				 YG_v(Yws, Yw1, Yw2) = YS_v(Yws, Yw1, Yw2) * pfact;}}}
}/*M2_GEQPOS */ 
void YgradEQPOstate_v(double pfact, char *codop)
{
	 for (int Yw1=0; Yw1<YA1_S_euler; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler;  ++Yw2){
			 for (int Yws=0; Yws<YNBS_v; ++Yws){
				 YG_v(Yws, Yw1, Yw2) = YS_v(Yws, Yw1, Yw2) * (pfact*(1+Ygetval(codop)));}}}
}/*M2_STOTB   */ 
void YstateTOtab_v(YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_euler; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler;  ++Yw2){
			 for (int Yws=0; Yws<YNBS_v; ++Yws){
				 tab[Y3windice] = YS_v(Yws, Yw1, Yw2); ++Y3windice;}}}
}/*M2_GTOTB   */ 
void YgradTOtab_v(YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_euler; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler;  ++Yw2){
			 for (int Yws=0; Yws<YNBS_v; ++Yws){
				 tab[Y3windice] = YG_v(Yws, Yw1, Yw2); ++Y3windice;}}}
}/*M2_RGRAD*/ 
void Yrazgrad_v()
{
	 for (int Yw1=0; Yw1<YA1_S_euler; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler; ++Yw2){
			 for (int Yws=0; Yws<YNBS_v; ++Yws){
				 Yv(Yw1, Yw2)->Ygrad[Yws] = 0.0;}}}
}/*M2_SETA*/ 
void Ysetstate_v(double val)
{
	 for (int Yw1=0; Yw1<YA1_S_euler; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler;  ++Yw2){
			 for (int Yws=0; Yws<YNBS_v; ++Yws){
				 Yv(Yw1, Yw2)->Ystate[Yws] = val;}}}
}/*M2_YIOTU*/ 
void Yio_v(int Yws, int Yw1, int Yw2, int Yw3, int Ywt, double val)	
{  Yiorwf ("v", Yws, Yw1, Yw2, Yw3, Ywt, &YS_v(Yws, Yw1, Yw2), val);
}/*M2_CREY*/ 
void Ycreate_v()  
{
	 for (int Yw1=0; Yw1<YA1_S_euler; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler; ++Yw2){
			 if ((Yv(Yw1, Yw2)= new Yaov ())==NULL) Ycreate_errnew();}}
}/*M3_TBTOGT  */ 
void YtabTOgrad_locate(int frompdt, int topdt, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_eulerlocate; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_eulerlocate; ++Yw2){
			 for (int Yw3=0; Yw3<YA3_S_eulerlocate; ++Yw3){
				 for (int Ywt=frompdt; Ywt<topdt; ++Ywt){
					 for (int Yws=0; Yws<YNBS_locate; ++Yws){
						 YG_locate(Yws, Yw1, Yw2, Yw3, Ywt) = tab[Y3windice]; ++Y3windice;}}}}}
}/*M3_GEQPDT*/ 
void YgradEQPdelta_locate(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_eulerlocate; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_eulerlocate; ++Yw2){
			 for (int Yw3=0; Yw3<YA3_S_eulerlocate; ++Yw3){
				 for (int Ywt=frompdt; Ywt<topdt; ++Ywt){
					 for (int Yws=0; Yws<YNBS_locate; ++Yws){
						 YG_locate(Yws, Yw1, Yw2, Yw3, Ywt) = YD_locate(Yws, Yw1, Yw2, Yw3, Ywt) * pfact;}}}}}
}/*M3_DEQPST*/ 
void YdeltaEQPstate_locate(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_eulerlocate; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_eulerlocate; ++Yw2){
			 for (int Yw3=0; Yw3<YA3_S_eulerlocate; ++Yw3){
				 for (int Ywt=frompdt; Ywt<topdt; ++Ywt){
					 for (int Yws=0; Yws<YNBS_locate; ++Yws){
						 YD_locate(Yws, Yw1, Yw2, Yw3, Ywt) = YS_locate(Yws, Yw1, Yw2, Yw3, Ywt) * pfact;}}}}}
}/*M3_DEQPGT*/ 
void YdeltaEQPgrad_locate(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_eulerlocate; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_eulerlocate; ++Yw2){
			 for (int Yw3=0; Yw3<YA3_S_eulerlocate; ++Yw3){
				 for (int Ywt=frompdt; Ywt<topdt; ++Ywt){
					 for (int Yws=0; Yws<YNBS_locate; ++Yws){
						 YD_locate(Yws, Yw1, Yw2, Yw3, Ywt) = YG_locate(Yws, Yw1, Yw2, Yw3, Ywt) * pfact;}}}}}
}/*M3_SEQPDT*/ 
void YstateEQPdelta_locate(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_eulerlocate; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_eulerlocate; ++Yw2){
			 for (int Yw3=0; Yw3<YA3_S_eulerlocate; ++Yw3){
				 for (int Ywt=frompdt; Ywt<topdt; ++Ywt){
					 for (int Yws=0; Yws<YNBS_locate; ++Yws){
						 YS_locate(Yws, Yw1, Yw2, Yw3, Ywt) = YD_locate(Yws, Yw1, Yw2, Yw3, Ywt) * pfact;}}}}}
}/*M3_SEQPODT*/ 
void YstateEQPOdelta_locate(int frompdt, int topdt, double pfact, char *codop)
{
	 for (int Yw1=0; Yw1<YA1_S_eulerlocate; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_eulerlocate; ++Yw2){
			 for (int Yw3=0; Yw3<YA3_S_eulerlocate; ++Yw3){
				 for (int Ywt=frompdt; Ywt<topdt; ++Ywt){
					 for (int Yws=0; Yws<YNBS_locate; ++Yws){
						 YS_locate(Yws, Yw1, Yw2, Yw3, Ywt) = YD_locate(Yws, Yw1, Yw2, Yw3, Ywt) * (pfact*(1+Ygetval(codop)));}}}}}
}/*M3_SEQAPTDT*/ 
void YstateEQAPTdelta_locate(int frompdt, int topdt, double pfact, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_eulerlocate; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_eulerlocate; ++Yw2){
			 for (int Yw3=0; Yw3<YA3_S_eulerlocate; ++Yw3){
				 for (int Ywt=frompdt; Ywt<topdt; ++Ywt){
					 for (int Yws=0; Yws<YNBS_locate; ++Yws){
						 YS_locate(Yws, Yw1, Yw2, Yw3, Ywt) = YD_locate(Yws, Yw1, Yw2, Yw3, Ywt) + pfact * tab[Y3windice++];}}}}}
}/*M3_GEQPST*/ 
void YgradEQPstate_locate(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_eulerlocate; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_eulerlocate; ++Yw2){
			 for (int Yw3=0; Yw3<YA3_S_eulerlocate; ++Yw3){
				 for (int Ywt=frompdt; Ywt<topdt; ++Ywt){
					 for (int Yws=0; Yws<YNBS_locate; ++Yws){
						 YG_locate(Yws, Yw1, Yw2, Yw3, Ywt) = YS_locate(Yws, Yw1, Yw2, Yw3, Ywt) * pfact;}}}}}
}/*M3_GEQPOST*/ 
void YgradEQPOstate_locate(int frompdt, int topdt, double pfact, char *codop)
{
	 for (int Yw1=0; Yw1<YA1_S_eulerlocate; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_eulerlocate; ++Yw2){
			 for (int Yw3=0; Yw3<YA3_S_eulerlocate; ++Yw3){
				 for (int Ywt=frompdt; Ywt<topdt; ++Ywt){
					 for (int Yws=0; Yws<YNBS_locate; ++Yws){
						 YG_locate(Yws, Yw1, Yw2, Yw3, Ywt) = YS_locate(Yws, Yw1, Yw2, Yw3, Ywt) * (pfact*(1+Ygetval(codop)));}}}}}
}/*M3_STOTBT  */ 
void YstateTOtab_locate(int frompdt, int topdt, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_eulerlocate; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_eulerlocate; ++Yw2){
			 for (int Yw3=0; Yw3<YA3_S_eulerlocate; ++Yw3){
				 for (int Ywt=frompdt; Ywt<topdt; ++Ywt){
					 for (int Yws=0; Yws<YNBS_locate; ++Yws){
						 tab[Y3windice] = YS_locate(Yws, Yw1, Yw2, Yw3, Ywt); ++Y3windice;}}}}}
}/*M3_GTOTBT  */ 
void YgradTOtab_locate(int frompdt, int topdt, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_eulerlocate; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_eulerlocate; ++Yw2){
			 for (int Yw3=0; Yw3<YA3_S_eulerlocate; ++Yw3){
				 for (int Ywt=frompdt; Ywt<topdt; ++Ywt){
					 for (int Yws=0; Yws<YNBS_locate; ++Yws){
						 tab[Y3windice] = YG_locate(Yws, Yw1, Yw2, Yw3, Ywt); ++Y3windice;}}}}}
}/*M3_RGRADT*/ 
void Yrazgrad_locate()
{
	 for (int Yw1=0; Yw1<YA1_S_eulerlocate; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_eulerlocate; ++Yw2){
			 for (int Yw3=0; Yw3<YA3_S_eulerlocate; ++Yw3){
				 for (int Ywt=0; Ywt<YNBALLTIME_T_float; ++Ywt){
					 for (int Yws=0; Yws<YNBS_locate; ++Yws){
						 Ylocate(Yw1, Yw2, Yw3)->Ygrad[Ywt][Yws] = 0.0;}}}}}
}/*M3_SETAT*/ 
void Ysetstate_locate(double val)
{
	 for (int Yw1=0; Yw1<YA1_S_eulerlocate; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_eulerlocate; ++Yw2){
			 for (int Yw3=0; Yw3<YA3_S_eulerlocate; ++Yw3){
				 for (int Ywt=0; Ywt<YNBALLTIME_T_float; ++Ywt){
					 for (int Yws=0; Yws<YNBS_locate; ++Yws){
						 Ylocate(Yw1, Yw2, Yw3)->Ystate[Ywt][Yws] = val;}}}}}
}/*M3_YIOUT*/ 
void Yio_locate(int Yws, int Yw1, int Yw2, int Yw3, int Ywt, double val)	
{  Yiorwf ("locate", Yws, Yw1, Yw2, Yw3, Ywt, &YS_locate(Yws, Yw1, Yw2, Yw3, Ywt), val);
}/*M00_GTEST */
int Ytestdf_locate (int modop, int KeKo, float pdx, float ptol, int ynbi, int ynbs, char *nmmod, YREAL yinput[], YREAL ystate[], locate *Modul)
{	 int wk, wj, nbKo=0; YREAL wSj[YMAX_NBS],wdxk,wJjk,wyj,wDjk,wdifjk;
	 YREAL svx; int resKo; char okko[2][3]={"ok", "KO"}; 
	 if (!Ydftestijkt(Yid_locate)) return(0); 
	 memset (Yjac, 0, YMAX_JAC_NBI*YMAX_JAC_NBS*sizeof(YREAL));
/*M01_GTESTT*/
	 Modul->forward( Yting[0], Yting[1] );
	 for (wj=0; wj<ynbs; ++wj) wSj[wj] = ystate[YTemps*ynbs + wj];
	 Modul->backward( Yting[0], Yting[1] );
	 for (wk=0; wk<ynbi; ++wk)
	 {	 svx=yinput[wk];
			 wdxk = svx*pdx;
			 if (modop||svx==0.0) wdxk = pdx; 
			 yinput[wk] += wdxk;
			 Modul->forward( Yting[0], Yting[1] );
			 for (wj=0; wj<ynbs; ++wj)
			 { wJjk = Yjac[wj][wk];
				 wyj = ystate[YTemps*ynbs + wj];
				 wDjk = (wyj - wSj[wj]) / wdxk;
				 wdifjk = fabs(wDjk - wJjk);
				 resKo=(wdifjk/fabs(wJjk))>ptol;  nbKo+=resKo;
				 if(resKo || !KeKo) 
				 printf(" %7s | J[%2i,%2i] |% 12.6e |% 12.6e |% 12.6e :  %s |% 12.6e | (%i %i %i %i)\n",
								 nmmod, wj+1, wk+1, wJjk, wDjk, wdifjk, okko[resKo], svx, Yi+1, Yj+1, Yk+1, Yt+1);
				} 
			  yinput[wk] = svx;
	 }
	 return(nbKo);
}/*M3_CREU*/ 
void Ycreate_locate()  
{
	 for (int Yw1=0; Yw1<YA1_S_eulerlocate; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_eulerlocate; ++Yw2){
			 for (int Yw3=0; Yw3<YA3_S_eulerlocate; ++Yw3){
				 if ((Ylocate(Yw1, Yw2, Yw3) = new locate ())==NULL) Ycreate_errnew();}}}
}/*M1_TBTOGT  */ 
void YtabTOgrad_r_float_d(int frompdt, int topdt, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_r_float_d; ++Yws){
				 YG_r_float_d(Yws, Yw1, Ywt) = tab[Y3windice]; ++Y3windice;}}}
}/*M1_GEQPDT*/ 
void YgradEQPdelta_r_float_d(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_r_float_d; ++Yws){
				 YG_r_float_d(Yws, Yw1, Ywt) = YD_r_float_d(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_DEQPST*/ 
void YdeltaEQPstate_r_float_d(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_r_float_d; ++Yws){
				 YD_r_float_d(Yws, Yw1, Ywt) = YS_r_float_d(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_DEQPGT*/ 
void YdeltaEQPgrad_r_float_d(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_r_float_d; ++Yws){
				 YD_r_float_d(Yws, Yw1, Ywt) = YG_r_float_d(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_SEQPDT*/ 
void YstateEQPdelta_r_float_d(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_r_float_d; ++Yws){
				 YS_r_float_d(Yws, Yw1, Ywt) = YD_r_float_d(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_SEQPODT*/ 
void YstateEQPOdelta_r_float_d(int frompdt, int topdt, double pfact, char *codop)
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_r_float_d; ++Yws){
				 YS_r_float_d(Yws, Yw1, Ywt) = YD_r_float_d(Yws, Yw1, Ywt) * (pfact*(1+Ygetval(codop)));}}}
}/*M1_SEQAPTDT*/ 
void YstateEQAPTdelta_r_float_d(int frompdt, int topdt, double pfact, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_r_float_d; ++Yws){
				 YS_r_float_d(Yws, Yw1, Ywt) = YD_r_float_d(Yws, Yw1, Ywt) + pfact * tab[Y3windice++];}}}
}/*M1_GEQPST*/ 
void YgradEQPstate_r_float_d(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_r_float_d; ++Yws){
				 YG_r_float_d(Yws, Yw1, Ywt) = YS_r_float_d(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_GEQPOST*/ 
void YgradEQPOstate_r_float_d(int frompdt, int topdt, double pfact, char *codop)
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_r_float_d; ++Yws){
				 YG_r_float_d(Yws, Yw1, Ywt) = YS_r_float_d(Yws, Yw1, Ywt) * (pfact*(1+Ygetval(codop)));}}}
}/*M1_STOTBT  */ 
void YstateTOtab_r_float_d(int frompdt, int topdt, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_r_float_d; ++Yws){
				 tab[Y3windice] = YS_r_float_d(Yws, Yw1, Ywt); ++Y3windice;}}}
}/*M1_GTOTBT  */ 
void YgradTOtab_r_float_d(int frompdt, int topdt, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_r_float_d; ++Yws){
				 tab[Y3windice] = YG_r_float_d(Yws, Yw1, Ywt); ++Y3windice;}}}
}/*M1_RGRADT*/ 
void Yrazgrad_r_float_d()
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=0; Ywt<YNBALLTIME_T_float_d; ++Ywt){
			 for (int Yws=0; Yws<YNBS_r_float_d; ++Yws){
				 Yr_float_d(Yw1)->Ygrad[Ywt][Yws] = 0.0;}}}
}/*M1_SETAT*/ 
void Ysetstate_r_float_d(double val)
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=0; Ywt<YNBALLTIME_T_float_d;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_r_float_d; ++Yws){
				 Yr_float_d(Yw1)->Ystate[Ywt][Yws] = val;}}}
}/*M1_YIOUT*/ 
void Yio_r_float_d(int Yws, int Yw1, int Yw2, int Yw3, int Ywt, double val)	
{  Yiorwf ("r_float_d", Yws, Yw1, Yw2, Yw3, Ywt, &YS_r_float_d(Yws, Yw1, Ywt), val);
}/*M00_GTEST */
int Ytestdf_r_float_d (int modop, int KeKo, float pdx, float ptol, int ynbi, int ynbs, char *nmmod, YREAL yinput[], YREAL ystate[], r_float_d *Modul)
{	 int wk, wj, nbKo=0; YREAL wSj[YMAX_NBS],wdxk,wJjk,wyj,wDjk,wdifjk;
	 YREAL svx; int resKo; char okko[2][3]={"ok", "KO"}; 
	 if (!Ydftestijkt(Yid_r_float_d)) return(0); 
	 memset (Yjac, 0, YMAX_JAC_NBI*YMAX_JAC_NBS*sizeof(YREAL));
/*M01_GTESTT*/
	 Modul->forward( Yting[0], Yting[1], Yting[2], Yting[3] );
	 for (wj=0; wj<ynbs; ++wj) wSj[wj] = ystate[YTemps*ynbs + wj];
	 Modul->backward( Yting[0], Yting[1], Yting[2], Yting[3] );
	 for (wk=0; wk<ynbi; ++wk)
	 {	 svx=yinput[wk];
			 wdxk = svx*pdx;
			 if (modop||svx==0.0) wdxk = pdx; 
			 yinput[wk] += wdxk;
			 Modul->forward( Yting[0], Yting[1], Yting[2], Yting[3] );
			 for (wj=0; wj<ynbs; ++wj)
			 { wJjk = Yjac[wj][wk];
				 wyj = ystate[YTemps*ynbs + wj];
				 wDjk = (wyj - wSj[wj]) / wdxk;
				 wdifjk = fabs(wDjk - wJjk);
				 resKo=(wdifjk/fabs(wJjk))>ptol;  nbKo+=resKo;
				 if(resKo || !KeKo) 
				 printf(" %7s | J[%2i,%2i] |% 12.6e |% 12.6e |% 12.6e :  %s |% 12.6e | (%i %i %i %i)\n",
								 nmmod, wj+1, wk+1, wJjk, wDjk, wdifjk, okko[resKo], svx, Yi+1, Yj+1, Yk+1, Yt+1);
				} 
			  yinput[wk] = svx;
	 }
	 return(nbKo);
}/*M1_CREU*/ 
void Ycreate_r_float_d()  
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 if ((Yr_float_d(Yw1) = new r_float_d ())==NULL) Ycreate_errnew();}
}/*M1_TBTOGT  */ 
void YtabTOgrad_ur_d(int frompdt, int topdt, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_ur_d; ++Yws){
				 YG_ur_d(Yws, Yw1, Ywt) = tab[Y3windice]; ++Y3windice;}}}
}/*M1_GEQPDT*/ 
void YgradEQPdelta_ur_d(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_ur_d; ++Yws){
				 YG_ur_d(Yws, Yw1, Ywt) = YD_ur_d(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_DEQPST*/ 
void YdeltaEQPstate_ur_d(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_ur_d; ++Yws){
				 YD_ur_d(Yws, Yw1, Ywt) = YS_ur_d(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_DEQPGT*/ 
void YdeltaEQPgrad_ur_d(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_ur_d; ++Yws){
				 YD_ur_d(Yws, Yw1, Ywt) = YG_ur_d(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_SEQPDT*/ 
void YstateEQPdelta_ur_d(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_ur_d; ++Yws){
				 YS_ur_d(Yws, Yw1, Ywt) = YD_ur_d(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_SEQPODT*/ 
void YstateEQPOdelta_ur_d(int frompdt, int topdt, double pfact, char *codop)
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_ur_d; ++Yws){
				 YS_ur_d(Yws, Yw1, Ywt) = YD_ur_d(Yws, Yw1, Ywt) * (pfact*(1+Ygetval(codop)));}}}
}/*M1_SEQAPTDT*/ 
void YstateEQAPTdelta_ur_d(int frompdt, int topdt, double pfact, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_ur_d; ++Yws){
				 YS_ur_d(Yws, Yw1, Ywt) = YD_ur_d(Yws, Yw1, Ywt) + pfact * tab[Y3windice++];}}}
}/*M1_GEQPST*/ 
void YgradEQPstate_ur_d(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_ur_d; ++Yws){
				 YG_ur_d(Yws, Yw1, Ywt) = YS_ur_d(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_GEQPOST*/ 
void YgradEQPOstate_ur_d(int frompdt, int topdt, double pfact, char *codop)
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_ur_d; ++Yws){
				 YG_ur_d(Yws, Yw1, Ywt) = YS_ur_d(Yws, Yw1, Ywt) * (pfact*(1+Ygetval(codop)));}}}
}/*M1_STOTBT  */ 
void YstateTOtab_ur_d(int frompdt, int topdt, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_ur_d; ++Yws){
				 tab[Y3windice] = YS_ur_d(Yws, Yw1, Ywt); ++Y3windice;}}}
}/*M1_GTOTBT  */ 
void YgradTOtab_ur_d(int frompdt, int topdt, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_ur_d; ++Yws){
				 tab[Y3windice] = YG_ur_d(Yws, Yw1, Ywt); ++Y3windice;}}}
}/*M1_RGRADT*/ 
void Yrazgrad_ur_d()
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=0; Ywt<YNBALLTIME_T_float_d; ++Ywt){
			 for (int Yws=0; Yws<YNBS_ur_d; ++Yws){
				 Yur_d(Yw1)->Ygrad[Ywt][Yws] = 0.0;}}}
}/*M1_SETAT*/ 
void Ysetstate_ur_d(double val)
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=0; Ywt<YNBALLTIME_T_float_d;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_ur_d; ++Yws){
				 Yur_d(Yw1)->Ystate[Ywt][Yws] = val;}}}
}/*M1_YIOUT*/ 
void Yio_ur_d(int Yws, int Yw1, int Yw2, int Yw3, int Ywt, double val)	
{  Yiorwf ("ur_d", Yws, Yw1, Yw2, Yw3, Ywt, &YS_ur_d(Yws, Yw1, Ywt), val);
}/*M00_GTEST */
int Ytestdf_ur_d (int modop, int KeKo, float pdx, float ptol, int ynbi, int ynbs, char *nmmod, YREAL yinput[], YREAL ystate[], ur_d *Modul)
{	 int wk, wj, nbKo=0; YREAL wSj[YMAX_NBS],wdxk,wJjk,wyj,wDjk,wdifjk;
	 YREAL svx; int resKo; char okko[2][3]={"ok", "KO"}; 
	 if (!Ydftestijkt(Yid_ur_d)) return(0); 
	 memset (Yjac, 0, YMAX_JAC_NBI*YMAX_JAC_NBS*sizeof(YREAL));
/*M01_GTESTT*/
	 Modul->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4], Yting[5], Yting[6], Yting[7], Yting[8], Yting[9], Yting[10], Yting[11], Yting[12], Yting[13], Yting[14], Yting[15], Yting[16], Yting[17] );
	 for (wj=0; wj<ynbs; ++wj) wSj[wj] = ystate[YTemps*ynbs + wj];
	 Modul->backward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4], Yting[5], Yting[6], Yting[7], Yting[8], Yting[9], Yting[10], Yting[11], Yting[12], Yting[13], Yting[14], Yting[15], Yting[16], Yting[17] );
	 for (wk=0; wk<ynbi; ++wk)
	 {	 svx=yinput[wk];
			 wdxk = svx*pdx;
			 if (modop||svx==0.0) wdxk = pdx; 
			 yinput[wk] += wdxk;
			 Modul->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4], Yting[5], Yting[6], Yting[7], Yting[8], Yting[9], Yting[10], Yting[11], Yting[12], Yting[13], Yting[14], Yting[15], Yting[16], Yting[17] );
			 for (wj=0; wj<ynbs; ++wj)
			 { wJjk = Yjac[wj][wk];
				 wyj = ystate[YTemps*ynbs + wj];
				 wDjk = (wyj - wSj[wj]) / wdxk;
				 wdifjk = fabs(wDjk - wJjk);
				 resKo=(wdifjk/fabs(wJjk))>ptol;  nbKo+=resKo;
				 if(resKo || !KeKo) 
				 printf(" %7s | J[%2i,%2i] |% 12.6e |% 12.6e |% 12.6e :  %s |% 12.6e | (%i %i %i %i)\n",
								 nmmod, wj+1, wk+1, wJjk, wDjk, wdifjk, okko[resKo], svx, Yi+1, Yj+1, Yk+1, Yt+1);
				} 
			  yinput[wk] = svx;
	 }
	 return(nbKo);
}/*M1_CREU*/ 
void Ycreate_ur_d()  
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 if ((Yur_d(Yw1) = new ur_d ())==NULL) Ycreate_errnew();}
}/*M1_TBTOGT  */ 
void YtabTOgrad_uinter_d(int frompdt, int topdt, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_uinter_d; ++Yws){
				 YG_uinter_d(Yws, Yw1, Ywt) = tab[Y3windice]; ++Y3windice;}}}
}/*M1_GEQPDT*/ 
void YgradEQPdelta_uinter_d(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_uinter_d; ++Yws){
				 YG_uinter_d(Yws, Yw1, Ywt) = YD_uinter_d(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_DEQPST*/ 
void YdeltaEQPstate_uinter_d(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_uinter_d; ++Yws){
				 YD_uinter_d(Yws, Yw1, Ywt) = YS_uinter_d(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_DEQPGT*/ 
void YdeltaEQPgrad_uinter_d(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_uinter_d; ++Yws){
				 YD_uinter_d(Yws, Yw1, Ywt) = YG_uinter_d(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_SEQPDT*/ 
void YstateEQPdelta_uinter_d(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_uinter_d; ++Yws){
				 YS_uinter_d(Yws, Yw1, Ywt) = YD_uinter_d(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_SEQPODT*/ 
void YstateEQPOdelta_uinter_d(int frompdt, int topdt, double pfact, char *codop)
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_uinter_d; ++Yws){
				 YS_uinter_d(Yws, Yw1, Ywt) = YD_uinter_d(Yws, Yw1, Ywt) * (pfact*(1+Ygetval(codop)));}}}
}/*M1_SEQAPTDT*/ 
void YstateEQAPTdelta_uinter_d(int frompdt, int topdt, double pfact, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_uinter_d; ++Yws){
				 YS_uinter_d(Yws, Yw1, Ywt) = YD_uinter_d(Yws, Yw1, Ywt) + pfact * tab[Y3windice++];}}}
}/*M1_GEQPST*/ 
void YgradEQPstate_uinter_d(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_uinter_d; ++Yws){
				 YG_uinter_d(Yws, Yw1, Ywt) = YS_uinter_d(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_GEQPOST*/ 
void YgradEQPOstate_uinter_d(int frompdt, int topdt, double pfact, char *codop)
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_uinter_d; ++Yws){
				 YG_uinter_d(Yws, Yw1, Ywt) = YS_uinter_d(Yws, Yw1, Ywt) * (pfact*(1+Ygetval(codop)));}}}
}/*M1_STOTBT  */ 
void YstateTOtab_uinter_d(int frompdt, int topdt, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_uinter_d; ++Yws){
				 tab[Y3windice] = YS_uinter_d(Yws, Yw1, Ywt); ++Y3windice;}}}
}/*M1_GTOTBT  */ 
void YgradTOtab_uinter_d(int frompdt, int topdt, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_uinter_d; ++Yws){
				 tab[Y3windice] = YG_uinter_d(Yws, Yw1, Ywt); ++Y3windice;}}}
}/*M1_RGRADT*/ 
void Yrazgrad_uinter_d()
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=0; Ywt<YNBALLTIME_T_float_d; ++Ywt){
			 for (int Yws=0; Yws<YNBS_uinter_d; ++Yws){
				 Yuinter_d(Yw1)->Ygrad[Ywt][Yws] = 0.0;}}}
}/*M1_SETAT*/ 
void Ysetstate_uinter_d(double val)
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=0; Ywt<YNBALLTIME_T_float_d;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_uinter_d; ++Yws){
				 Yuinter_d(Yw1)->Ystate[Ywt][Yws] = val;}}}
}/*M1_YIOUT*/ 
void Yio_uinter_d(int Yws, int Yw1, int Yw2, int Yw3, int Ywt, double val)	
{  Yiorwf ("uinter_d", Yws, Yw1, Yw2, Yw3, Ywt, &YS_uinter_d(Yws, Yw1, Ywt), val);
}/*M1_CREY*/ 
void Ycreate_uinter_d() 	
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 if ((Yuinter_d(Yw1) = new Yaouinter_d ())==NULL) Ycreate_errnew();}
}/*M2_ADJUS*/ 
void Yadjust_u_d()
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1) 
	 for (int Yw2=0; Yw2<YA2_S_euler_d; ++Yw2) 
	 for (int Yws=0; Yws<YNBS_u_d; ++Yws) 
	 { Ystate_u_d(Yws, Yw1, Yw2) -= Yepsi_u_d(Yws, Yw1, Yw2) * Ygrad_u_d(Yws, Yw1, Yw2);
	}
}/*M2_GCTOTB  */ 
void YgradCTOtab_u_d(YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d;  ++Yw2){
			 for (int Yws=0; Yws<YNBS_u_d; ++Yws){
				 tab[Y3windice] += YG_u_d(Yws, Yw1, Yw2); ++Y3windice;}}}
}/*M2_VSTA*/ 
void Y3valstate_u_d()
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d;  ++Yw2){
			 for (int Yws=0; Yws<YNBS_u_d; ++Yws){
				 Y3x[Y3windice] = Yu_d(Yw1, Yw2)->Ystate[Yws]; ++Y3windice;}}}
}/*M2_GSTA*/ 
void Y3getstate_u_d(float x[])
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d;  ++Yw2){
			 for (int Yws=0; Yws<YNBS_u_d; ++Yws){
				 Yu_d(Yw1, Yw2)->Ystate[Yws] = x[Y3windice]; ++Y3windice;}}}
}/*M2_VGRAD*/ 
void Y3valgrad_u_d(float g[])
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d;  ++Yw2){
			 for (int Yws=0; Yws<YNBS_u_d; ++Yws){
				 g[Y3windice] = Yu_d(Yw1, Yw2)->Ygrad[Yws]; ++Y3windice;}}}
}/*M2_TBTOG   */ 
void YtabTOgrad_u_d(YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d;  ++Yw2){
			 for (int Yws=0; Yws<YNBS_u_d; ++Yws){
				 YG_u_d(Yws, Yw1, Yw2) = tab[Y3windice]; ++Y3windice;}}}
}/*M2_GEQPD */ 
void YgradEQPdelta_u_d(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d;  ++Yw2){
			 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_u_d; ++Yws){
				 YG_u_d(Yws, Yw1, Yw2) = YD_u_d(Yws, Yw1, Yw2, Ywt) * pfact;}}}}
}/*M2_DEQPS */ 
void YdeltaEQPstate_u_d(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d;  ++Yw2){
			 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_u_d; ++Yws){
				 YD_u_d(Yws, Yw1, Yw2, Ywt) = YS_u_d(Yws, Yw1, Yw2) * pfact;}}}}
}/*M2_DEQPG */ 
void YdeltaEQPgrad_u_d(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d;  ++Yw2){
			 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_u_d; ++Yws){
				 YD_u_d(Yws, Yw1, Yw2, Ywt) = YG_u_d(Yws, Yw1, Yw2) * pfact;}}}}
}/*M2_SEQPD */ 
void YstateEQPdelta_u_d(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d;  ++Yw2){
			 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_u_d; ++Yws){
				 YS_u_d(Yws, Yw1, Yw2) = YD_u_d(Yws, Yw1, Yw2, Ywt) * pfact;}}}}
}/*M2_SEQPOD */ 
void YstateEQPOdelta_u_d(int frompdt, int topdt, double pfact, char *codop)
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d;  ++Yw2){
			 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_u_d; ++Yws){
				 YS_u_d(Yws, Yw1, Yw2) = YD_u_d(Yws, Yw1, Yw2, Ywt) * (pfact*(1+Ygetval(codop)));}}}}
}/*M2_SEQAPTD */ 
void YstateEQAPTdelta_u_d(int frompdt, int topdt, double pfact, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d;  ++Yw2){
			 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_u_d; ++Yws){
				 YS_u_d(Yws, Yw1, Yw2) = YD_u_d(Yws, Yw1, Yw2, Ywt) + pfact * tab[Y3windice++];}}}}
}/*M2_GEQPS */ 
void YgradEQPstate_u_d(double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d;  ++Yw2){
			 for (int Yws=0; Yws<YNBS_u_d; ++Yws){
				 YG_u_d(Yws, Yw1, Yw2) = YS_u_d(Yws, Yw1, Yw2) * pfact;}}}
}/*M2_GEQPOS */ 
void YgradEQPOstate_u_d(double pfact, char *codop)
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d;  ++Yw2){
			 for (int Yws=0; Yws<YNBS_u_d; ++Yws){
				 YG_u_d(Yws, Yw1, Yw2) = YS_u_d(Yws, Yw1, Yw2) * (pfact*(1+Ygetval(codop)));}}}
}/*M2_STOTB   */ 
void YstateTOtab_u_d(YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d;  ++Yw2){
			 for (int Yws=0; Yws<YNBS_u_d; ++Yws){
				 tab[Y3windice] = YS_u_d(Yws, Yw1, Yw2); ++Y3windice;}}}
}/*M2_GTOTB   */ 
void YgradTOtab_u_d(YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d;  ++Yw2){
			 for (int Yws=0; Yws<YNBS_u_d; ++Yws){
				 tab[Y3windice] = YG_u_d(Yws, Yw1, Yw2); ++Y3windice;}}}
}/*M2_RGRAD*/ 
void Yrazgrad_u_d()
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d; ++Yw2){
			 for (int Yws=0; Yws<YNBS_u_d; ++Yws){
				 Yu_d(Yw1, Yw2)->Ygrad[Yws] = 0.0;}}}
}/*M2_SETA*/ 
void Ysetstate_u_d(double val)
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d;  ++Yw2){
			 for (int Yws=0; Yws<YNBS_u_d; ++Yws){
				 Yu_d(Yw1, Yw2)->Ystate[Yws] = val;}}}
}/*M2_YIOTU*/ 
void Yio_u_d(int Yws, int Yw1, int Yw2, int Yw3, int Ywt, double val)	
{  Yiorwf ("u_d", Yws, Yw1, Yw2, Yw3, Ywt, &YS_u_d(Yws, Yw1, Yw2), val);
}/*M0_OUTOB*/ 
int Youtoobs_u_d(YioKind yiokind, int Yws, int Ywt, int arbpdt)
{
	 Yaobs.time = arbpdt;
	 Yaobs.imod = Yimod("u_d");
	 Yaobs.smod = Yws;/*M2_OUTOB*/ 
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1)
	 for (int Yw2=0; Yw2<YA2_S_euler_d; ++Yw2)
	 {	if(select_io(yiokind, "u_d", Yws, Yw1, Yw2, -1, Ywt, &Ystate_u_d(Yws, Yw1, Yw2)))
			{	Yaobs.iaxe = Yw1; Yaobs.jaxe = Yw2; Yaobs.kaxe = -1;
				Yaobs.vobs = Ystate_u_d(Yws, Yw1, Yw2);
				if(!Yobs_insert (&Yaobs)) return(0);
			}
	 }/*M9_OUTOB*/  
	 return(1);
}/*M2_SEPS*/ 
void Ysetepsi_u_d(double val)
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d;  ++Yw2){
			 for (int Yws=0; Yws<YNBS_u_d; ++Yws){
				 Yu_d(Yw1, Yw2)->Yepsi[Yws] = val;}}}
}/*M2_CREY*/ 
void Ycreate_u_d()  
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d; ++Yw2){
			 if ((Yu_d(Yw1, Yw2)= new Yaou_d ())==NULL) Ycreate_errnew();}}
}/*M2_SWISH*/ 
void Ysetwish_u_d(double val)
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d;  ++Yw2){
			 for (int Yws=0; Yws<YNBS_u_d; ++Yws){
				 Ywish_u_d(Yws, Yw1, Yw2) = val;}}}
}/*M2_ADJUS*/ 
void Yadjust_v_d()
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1) 
	 for (int Yw2=0; Yw2<YA2_S_euler_d; ++Yw2) 
	 for (int Yws=0; Yws<YNBS_v_d; ++Yws) 
	 { Ystate_v_d(Yws, Yw1, Yw2) -= Yepsi_v_d(Yws, Yw1, Yw2) * Ygrad_v_d(Yws, Yw1, Yw2);
	}
}/*M2_GCTOTB  */ 
void YgradCTOtab_v_d(YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d;  ++Yw2){
			 for (int Yws=0; Yws<YNBS_v_d; ++Yws){
				 tab[Y3windice] += YG_v_d(Yws, Yw1, Yw2); ++Y3windice;}}}
}/*M2_VSTA*/ 
void Y3valstate_v_d()
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d;  ++Yw2){
			 for (int Yws=0; Yws<YNBS_v_d; ++Yws){
				 Y3x[Y3windice] = Yv_d(Yw1, Yw2)->Ystate[Yws]; ++Y3windice;}}}
}/*M2_GSTA*/ 
void Y3getstate_v_d(float x[])
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d;  ++Yw2){
			 for (int Yws=0; Yws<YNBS_v_d; ++Yws){
				 Yv_d(Yw1, Yw2)->Ystate[Yws] = x[Y3windice]; ++Y3windice;}}}
}/*M2_VGRAD*/ 
void Y3valgrad_v_d(float g[])
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d;  ++Yw2){
			 for (int Yws=0; Yws<YNBS_v_d; ++Yws){
				 g[Y3windice] = Yv_d(Yw1, Yw2)->Ygrad[Yws]; ++Y3windice;}}}
}/*M2_TBTOG   */ 
void YtabTOgrad_v_d(YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d;  ++Yw2){
			 for (int Yws=0; Yws<YNBS_v_d; ++Yws){
				 YG_v_d(Yws, Yw1, Yw2) = tab[Y3windice]; ++Y3windice;}}}
}/*M2_GEQPD */ 
void YgradEQPdelta_v_d(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d;  ++Yw2){
			 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_v_d; ++Yws){
				 YG_v_d(Yws, Yw1, Yw2) = YD_v_d(Yws, Yw1, Yw2, Ywt) * pfact;}}}}
}/*M2_DEQPS */ 
void YdeltaEQPstate_v_d(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d;  ++Yw2){
			 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_v_d; ++Yws){
				 YD_v_d(Yws, Yw1, Yw2, Ywt) = YS_v_d(Yws, Yw1, Yw2) * pfact;}}}}
}/*M2_DEQPG */ 
void YdeltaEQPgrad_v_d(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d;  ++Yw2){
			 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_v_d; ++Yws){
				 YD_v_d(Yws, Yw1, Yw2, Ywt) = YG_v_d(Yws, Yw1, Yw2) * pfact;}}}}
}/*M2_SEQPD */ 
void YstateEQPdelta_v_d(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d;  ++Yw2){
			 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_v_d; ++Yws){
				 YS_v_d(Yws, Yw1, Yw2) = YD_v_d(Yws, Yw1, Yw2, Ywt) * pfact;}}}}
}/*M2_SEQPOD */ 
void YstateEQPOdelta_v_d(int frompdt, int topdt, double pfact, char *codop)
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d;  ++Yw2){
			 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_v_d; ++Yws){
				 YS_v_d(Yws, Yw1, Yw2) = YD_v_d(Yws, Yw1, Yw2, Ywt) * (pfact*(1+Ygetval(codop)));}}}}
}/*M2_SEQAPTD */ 
void YstateEQAPTdelta_v_d(int frompdt, int topdt, double pfact, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d;  ++Yw2){
			 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_v_d; ++Yws){
				 YS_v_d(Yws, Yw1, Yw2) = YD_v_d(Yws, Yw1, Yw2, Ywt) + pfact * tab[Y3windice++];}}}}
}/*M2_GEQPS */ 
void YgradEQPstate_v_d(double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d;  ++Yw2){
			 for (int Yws=0; Yws<YNBS_v_d; ++Yws){
				 YG_v_d(Yws, Yw1, Yw2) = YS_v_d(Yws, Yw1, Yw2) * pfact;}}}
}/*M2_GEQPOS */ 
void YgradEQPOstate_v_d(double pfact, char *codop)
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d;  ++Yw2){
			 for (int Yws=0; Yws<YNBS_v_d; ++Yws){
				 YG_v_d(Yws, Yw1, Yw2) = YS_v_d(Yws, Yw1, Yw2) * (pfact*(1+Ygetval(codop)));}}}
}/*M2_STOTB   */ 
void YstateTOtab_v_d(YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d;  ++Yw2){
			 for (int Yws=0; Yws<YNBS_v_d; ++Yws){
				 tab[Y3windice] = YS_v_d(Yws, Yw1, Yw2); ++Y3windice;}}}
}/*M2_GTOTB   */ 
void YgradTOtab_v_d(YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d;  ++Yw2){
			 for (int Yws=0; Yws<YNBS_v_d; ++Yws){
				 tab[Y3windice] = YG_v_d(Yws, Yw1, Yw2); ++Y3windice;}}}
}/*M2_RGRAD*/ 
void Yrazgrad_v_d()
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d; ++Yw2){
			 for (int Yws=0; Yws<YNBS_v_d; ++Yws){
				 Yv_d(Yw1, Yw2)->Ygrad[Yws] = 0.0;}}}
}/*M2_SETA*/ 
void Ysetstate_v_d(double val)
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d;  ++Yw2){
			 for (int Yws=0; Yws<YNBS_v_d; ++Yws){
				 Yv_d(Yw1, Yw2)->Ystate[Yws] = val;}}}
}/*M2_YIOTU*/ 
void Yio_v_d(int Yws, int Yw1, int Yw2, int Yw3, int Ywt, double val)	
{  Yiorwf ("v_d", Yws, Yw1, Yw2, Yw3, Ywt, &YS_v_d(Yws, Yw1, Yw2), val);
}/*M0_OUTOB*/ 
int Youtoobs_v_d(YioKind yiokind, int Yws, int Ywt, int arbpdt)
{
	 Yaobs.time = arbpdt;
	 Yaobs.imod = Yimod("v_d");
	 Yaobs.smod = Yws;/*M2_OUTOB*/ 
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1)
	 for (int Yw2=0; Yw2<YA2_S_euler_d; ++Yw2)
	 {	if(select_io(yiokind, "v_d", Yws, Yw1, Yw2, -1, Ywt, &Ystate_v_d(Yws, Yw1, Yw2)))
			{	Yaobs.iaxe = Yw1; Yaobs.jaxe = Yw2; Yaobs.kaxe = -1;
				Yaobs.vobs = Ystate_v_d(Yws, Yw1, Yw2);
				if(!Yobs_insert (&Yaobs)) return(0);
			}
	 }/*M9_OUTOB*/  
	 return(1);
}/*M2_SEPS*/ 
void Ysetepsi_v_d(double val)
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d;  ++Yw2){
			 for (int Yws=0; Yws<YNBS_v_d; ++Yws){
				 Yv_d(Yw1, Yw2)->Yepsi[Yws] = val;}}}
}/*M2_CREY*/ 
void Ycreate_v_d()  
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d; ++Yw2){
			 if ((Yv_d(Yw1, Yw2)= new Yaov_d ())==NULL) Ycreate_errnew();}}
}/*M2_SWISH*/ 
void Ysetwish_v_d(double val)
{
	 for (int Yw1=0; Yw1<YA1_S_euler_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_euler_d;  ++Yw2){
			 for (int Yws=0; Yws<YNBS_v_d; ++Yws){
				 Ywish_v_d(Yws, Yw1, Yw2) = val;}}}
}/*M3_TBTOGT  */ 
void YtabTOgrad_locate_d(int frompdt, int topdt, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_eulerlocate_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_eulerlocate_d; ++Yw2){
			 for (int Yw3=0; Yw3<YA3_S_eulerlocate_d; ++Yw3){
				 for (int Ywt=frompdt; Ywt<topdt; ++Ywt){
					 for (int Yws=0; Yws<YNBS_locate_d; ++Yws){
						 YG_locate_d(Yws, Yw1, Yw2, Yw3, Ywt) = tab[Y3windice]; ++Y3windice;}}}}}
}/*M3_GEQPDT*/ 
void YgradEQPdelta_locate_d(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_eulerlocate_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_eulerlocate_d; ++Yw2){
			 for (int Yw3=0; Yw3<YA3_S_eulerlocate_d; ++Yw3){
				 for (int Ywt=frompdt; Ywt<topdt; ++Ywt){
					 for (int Yws=0; Yws<YNBS_locate_d; ++Yws){
						 YG_locate_d(Yws, Yw1, Yw2, Yw3, Ywt) = YD_locate_d(Yws, Yw1, Yw2, Yw3, Ywt) * pfact;}}}}}
}/*M3_DEQPST*/ 
void YdeltaEQPstate_locate_d(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_eulerlocate_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_eulerlocate_d; ++Yw2){
			 for (int Yw3=0; Yw3<YA3_S_eulerlocate_d; ++Yw3){
				 for (int Ywt=frompdt; Ywt<topdt; ++Ywt){
					 for (int Yws=0; Yws<YNBS_locate_d; ++Yws){
						 YD_locate_d(Yws, Yw1, Yw2, Yw3, Ywt) = YS_locate_d(Yws, Yw1, Yw2, Yw3, Ywt) * pfact;}}}}}
}/*M3_DEQPGT*/ 
void YdeltaEQPgrad_locate_d(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_eulerlocate_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_eulerlocate_d; ++Yw2){
			 for (int Yw3=0; Yw3<YA3_S_eulerlocate_d; ++Yw3){
				 for (int Ywt=frompdt; Ywt<topdt; ++Ywt){
					 for (int Yws=0; Yws<YNBS_locate_d; ++Yws){
						 YD_locate_d(Yws, Yw1, Yw2, Yw3, Ywt) = YG_locate_d(Yws, Yw1, Yw2, Yw3, Ywt) * pfact;}}}}}
}/*M3_SEQPDT*/ 
void YstateEQPdelta_locate_d(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_eulerlocate_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_eulerlocate_d; ++Yw2){
			 for (int Yw3=0; Yw3<YA3_S_eulerlocate_d; ++Yw3){
				 for (int Ywt=frompdt; Ywt<topdt; ++Ywt){
					 for (int Yws=0; Yws<YNBS_locate_d; ++Yws){
						 YS_locate_d(Yws, Yw1, Yw2, Yw3, Ywt) = YD_locate_d(Yws, Yw1, Yw2, Yw3, Ywt) * pfact;}}}}}
}/*M3_SEQPODT*/ 
void YstateEQPOdelta_locate_d(int frompdt, int topdt, double pfact, char *codop)
{
	 for (int Yw1=0; Yw1<YA1_S_eulerlocate_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_eulerlocate_d; ++Yw2){
			 for (int Yw3=0; Yw3<YA3_S_eulerlocate_d; ++Yw3){
				 for (int Ywt=frompdt; Ywt<topdt; ++Ywt){
					 for (int Yws=0; Yws<YNBS_locate_d; ++Yws){
						 YS_locate_d(Yws, Yw1, Yw2, Yw3, Ywt) = YD_locate_d(Yws, Yw1, Yw2, Yw3, Ywt) * (pfact*(1+Ygetval(codop)));}}}}}
}/*M3_SEQAPTDT*/ 
void YstateEQAPTdelta_locate_d(int frompdt, int topdt, double pfact, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_eulerlocate_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_eulerlocate_d; ++Yw2){
			 for (int Yw3=0; Yw3<YA3_S_eulerlocate_d; ++Yw3){
				 for (int Ywt=frompdt; Ywt<topdt; ++Ywt){
					 for (int Yws=0; Yws<YNBS_locate_d; ++Yws){
						 YS_locate_d(Yws, Yw1, Yw2, Yw3, Ywt) = YD_locate_d(Yws, Yw1, Yw2, Yw3, Ywt) + pfact * tab[Y3windice++];}}}}}
}/*M3_GEQPST*/ 
void YgradEQPstate_locate_d(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_eulerlocate_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_eulerlocate_d; ++Yw2){
			 for (int Yw3=0; Yw3<YA3_S_eulerlocate_d; ++Yw3){
				 for (int Ywt=frompdt; Ywt<topdt; ++Ywt){
					 for (int Yws=0; Yws<YNBS_locate_d; ++Yws){
						 YG_locate_d(Yws, Yw1, Yw2, Yw3, Ywt) = YS_locate_d(Yws, Yw1, Yw2, Yw3, Ywt) * pfact;}}}}}
}/*M3_GEQPOST*/ 
void YgradEQPOstate_locate_d(int frompdt, int topdt, double pfact, char *codop)
{
	 for (int Yw1=0; Yw1<YA1_S_eulerlocate_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_eulerlocate_d; ++Yw2){
			 for (int Yw3=0; Yw3<YA3_S_eulerlocate_d; ++Yw3){
				 for (int Ywt=frompdt; Ywt<topdt; ++Ywt){
					 for (int Yws=0; Yws<YNBS_locate_d; ++Yws){
						 YG_locate_d(Yws, Yw1, Yw2, Yw3, Ywt) = YS_locate_d(Yws, Yw1, Yw2, Yw3, Ywt) * (pfact*(1+Ygetval(codop)));}}}}}
}/*M3_STOTBT  */ 
void YstateTOtab_locate_d(int frompdt, int topdt, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_eulerlocate_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_eulerlocate_d; ++Yw2){
			 for (int Yw3=0; Yw3<YA3_S_eulerlocate_d; ++Yw3){
				 for (int Ywt=frompdt; Ywt<topdt; ++Ywt){
					 for (int Yws=0; Yws<YNBS_locate_d; ++Yws){
						 tab[Y3windice] = YS_locate_d(Yws, Yw1, Yw2, Yw3, Ywt); ++Y3windice;}}}}}
}/*M3_GTOTBT  */ 
void YgradTOtab_locate_d(int frompdt, int topdt, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_eulerlocate_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_eulerlocate_d; ++Yw2){
			 for (int Yw3=0; Yw3<YA3_S_eulerlocate_d; ++Yw3){
				 for (int Ywt=frompdt; Ywt<topdt; ++Ywt){
					 for (int Yws=0; Yws<YNBS_locate_d; ++Yws){
						 tab[Y3windice] = YG_locate_d(Yws, Yw1, Yw2, Yw3, Ywt); ++Y3windice;}}}}}
}/*M3_RGRADT*/ 
void Yrazgrad_locate_d()
{
	 for (int Yw1=0; Yw1<YA1_S_eulerlocate_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_eulerlocate_d; ++Yw2){
			 for (int Yw3=0; Yw3<YA3_S_eulerlocate_d; ++Yw3){
				 for (int Ywt=0; Ywt<YNBALLTIME_T_float_d; ++Ywt){
					 for (int Yws=0; Yws<YNBS_locate_d; ++Yws){
						 Ylocate_d(Yw1, Yw2, Yw3)->Ygrad[Ywt][Yws] = 0.0;}}}}}
}/*M3_SETAT*/ 
void Ysetstate_locate_d(double val)
{
	 for (int Yw1=0; Yw1<YA1_S_eulerlocate_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_eulerlocate_d; ++Yw2){
			 for (int Yw3=0; Yw3<YA3_S_eulerlocate_d; ++Yw3){
				 for (int Ywt=0; Ywt<YNBALLTIME_T_float_d; ++Ywt){
					 for (int Yws=0; Yws<YNBS_locate_d; ++Yws){
						 Ylocate_d(Yw1, Yw2, Yw3)->Ystate[Ywt][Yws] = val;}}}}}
}/*M3_YIOUT*/ 
void Yio_locate_d(int Yws, int Yw1, int Yw2, int Yw3, int Ywt, double val)	
{  Yiorwf ("locate_d", Yws, Yw1, Yw2, Yw3, Ywt, &YS_locate_d(Yws, Yw1, Yw2, Yw3, Ywt), val);
}/*M00_GTEST */
int Ytestdf_locate_d (int modop, int KeKo, float pdx, float ptol, int ynbi, int ynbs, char *nmmod, YREAL yinput[], YREAL ystate[], locate_d *Modul)
{	 int wk, wj, nbKo=0; YREAL wSj[YMAX_NBS],wdxk,wJjk,wyj,wDjk,wdifjk;
	 YREAL svx; int resKo; char okko[2][3]={"ok", "KO"}; 
	 if (!Ydftestijkt(Yid_locate_d)) return(0); 
	 memset (Yjac, 0, YMAX_JAC_NBI*YMAX_JAC_NBS*sizeof(YREAL));
/*M01_GTESTT*/
	 Modul->forward( Yting[0], Yting[1] );
	 for (wj=0; wj<ynbs; ++wj) wSj[wj] = ystate[YTemps*ynbs + wj];
	 Modul->backward( Yting[0], Yting[1] );
	 for (wk=0; wk<ynbi; ++wk)
	 {	 svx=yinput[wk];
			 wdxk = svx*pdx;
			 if (modop||svx==0.0) wdxk = pdx; 
			 yinput[wk] += wdxk;
			 Modul->forward( Yting[0], Yting[1] );
			 for (wj=0; wj<ynbs; ++wj)
			 { wJjk = Yjac[wj][wk];
				 wyj = ystate[YTemps*ynbs + wj];
				 wDjk = (wyj - wSj[wj]) / wdxk;
				 wdifjk = fabs(wDjk - wJjk);
				 resKo=(wdifjk/fabs(wJjk))>ptol;  nbKo+=resKo;
				 if(resKo || !KeKo) 
				 printf(" %7s | J[%2i,%2i] |% 12.6e |% 12.6e |% 12.6e :  %s |% 12.6e | (%i %i %i %i)\n",
								 nmmod, wj+1, wk+1, wJjk, wDjk, wdifjk, okko[resKo], svx, Yi+1, Yj+1, Yk+1, Yt+1);
				} 
			  yinput[wk] = svx;
	 }
	 return(nbKo);
}/*M3_CREU*/ 
void Ycreate_locate_d()  
{
	 for (int Yw1=0; Yw1<YA1_S_eulerlocate_d; ++Yw1){
		 for (int Yw2=0; Yw2<YA2_S_eulerlocate_d; ++Yw2){
			 for (int Yw3=0; Yw3<YA3_S_eulerlocate_d; ++Yw3){
				 if ((Ylocate_d(Yw1, Yw2, Yw3) = new locate_d ())==NULL) Ycreate_errnew();}}}
}/*M1_TBTOGT  */ 
void YtabTOgrad_r_cout_d(int frompdt, int topdt, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_r_cout_d; ++Yws){
				 YG_r_cout_d(Yws, Yw1, Ywt) = tab[Y3windice]; ++Y3windice;}}}
}/*M1_GEQPDT*/ 
void YgradEQPdelta_r_cout_d(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_r_cout_d; ++Yws){
				 YG_r_cout_d(Yws, Yw1, Ywt) = YD_r_cout_d(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_DEQPST*/ 
void YdeltaEQPstate_r_cout_d(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_r_cout_d; ++Yws){
				 YD_r_cout_d(Yws, Yw1, Ywt) = YS_r_cout_d(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_DEQPGT*/ 
void YdeltaEQPgrad_r_cout_d(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_r_cout_d; ++Yws){
				 YD_r_cout_d(Yws, Yw1, Ywt) = YG_r_cout_d(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_SEQPDT*/ 
void YstateEQPdelta_r_cout_d(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_r_cout_d; ++Yws){
				 YS_r_cout_d(Yws, Yw1, Ywt) = YD_r_cout_d(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_SEQPODT*/ 
void YstateEQPOdelta_r_cout_d(int frompdt, int topdt, double pfact, char *codop)
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_r_cout_d; ++Yws){
				 YS_r_cout_d(Yws, Yw1, Ywt) = YD_r_cout_d(Yws, Yw1, Ywt) * (pfact*(1+Ygetval(codop)));}}}
}/*M1_SEQAPTDT*/ 
void YstateEQAPTdelta_r_cout_d(int frompdt, int topdt, double pfact, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_r_cout_d; ++Yws){
				 YS_r_cout_d(Yws, Yw1, Ywt) = YD_r_cout_d(Yws, Yw1, Ywt) + pfact * tab[Y3windice++];}}}
}/*M1_GEQPST*/ 
void YgradEQPstate_r_cout_d(int frompdt, int topdt, double pfact)
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_r_cout_d; ++Yws){
				 YG_r_cout_d(Yws, Yw1, Ywt) = YS_r_cout_d(Yws, Yw1, Ywt) * pfact;}}}
}/*M1_GEQPOST*/ 
void YgradEQPOstate_r_cout_d(int frompdt, int topdt, double pfact, char *codop)
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_r_cout_d; ++Yws){
				 YG_r_cout_d(Yws, Yw1, Ywt) = YS_r_cout_d(Yws, Yw1, Ywt) * (pfact*(1+Ygetval(codop)));}}}
}/*M1_STOTBT  */ 
void YstateTOtab_r_cout_d(int frompdt, int topdt, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_r_cout_d; ++Yws){
				 tab[Y3windice] = YS_r_cout_d(Yws, Yw1, Ywt); ++Y3windice;}}}
}/*M1_GTOTBT  */ 
void YgradTOtab_r_cout_d(int frompdt, int topdt, YREAL tab[])
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=frompdt; Ywt<topdt;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_r_cout_d; ++Yws){
				 tab[Y3windice] = YG_r_cout_d(Yws, Yw1, Ywt); ++Y3windice;}}}
}/*M1_RGRADT*/ 
void Yrazgrad_r_cout_d()
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=0; Ywt<YNBALLTIME_T_float_d; ++Ywt){
			 for (int Yws=0; Yws<YNBS_r_cout_d; ++Yws){
				 Yr_cout_d(Yw1)->Ygrad[Ywt][Yws] = 0.0;}}}
}/*M1_SETAT*/ 
void Ysetstate_r_cout_d(double val)
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Ywt=0; Ywt<YNBALLTIME_T_float_d;  ++Ywt){
			 for (int Yws=0; Yws<YNBS_r_cout_d; ++Yws){
				 Yr_cout_d(Yw1)->Ystate[Ywt][Yws] = val;}}}
}/*M1_YIOUT*/ 
void Yio_r_cout_d(int Yws, int Yw1, int Yw2, int Yw3, int Ywt, double val)	
{  Yiorwf ("r_cout_d", Yws, Yw1, Yw2, Yw3, Ywt, &YS_r_cout_d(Yws, Yw1, Ywt), val);
}/*M00_GTEST */
int Ytestdf_r_cout_d (int modop, int KeKo, float pdx, float ptol, int ynbi, int ynbs, char *nmmod, YREAL yinput[], YREAL ystate[], r_cout_d *Modul)
{	 int wk, wj, nbKo=0; YREAL wSj[YMAX_NBS],wdxk,wJjk,wyj,wDjk,wdifjk;
	 YREAL svx; int resKo; char okko[2][3]={"ok", "KO"}; 
	 if (!Ydftestijkt(Yid_r_cout_d)) return(0); 
	 memset (Yjac, 0, YMAX_JAC_NBI*YMAX_JAC_NBS*sizeof(YREAL));
/*M01_GTESTT*/
	 Modul->forward( Yting[0], Yting[1] );
	 for (wj=0; wj<ynbs; ++wj) wSj[wj] = ystate[YTemps*ynbs + wj];
	 Modul->backward( Yting[0], Yting[1] );
	 for (wk=0; wk<ynbi; ++wk)
	 {	 svx=yinput[wk];
			 wdxk = svx*pdx;
			 if (modop||svx==0.0) wdxk = pdx; 
			 yinput[wk] += wdxk;
			 Modul->forward( Yting[0], Yting[1] );
			 for (wj=0; wj<ynbs; ++wj)
			 { wJjk = Yjac[wj][wk];
				 wyj = ystate[YTemps*ynbs + wj];
				 wDjk = (wyj - wSj[wj]) / wdxk;
				 wdifjk = fabs(wDjk - wJjk);
				 resKo=(wdifjk/fabs(wJjk))>ptol;  nbKo+=resKo;
				 if(resKo || !KeKo) 
				 printf(" %7s | J[%2i,%2i] |% 12.6e |% 12.6e |% 12.6e :  %s |% 12.6e | (%i %i %i %i)\n",
								 nmmod, wj+1, wk+1, wJjk, wDjk, wdifjk, okko[resKo], svx, Yi+1, Yj+1, Yk+1, Yt+1);
				} 
			  yinput[wk] = svx;
	 }
	 return(nbKo);
}/*M0_OUTOB*/ 
int Youtoobs_r_cout_d(YioKind yiokind, int Yws, int Ywt, int arbpdt)
{
	 Yaobs.time = arbpdt;
	 Yaobs.imod = Yimod("r_cout_d");
	 Yaobs.smod = Yws;/*M1_OUTOBT*/
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1)
	 {	if(select_io(yiokind, "r_cout_d", Yws, Yw1, -1, -1, Ywt, &Ystate_r_cout_d(Yws, Yw1, Ywt)))
			{	Yaobs.iaxe = Yw1; Yaobs.jaxe = -1; Yaobs.kaxe = -1;
				Yaobs.vobs = Ystate_r_cout_d(Yws, Yw1, Ywt);
				if(!Yobs_insert (&Yaobs)) return(0);
			}
	 }/*M9_OUTOB*/  
	 return(1);
}/*M1_CREU*/ 
void Ycreate_r_cout_d()  
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 if ((Yr_cout_d(Yw1) = new r_cout_d ())==NULL) Ycreate_errnew();}
}/*M1_SWISH*/ 
void Ysetwish_r_cout_d(double val)
{
	 for (int Yw1=0; Yw1<YA1_S_float_d; ++Yw1){
		 for (int Yws=0; Yws<YNBS_r_cout_d; ++Yws){
			 Ywish_r_cout_d(Yws, Yw1) = val;}}
}
void Ycreate_all () 
{
	 Ycreate_r_float ();
	 Ycreate_ur ();
	 Ycreate_uinter ();
	 Ycreate_xinter ();
	 Ycreate_u ();
	 Ycreate_v ();
	 Ycreate_locate ();
	 Ycreate_r_float_d ();
	 Ycreate_ur_d ();
	 Ycreate_uinter_d ();
	 Ycreate_u_d ();
	 Ycreate_v_d ();
	 Ycreate_locate_d ();
	 Ycreate_r_cout_d ();
}
void Ywishdiff_all (char *nmmod, int Yws, int Yw1, int Yw2, int Yw3, int Ywt, YREAL vobs) 
{	
	 if(1==0);
	 else if (strcmp(nmmod, "u_d") == 0)
				 YW_u_d(Yws, Yw1, Yw2) = YS_u_d(Yws, Yw1, Yw2) - vobs;
	 else if (strcmp(nmmod, "v_d") == 0)
				 YW_v_d(Yws, Yw1, Yw2) = YS_v_d(Yws, Yw1, Yw2) - vobs;
	 else if (strcmp(nmmod, "r_cout_d") == 0)
				 YW_r_cout_d(Yws, Yw1) = YS_r_cout_d(Yws, Yw1, Ywt) - vobs;
}
void Ycostwishdiff_all (char *nmmod, int Yws, int Yw1, int Yw2, int Yw3, int Ywt) 
{	
	 if(1==0);
	 else if (strcmp(nmmod, "u_d") == 0)
				 YG_u_d(Yws, Yw1, Yw2) += Ycostdiff(YW_u_d(Yws, Yw1, Yw2), YW_u_d(Yws, Yw1, Yw2));
	 else if (strcmp(nmmod, "v_d") == 0)
				 YG_v_d(Yws, Yw1, Yw2) += Ycostdiff(YW_v_d(Yws, Yw1, Yw2), YW_v_d(Yws, Yw1, Yw2));
	 else if (strcmp(nmmod, "r_cout_d") == 0)
				 YG_r_cout_d(Yws, Yw1, Ywt) += Ycostdiff(YW_r_cout_d(Yws, Yw1), YW_r_cout_d(Yws, Yw1));
}
void Youtoobs_mod(YioKind yiokind, char *nmmod, int numout, int pdt, int arbpdt)
{
 	 if(1==0);
 	 else if (strcmp(nmmod, "u_d") == 0) Youtoobs_u_d(yiokind, numout, pdt, arbpdt);
 	 else if (strcmp(nmmod, "v_d") == 0) Youtoobs_v_d(yiokind, numout, pdt, arbpdt);
 	 else if (strcmp(nmmod, "r_cout_d") == 0) Youtoobs_r_cout_d(yiokind, numout, pdt, arbpdt);
}
void Yadjust_all () 
{
	 Yadjust_u_d ();
	 Yadjust_v_d ();
}
void Yrazgrad_all () 
{
	 Yrazgrad_r_float ();
	 Yrazgrad_ur ();
	 Yrazgrad_uinter ();
	 Yrazgrad_xinter ();
	 Yrazgrad_u ();
	 Yrazgrad_v ();
	 Yrazgrad_locate ();
	 Yrazgrad_r_float_d ();
	 Yrazgrad_ur_d ();
	 Yrazgrad_uinter_d ();
	 Yrazgrad_u_d ();
	 Yrazgrad_v_d ();
	 Yrazgrad_locate_d ();
	 Yrazgrad_r_cout_d ();
}
void Yrazgrad_only (int itraj) 
{
	if (itraj==Yid_T_float)
	{}
	if (itraj==Yid_T_euler)
	{}
	if (itraj==Yid_T_float_d)
	{}
	if (itraj==Yid_T_euler_d)
	{}
}
int Ysetstate_mod (char *nmmod, YREAL val) 
{	int codret=0; int all=0;
	if (!strcmp(nmmod, "Y#A")) all=1;
	if (!strcmp(nmmod, "r_float") || all)
	{  Ysetstate_r_float(val);codret=1;}
	if (!strcmp(nmmod, "ur") || all)
	{  Ysetstate_ur(val);codret=1;}
	if (!strcmp(nmmod, "uinter") || all)
	{  Ysetstate_uinter(val);codret=1;}
	if (!strcmp(nmmod, "xinter") || all)
	{  Ysetstate_xinter(val);codret=1;}
	if (!strcmp(nmmod, "u") || all)
	{  Ysetstate_u(val);codret=1;}
	if (!strcmp(nmmod, "v") || all)
	{  Ysetstate_v(val);codret=1;}
	if (!strcmp(nmmod, "locate") || all)
	{  Ysetstate_locate(val);codret=1;}
	if (!strcmp(nmmod, "r_float_d") || all)
	{  Ysetstate_r_float_d(val);codret=1;}
	if (!strcmp(nmmod, "ur_d") || all)
	{  Ysetstate_ur_d(val);codret=1;}
	if (!strcmp(nmmod, "uinter_d") || all)
	{  Ysetstate_uinter_d(val);codret=1;}
	if (!strcmp(nmmod, "u_d") || all)
	{  Ysetstate_u_d(val);codret=1;}
	if (!strcmp(nmmod, "v_d") || all)
	{  Ysetstate_v_d(val);codret=1;}
	if (!strcmp(nmmod, "locate_d") || all)
	{  Ysetstate_locate_d(val);codret=1;}
	if (!strcmp(nmmod, "r_cout_d") || all)
	{  Ysetstate_r_cout_d(val);codret=1;}
	return(codret);
}
void Ysetwish_mod (int imod, YREAL val) 
{	if(1==0);
	else if (imod==10) Ysetwish_u_d(val);
	else if (imod==11) Ysetwish_v_d(val);
	else if (imod==13) Ysetwish_r_cout_d(val);
}
void  Yio_mod (char *nmmod, int Yws, int Yw1, int Yw2, int Yw3, int Ywt, double val)
{
	 if(1==0);
	 else if (strcmp(nmmod, "r_float") == 0) Yio_r_float(Yws, Yw1, Yw2, Yw3, Ywt, val);
	 else if (strcmp(nmmod, "ur") == 0) Yio_ur(Yws, Yw1, Yw2, Yw3, Ywt, val);
	 else if (strcmp(nmmod, "uinter") == 0) Yio_uinter(Yws, Yw1, Yw2, Yw3, Ywt, val);
	 else if (strcmp(nmmod, "xinter") == 0) Yio_xinter(Yws, Yw1, Yw2, Yw3, Ywt, val);
	 else if (strcmp(nmmod, "u") == 0) Yio_u(Yws, Yw1, Yw2, Yw3, Ywt, val);
	 else if (strcmp(nmmod, "v") == 0) Yio_v(Yws, Yw1, Yw2, Yw3, Ywt, val);
	 else if (strcmp(nmmod, "locate") == 0) Yio_locate(Yws, Yw1, Yw2, Yw3, Ywt, val);
	 else if (strcmp(nmmod, "r_float_d") == 0) Yio_r_float_d(Yws, Yw1, Yw2, Yw3, Ywt, val);
	 else if (strcmp(nmmod, "ur_d") == 0) Yio_ur_d(Yws, Yw1, Yw2, Yw3, Ywt, val);
	 else if (strcmp(nmmod, "uinter_d") == 0) Yio_uinter_d(Yws, Yw1, Yw2, Yw3, Ywt, val);
	 else if (strcmp(nmmod, "u_d") == 0) Yio_u_d(Yws, Yw1, Yw2, Yw3, Ywt, val);
	 else if (strcmp(nmmod, "v_d") == 0) Yio_v_d(Yws, Yw1, Yw2, Yw3, Ywt, val);
	 else if (strcmp(nmmod, "locate_d") == 0) Yio_locate_d(Yws, Yw1, Yw2, Yw3, Ywt, val);
	 else if (strcmp(nmmod, "r_cout_d") == 0) Yio_r_cout_d(Yws, Yw1, Yw2, Yw3, Ywt, val);
}
void YgradCTOtab_target (YREAL tab[]) 
{	//Y3windice = 0;
	 YgradCTOtab_u_d (tab);
	 YgradCTOtab_v_d (tab);
}
void YtabTOgrad_target (YREAL tab[]) 
{	//Y3windice = 0;
	 YtabTOgrad_u_d (tab);
	 YtabTOgrad_v_d (tab);
}
void Y3getstate_all (float x[]) 
{	 Y3windice = 0;
	 Y3getstate_u_d (x);
	 Y3getstate_v_d (x);
}
void Y3valstate_all () 
{	 Y3windice = 0;
	 Y3valstate_u_d ();
	 Y3valstate_v_d ();
}
void Y3valgrad_all (float g[]) 
{	 Y3windice = 0;
	 Y3valgrad_u_d (g);
	 Y3valgrad_v_d (g);
}
void YdeltaEQPstate_traj(int itraj, char *norkmod, int frompdt, int topdt, double pfact) 
{
	if (1==0){}
	else if (itraj==Yid_T_float)
	{
		if ( !strcmp(norkmod, "r_float") || !strcmp(norkmod, "Y#A") )
	 YdeltaEQPstate_r_float (frompdt, topdt, pfact);
		if ( !strcmp(norkmod, "ur") || !strcmp(norkmod, "Y#A") )
	 YdeltaEQPstate_ur (frompdt, topdt, pfact);
		if ( !strcmp(norkmod, "uinter") || !strcmp(norkmod, "Y#A") )
	 YdeltaEQPstate_uinter (frompdt, topdt, pfact);
		if ( !strcmp(norkmod, "xinter") || !strcmp(norkmod, "Y#A") )
	 YdeltaEQPstate_xinter (frompdt, topdt, pfact);
		if ( !strcmp(norkmod, "locate") || !strcmp(norkmod, "Y#A") )
	 YdeltaEQPstate_locate (frompdt, topdt, pfact);
	}
	else if (itraj==Yid_T_euler)
	{
		if ( !strcmp(norkmod, "u") || !strcmp(norkmod, "Y#A") )
	 YdeltaEQPstate_u (0, 1, pfact);
		if ( !strcmp(norkmod, "v") || !strcmp(norkmod, "Y#A") )
	 YdeltaEQPstate_v (0, 1, pfact);
	}
	else if (itraj==Yid_T_float_d)
	{
		if ( !strcmp(norkmod, "r_float_d") || !strcmp(norkmod, "Y#A") )
	 YdeltaEQPstate_r_float_d (frompdt, topdt, pfact);
		if ( !strcmp(norkmod, "ur_d") || !strcmp(norkmod, "Y#A") )
	 YdeltaEQPstate_ur_d (frompdt, topdt, pfact);
		if ( !strcmp(norkmod, "uinter_d") || !strcmp(norkmod, "Y#A") )
	 YdeltaEQPstate_uinter_d (frompdt, topdt, pfact);
		if ( !strcmp(norkmod, "locate_d") || !strcmp(norkmod, "Y#A") )
	 YdeltaEQPstate_locate_d (frompdt, topdt, pfact);
		if ( !strcmp(norkmod, "r_cout_d") || !strcmp(norkmod, "Y#A") || !strcmp(norkmod, "Y#C") )
	 YdeltaEQPstate_r_cout_d (frompdt, topdt, pfact);
	}
	else if (itraj==Yid_T_euler_d)
	{
		if ( !strcmp(norkmod, "u_d") || !strcmp(norkmod, "Y#A") || !strcmp(norkmod, "Y#T") )
	 YdeltaEQPstate_u_d (0, 1, pfact);
		if ( !strcmp(norkmod, "v_d") || !strcmp(norkmod, "Y#A") || !strcmp(norkmod, "Y#T") )
	 YdeltaEQPstate_v_d (0, 1, pfact);
	}
}
void YdeltaEQPstate_target(double pfact) 
{
	 YdeltaEQPstate_u_d (0, 1, pfact);
	 YdeltaEQPstate_v_d (0, 1, pfact);
}
void YdeltaEQPgrad_target(double pfact) 
{
	 YdeltaEQPgrad_u_d (0, 1, pfact);
	 YdeltaEQPgrad_v_d (0, 1, pfact);
}
void YstateEQPdelta_traj(int itraj, char *norkmod, int frompdt, int topdt, double pfact) 
{
	if (1==0){}
	else if (itraj==Yid_T_float)
	{
		if ( !strcmp(norkmod, "r_float") || !strcmp(norkmod, "Y#A") )
	 YstateEQPdelta_r_float (frompdt, topdt, pfact);
		if ( !strcmp(norkmod, "ur") || !strcmp(norkmod, "Y#A") )
	 YstateEQPdelta_ur (frompdt, topdt, pfact);
		if ( !strcmp(norkmod, "uinter") || !strcmp(norkmod, "Y#A") )
	 YstateEQPdelta_uinter (frompdt, topdt, pfact);
		if ( !strcmp(norkmod, "xinter") || !strcmp(norkmod, "Y#A") )
	 YstateEQPdelta_xinter (frompdt, topdt, pfact);
		if ( !strcmp(norkmod, "locate") || !strcmp(norkmod, "Y#A") )
	 YstateEQPdelta_locate (frompdt, topdt, pfact);
	}
	else if (itraj==Yid_T_euler)
	{
		if ( !strcmp(norkmod, "u") || !strcmp(norkmod, "Y#A") )
	 YstateEQPdelta_u (0, 1, pfact);
		if ( !strcmp(norkmod, "v") || !strcmp(norkmod, "Y#A") )
	 YstateEQPdelta_v (0, 1, pfact);
	}
	else if (itraj==Yid_T_float_d)
	{
		if ( !strcmp(norkmod, "r_float_d") || !strcmp(norkmod, "Y#A") )
	 YstateEQPdelta_r_float_d (frompdt, topdt, pfact);
		if ( !strcmp(norkmod, "ur_d") || !strcmp(norkmod, "Y#A") )
	 YstateEQPdelta_ur_d (frompdt, topdt, pfact);
		if ( !strcmp(norkmod, "uinter_d") || !strcmp(norkmod, "Y#A") )
	 YstateEQPdelta_uinter_d (frompdt, topdt, pfact);
		if ( !strcmp(norkmod, "locate_d") || !strcmp(norkmod, "Y#A") )
	 YstateEQPdelta_locate_d (frompdt, topdt, pfact);
		if ( !strcmp(norkmod, "r_cout_d") || !strcmp(norkmod, "Y#A") || !strcmp(norkmod, "Y#C") )
	 YstateEQPdelta_r_cout_d (frompdt, topdt, pfact);
	}
	else if (itraj==Yid_T_euler_d)
	{
		if ( !strcmp(norkmod, "u_d") || !strcmp(norkmod, "Y#A") || !strcmp(norkmod, "Y#T") )
	 YstateEQPdelta_u_d (0, 1, pfact);
		if ( !strcmp(norkmod, "v_d") || !strcmp(norkmod, "Y#A") || !strcmp(norkmod, "Y#T") )
	 YstateEQPdelta_v_d (0, 1, pfact);
	}
}
void YstateEQPdelta_target(double pfact) 
{
	 YstateEQPdelta_u_d (0, 1, pfact);
	 YstateEQPdelta_v_d (0, 1, pfact);
}
void YstateEQPOdelta_target(double pfact, char *codop) 
{
	 YstateEQPOdelta_u_d (0, 1, pfact, codop);
	 YstateEQPOdelta_v_d (0, 1, pfact, codop);
}
void YstateEQAPTdelta_target(double pfact, YREAL tab[]) 
{
	 YstateEQAPTdelta_u_d (0, 1, pfact, tab);
	 YstateEQAPTdelta_v_d (0, 1, pfact, tab);
}
void YgradEQPstate_traj(int itraj, char *norkmod, int frompdt, int topdt, double pfact) 
{
	if (1==0){}
	else if (itraj==Yid_T_float)
	{
		if ( !strcmp(norkmod, "r_float") || !strcmp(norkmod, "Y#A") )
	 YgradEQPstate_r_float (frompdt, topdt, pfact);
		if ( !strcmp(norkmod, "ur") || !strcmp(norkmod, "Y#A") )
	 YgradEQPstate_ur (frompdt, topdt, pfact);
		if ( !strcmp(norkmod, "uinter") || !strcmp(norkmod, "Y#A") )
	 YgradEQPstate_uinter (frompdt, topdt, pfact);
		if ( !strcmp(norkmod, "xinter") || !strcmp(norkmod, "Y#A") )
	 YgradEQPstate_xinter (frompdt, topdt, pfact);
		if ( !strcmp(norkmod, "locate") || !strcmp(norkmod, "Y#A") )
	 YgradEQPstate_locate (frompdt, topdt, pfact);
	}
	else if (itraj==Yid_T_euler)
	{
		if ( !strcmp(norkmod, "u") || !strcmp(norkmod, "Y#A") )
	 YgradEQPstate_u(pfact);
		if ( !strcmp(norkmod, "v") || !strcmp(norkmod, "Y#A") )
	 YgradEQPstate_v(pfact);
	}
	else if (itraj==Yid_T_float_d)
	{
		if ( !strcmp(norkmod, "r_float_d") || !strcmp(norkmod, "Y#A") )
	 YgradEQPstate_r_float_d (frompdt, topdt, pfact);
		if ( !strcmp(norkmod, "ur_d") || !strcmp(norkmod, "Y#A") )
	 YgradEQPstate_ur_d (frompdt, topdt, pfact);
		if ( !strcmp(norkmod, "uinter_d") || !strcmp(norkmod, "Y#A") )
	 YgradEQPstate_uinter_d (frompdt, topdt, pfact);
		if ( !strcmp(norkmod, "locate_d") || !strcmp(norkmod, "Y#A") )
	 YgradEQPstate_locate_d (frompdt, topdt, pfact);
		if ( !strcmp(norkmod, "r_cout_d") || !strcmp(norkmod, "Y#A") || !strcmp(norkmod, "Y#C") )
	 YgradEQPstate_r_cout_d (frompdt, topdt, pfact);
	}
	else if (itraj==Yid_T_euler_d)
	{
		if ( !strcmp(norkmod, "u_d") || !strcmp(norkmod, "Y#A") || !strcmp(norkmod, "Y#T") )
	 YgradEQPstate_u_d(pfact);
		if ( !strcmp(norkmod, "v_d") || !strcmp(norkmod, "Y#A") || !strcmp(norkmod, "Y#T") )
	 YgradEQPstate_v_d(pfact);
	}
}
void YgradEQPstate_target(double pfact) 
{
	 YgradEQPstate_u_d (pfact);
	 YgradEQPstate_v_d (pfact);
}
void YgradEQPOstate_target(double pfact, char *codop) 
{
	 YgradEQPOstate_u_d (pfact, codop);
	 YgradEQPOstate_v_d (pfact, codop);
}
void YstateTOtab_traj (int itraj, char *norkmod, int frompdt, int topdt, YREAL tab[]) 
{	//Y3windice = 0;
	if (1==0){}
	else if (itraj==Yid_T_float)
	{
		if ( !strcmp(norkmod, "r_float") || !strcmp(norkmod, "Y#A") )
	 YstateTOtab_r_float (frompdt, topdt, tab);
		if ( !strcmp(norkmod, "ur") || !strcmp(norkmod, "Y#A") )
	 YstateTOtab_ur (frompdt, topdt, tab);
		if ( !strcmp(norkmod, "uinter") || !strcmp(norkmod, "Y#A") )
	 YstateTOtab_uinter (frompdt, topdt, tab);
		if ( !strcmp(norkmod, "xinter") || !strcmp(norkmod, "Y#A") )
	 YstateTOtab_xinter (frompdt, topdt, tab);
		if ( !strcmp(norkmod, "locate") || !strcmp(norkmod, "Y#A") )
	 YstateTOtab_locate (frompdt, topdt, tab);
	}
	else if (itraj==Yid_T_euler)
	{
		if ( !strcmp(norkmod, "u") || !strcmp(norkmod, "Y#A") )
	 YstateTOtab_u (tab);
		if ( !strcmp(norkmod, "v") || !strcmp(norkmod, "Y#A") )
	 YstateTOtab_v (tab);
	}
	else if (itraj==Yid_T_float_d)
	{
		if ( !strcmp(norkmod, "r_float_d") || !strcmp(norkmod, "Y#A") )
	 YstateTOtab_r_float_d (frompdt, topdt, tab);
		if ( !strcmp(norkmod, "ur_d") || !strcmp(norkmod, "Y#A") )
	 YstateTOtab_ur_d (frompdt, topdt, tab);
		if ( !strcmp(norkmod, "uinter_d") || !strcmp(norkmod, "Y#A") )
	 YstateTOtab_uinter_d (frompdt, topdt, tab);
		if ( !strcmp(norkmod, "locate_d") || !strcmp(norkmod, "Y#A") )
	 YstateTOtab_locate_d (frompdt, topdt, tab);
		if ( !strcmp(norkmod, "r_cout_d") || !strcmp(norkmod, "Y#A") || !strcmp(norkmod, "Y#C") )
	 YstateTOtab_r_cout_d (frompdt, topdt, tab);
	}
	else if (itraj==Yid_T_euler_d)
	{
		if ( !strcmp(norkmod, "u_d") || !strcmp(norkmod, "Y#A") || !strcmp(norkmod, "Y#T") )
	 YstateTOtab_u_d (tab);
		if ( !strcmp(norkmod, "v_d") || !strcmp(norkmod, "Y#A") || !strcmp(norkmod, "Y#T") )
	 YstateTOtab_v_d (tab);
	}
}
void YstateTOtab_target (YREAL tab[]) 
{	//Y3windice = 0;
	 YstateTOtab_u_d (tab);
	 YstateTOtab_v_d (tab);
}
void YgradTOtab_traj (int itraj, char *norkmod, int frompdt, int topdt, YREAL tab[]) 
{	//Y3windice = 0;
	if (1==0){}
	else if (itraj==Yid_T_float)
	{
		if ( !strcmp(norkmod, "r_float") || !strcmp(norkmod, "Y#A") )
	 YgradTOtab_r_float (frompdt, topdt, tab);
		if ( !strcmp(norkmod, "ur") || !strcmp(norkmod, "Y#A") )
	 YgradTOtab_ur (frompdt, topdt, tab);
		if ( !strcmp(norkmod, "uinter") || !strcmp(norkmod, "Y#A") )
	 YgradTOtab_uinter (frompdt, topdt, tab);
		if ( !strcmp(norkmod, "xinter") || !strcmp(norkmod, "Y#A") )
	 YgradTOtab_xinter (frompdt, topdt, tab);
		if ( !strcmp(norkmod, "locate") || !strcmp(norkmod, "Y#A") )
	 YgradTOtab_locate (frompdt, topdt, tab);
	}
	else if (itraj==Yid_T_euler)
	{
		if ( !strcmp(norkmod, "u") || !strcmp(norkmod, "Y#A") )
	 YgradTOtab_u (tab);
		if ( !strcmp(norkmod, "v") || !strcmp(norkmod, "Y#A") )
	 YgradTOtab_v (tab);
	}
	else if (itraj==Yid_T_float_d)
	{
		if ( !strcmp(norkmod, "r_float_d") || !strcmp(norkmod, "Y#A") )
	 YgradTOtab_r_float_d (frompdt, topdt, tab);
		if ( !strcmp(norkmod, "ur_d") || !strcmp(norkmod, "Y#A") )
	 YgradTOtab_ur_d (frompdt, topdt, tab);
		if ( !strcmp(norkmod, "uinter_d") || !strcmp(norkmod, "Y#A") )
	 YgradTOtab_uinter_d (frompdt, topdt, tab);
		if ( !strcmp(norkmod, "locate_d") || !strcmp(norkmod, "Y#A") )
	 YgradTOtab_locate_d (frompdt, topdt, tab);
		if ( !strcmp(norkmod, "r_cout_d") || !strcmp(norkmod, "Y#A") || !strcmp(norkmod, "Y#C") )
	 YgradTOtab_r_cout_d (frompdt, topdt, tab);
	}
	else if (itraj==Yid_T_euler_d)
	{
		if ( !strcmp(norkmod, "u_d") || !strcmp(norkmod, "Y#A") || !strcmp(norkmod, "Y#T") )
	 YgradTOtab_u_d (tab);
		if ( !strcmp(norkmod, "v_d") || !strcmp(norkmod, "Y#A") || !strcmp(norkmod, "Y#T") )
	 YgradTOtab_v_d (tab);
	}
}
void YgradTOtab_target (YREAL tab[]) 
{	//Y3windice = 0;
	 YgradTOtab_u_d (tab);
	 YgradTOtab_v_d (tab);
}
void YtabTOgrad_traj (int itraj, char *norkmod, int frompdt, int topdt, YREAL tab[]) 
{	//Y3windice = 0;
	if (1==0){}
	else if (itraj==Yid_T_float)
	{
		if ( !strcmp(norkmod, "r_float") || !strcmp(norkmod, "Y#A") )
			 YtabTOgrad_r_float (frompdt, topdt, tab);
		if ( !strcmp(norkmod, "ur") || !strcmp(norkmod, "Y#A") )
			 YtabTOgrad_ur (frompdt, topdt, tab);
		if ( !strcmp(norkmod, "uinter") || !strcmp(norkmod, "Y#A") )
			 YtabTOgrad_uinter (frompdt, topdt, tab);
		if ( !strcmp(norkmod, "xinter") || !strcmp(norkmod, "Y#A") )
			 YtabTOgrad_xinter (frompdt, topdt, tab);
		if ( !strcmp(norkmod, "locate") || !strcmp(norkmod, "Y#A") )
			 YtabTOgrad_locate (frompdt, topdt, tab);
	}
	else if (itraj==Yid_T_euler)
	{
		if ( !strcmp(norkmod, "u") || !strcmp(norkmod, "Y#A") )
			 YtabTOgrad_u (tab);
		if ( !strcmp(norkmod, "v") || !strcmp(norkmod, "Y#A") )
			 YtabTOgrad_v (tab);
	}
	else if (itraj==Yid_T_float_d)
	{
		if ( !strcmp(norkmod, "r_float_d") || !strcmp(norkmod, "Y#A") )
			 YtabTOgrad_r_float_d (frompdt, topdt, tab);
		if ( !strcmp(norkmod, "ur_d") || !strcmp(norkmod, "Y#A") )
			 YtabTOgrad_ur_d (frompdt, topdt, tab);
		if ( !strcmp(norkmod, "uinter_d") || !strcmp(norkmod, "Y#A") )
			 YtabTOgrad_uinter_d (frompdt, topdt, tab);
		if ( !strcmp(norkmod, "locate_d") || !strcmp(norkmod, "Y#A") )
			 YtabTOgrad_locate_d (frompdt, topdt, tab);
		if ( !strcmp(norkmod, "r_cout_d") || !strcmp(norkmod, "Y#A") || !strcmp(norkmod, "Y#C") )
			 YtabTOgrad_r_cout_d (frompdt, topdt, tab);
	}
	else if (itraj==Yid_T_euler_d)
	{
		if ( !strcmp(norkmod, "u_d") || !strcmp(norkmod, "Y#A") || !strcmp(norkmod, "Y#T") )
			 YtabTOgrad_u_d (tab);
		if ( !strcmp(norkmod, "v_d") || !strcmp(norkmod, "Y#A") || !strcmp(norkmod, "Y#T") )
			 YtabTOgrad_v_d (tab);
	}
}
 int Ydfward_all(int modop, char *nmmod, int All, int KeKo,  float pdx, float ptol)
{
	int nbko=0;
	if ( !strcmp(nmmod, "r_float") || All )
		nbko += Ytestdf_r_float(modop,KeKo,pdx,ptol,YNBI_r_float,YNBS_r_float,"r_float",Yting,Yr_float[Yi]->Ystate[0],Yr_float[Yi]);
	if ( !strcmp(nmmod, "ur") || All )
		nbko += Ytestdf_ur(modop,KeKo,pdx,ptol,YNBI_ur,YNBS_ur,"ur",Yting,Yur[Yi]->Ystate[0],Yur[Yi]);
	if ( !strcmp(nmmod, "locate") || All )
		nbko += Ytestdf_locate(modop,KeKo,pdx,ptol,YNBI_locate,YNBS_locate,"locate",Yting,Ylocate[Yi][Yj][Yk]->Ystate[0],Ylocate[Yi][Yj][Yk]);
	if ( !strcmp(nmmod, "r_float_d") || All )
		nbko += Ytestdf_r_float_d(modop,KeKo,pdx,ptol,YNBI_r_float_d,YNBS_r_float_d,"r_float_d",Yting,Yr_float_d[Yi]->Ystate[0],Yr_float_d[Yi]);
	if ( !strcmp(nmmod, "ur_d") || All )
		nbko += Ytestdf_ur_d(modop,KeKo,pdx,ptol,YNBI_ur_d,YNBS_ur_d,"ur_d",Yting,Yur_d[Yi]->Ystate[0],Yur_d[Yi]);
	if ( !strcmp(nmmod, "locate_d") || All )
		nbko += Ytestdf_locate_d(modop,KeKo,pdx,ptol,YNBI_locate_d,YNBS_locate_d,"locate_d",Yting,Ylocate_d[Yi][Yj][Yk]->Ystate[0],Ylocate_d[Yi][Yj][Yk]);
	if ( !strcmp(nmmod, "r_cout_d") || All )
		nbko += Ytestdf_r_cout_d(modop,KeKo,pdx,ptol,YNBI_r_cout_d,YNBS_r_cout_d,"r_cout_d",Yting,Yr_cout_d[Yi]->Ystate[0],Yr_cout_d[Yi]);
	return(nbko);
}

void YwishEQPstate_traj_tocov (int itraj, int pdt, double pfact) 
{
	if (1==0){}
	else if (itraj==Yid_T_float)
	{
	}
	else if (itraj==Yid_T_euler)
	{
	}
	else if (itraj==Yid_T_float_d)
	{
	}
	else if (itraj==Yid_T_euler_d)
	{
	}
}
int Yauto_call (int argc, char *argv[]) 
{
	 int codret=0; int all=0; double val;
	 if (1==0);
	 else if ( !strcmp(argv[0], "setepsi") || !strcmp(argv[0], "SETEPSI") 
					 ||!strcmp(argv[0], "setepsi_all") || !strcmp(argv[0], "SETEPSI_ALL"))
	 {
			 if (!strcmp(argv[0], "setepsi_all") || !strcmp(argv[0], "SETEPSI_ALL")){all=1; val=atof(argv[1]);} else val=atof(argv[2]);
			 if (!strcmp(argv[1], "u_d") || all)
					 {Ysetepsi_u_d(val);codret=1;}
			 if (!strcmp(argv[1], "v_d") || all)
					 {Ysetepsi_v_d(val);codret=1;}
	 }
	 else codret=0;
	 return(codret);
}

 int Yforward_space_S_float()
{
	 YA1=1; YA2=0; YA3=0;
Yi=-1; Yj=-1; Yk=-1; /* init des indices de maille: maj par la boucle si valide*/
for(Yi=0; Yi<YA1_S_float; ++Yi)
{
	 Yting[0]=YS1_uinter( Yi, YTemps);
	 Yting[1]=YS2_uinter( Yi, YTemps);
	 Yting[2]=YS3_uinter( Yi, YTemps);
	 Yting[3]=YS4_uinter( Yi, YTemps);
	 Yting[4]=YS5_uinter( Yi, YTemps);
	 Yting[5]=YS6_uinter( Yi, YTemps);
	 Yting[6]=YS7_uinter( Yi, YTemps);
	 Yting[7]=YS8_uinter( Yi, YTemps);
	 Yting[8]=YS1_xinter( Yi, YTemps);
	 Yting[9]=YS2_xinter( Yi, YTemps);
	 Yting[10]=YS3_xinter( Yi, YTemps);
	 Yting[11]=YS4_xinter( Yi, YTemps);
	 Yting[12]=YS5_xinter( Yi, YTemps);
	 Yting[13]=YS6_xinter( Yi, YTemps);
	 Yting[14]=YS7_xinter( Yi, YTemps);
	 Yting[15]=YS8_xinter( Yi, YTemps);
	 Yting[16]=YS1_r_float( Yi, YTemps-1);
	 Yting[17]=YS2_r_float( Yi, YTemps-1);
	 Yur(Yi)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4], Yting[5], Yting[6], Yting[7], Yting[8], Yting[9], Yting[10], Yting[11], Yting[12], Yting[13], Yting[14], Yting[15], Yting[16], Yting[17] );
	 Yting[0]=YS1_r_float( Yi, YTemps-1);
	 Yting[1]=YS2_r_float( Yi, YTemps-1);
	 Yting[2]=YS1_ur( Yi, YTemps);
	 Yting[3]=YS2_ur( Yi, YTemps);
	 Yr_float(Yi)->forward( Yting[0], Yting[1], Yting[2], Yting[3] );

}

	 return(0);
}

 int Ylinward_space_S_float()
{
	 YA1=1; YA2=0; YA3=0;
Yi=-1; Yj=-1; Yk=-1; /* init des indices de maille: maj par la boucle si valide*/
for(Yi=0; Yi<YA1_S_float; ++Yi)
{
	 Yting[0]=YS1_uinter( Yi, YTemps);
	 Yting[1]=YS2_uinter( Yi, YTemps);
	 Yting[2]=YS3_uinter( Yi, YTemps);
	 Yting[3]=YS4_uinter( Yi, YTemps);
	 Yting[4]=YS5_uinter( Yi, YTemps);
	 Yting[5]=YS6_uinter( Yi, YTemps);
	 Yting[6]=YS7_uinter( Yi, YTemps);
	 Yting[7]=YS8_uinter( Yi, YTemps);
	 Yting[8]=YS1_xinter( Yi, YTemps);
	 Yting[9]=YS2_xinter( Yi, YTemps);
	 Yting[10]=YS3_xinter( Yi, YTemps);
	 Yting[11]=YS4_xinter( Yi, YTemps);
	 Yting[12]=YS5_xinter( Yi, YTemps);
	 Yting[13]=YS6_xinter( Yi, YTemps);
	 Yting[14]=YS7_xinter( Yi, YTemps);
	 Yting[15]=YS8_xinter( Yi, YTemps);
	 Yting[16]=YS1_r_float( Yi, YTemps-1);
	 Yting[17]=YS2_r_float( Yi, YTemps-1);
	 Ytbeta[0]=YG1_uinter( Yi, YTemps);
	 Ytbeta[1]=YG2_uinter( Yi, YTemps);
	 Ytbeta[2]=YG3_uinter( Yi, YTemps);
	 Ytbeta[3]=YG4_uinter( Yi, YTemps);
	 Ytbeta[4]=YG5_uinter( Yi, YTemps);
	 Ytbeta[5]=YG6_uinter( Yi, YTemps);
	 Ytbeta[6]=YG7_uinter( Yi, YTemps);
	 Ytbeta[7]=YG8_uinter( Yi, YTemps);
	 Ytbeta[8]=YG1_xinter( Yi, YTemps);
	 Ytbeta[9]=YG2_xinter( Yi, YTemps);
	 Ytbeta[10]=YG3_xinter( Yi, YTemps);
	 Ytbeta[11]=YG4_xinter( Yi, YTemps);
	 Ytbeta[12]=YG5_xinter( Yi, YTemps);
	 Ytbeta[13]=YG6_xinter( Yi, YTemps);
	 Ytbeta[14]=YG7_xinter( Yi, YTemps);
	 Ytbeta[15]=YG8_xinter( Yi, YTemps);
	 Ytbeta[16]=YG1_r_float( Yi, YTemps-1);
	 Ytbeta[17]=YG2_r_float( Yi, YTemps-1);
	 memset(Yjac, 0, 36*sizeof(double));
	 Yur(Yi)->backward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4], Yting[5], Yting[6], Yting[7], Yting[8], Yting[9], Yting[10], Yting[11], Yting[12], Yting[13], Yting[14], Yting[15], Yting[16], Yting[17] );
	 Yvsmat (YNBS_ur, YNBI_ur, YMAX_JAC_NBI, &YG1Y_ur, Yjac[0], Ytbeta);
	 if (Ytestad_module)
	 {	memcpy(&YD1Y_ur, &YG1Y_ur, YNBS_ur*sizeof(YREAL));
	 }
	 Yting[0]=YS1_r_float( Yi, YTemps-1);
	 Yting[1]=YS2_r_float( Yi, YTemps-1);
	 Yting[2]=YS1_ur( Yi, YTemps);
	 Yting[3]=YS2_ur( Yi, YTemps);
	 Ytbeta[0]=YG1_r_float( Yi, YTemps-1);
	 Ytbeta[1]=YG2_r_float( Yi, YTemps-1);
	 Ytbeta[2]=YG1_ur( Yi, YTemps);
	 Ytbeta[3]=YG2_ur( Yi, YTemps);
	 memset(Yjac, 0, 36*sizeof(double));
	 Yr_float(Yi)->backward( Yting[0], Yting[1], Yting[2], Yting[3] );
	 Yvsmat (YNBS_r_float, YNBI_r_float, YMAX_JAC_NBI, &YG1Y_r_float, Yjac[0], Ytbeta);
	 if (Ytestad_module)
	 {	memcpy(&YD1Y_r_float, &YG1Y_r_float, YNBS_r_float*sizeof(YREAL));
	 }

}

	 return(0);
}

int Ydfward_space_S_float(int modop, char *nmmod, int All, int KeKo, int koleft, float pdx, float ptol, int yi, int yj, int yk)
{
	 int nbko=0;
	 YA1=1; YA2=0; YA3=0;
Yi=-1; Yj=-1; Yk=-1; /* init des indices de maille: maj par la boucle si valide*/
for(Yi=0; Yi<YA1_S_float; ++Yi)
{
	 Yting[0]=YS1_uinter( Yi, YTemps);
	 Yting[1]=YS2_uinter( Yi, YTemps);
	 Yting[2]=YS3_uinter( Yi, YTemps);
	 Yting[3]=YS4_uinter( Yi, YTemps);
	 Yting[4]=YS5_uinter( Yi, YTemps);
	 Yting[5]=YS6_uinter( Yi, YTemps);
	 Yting[6]=YS7_uinter( Yi, YTemps);
	 Yting[7]=YS8_uinter( Yi, YTemps);
	 Yting[8]=YS1_xinter( Yi, YTemps);
	 Yting[9]=YS2_xinter( Yi, YTemps);
	 Yting[10]=YS3_xinter( Yi, YTemps);
	 Yting[11]=YS4_xinter( Yi, YTemps);
	 Yting[12]=YS5_xinter( Yi, YTemps);
	 Yting[13]=YS6_xinter( Yi, YTemps);
	 Yting[14]=YS7_xinter( Yi, YTemps);
	 Yting[15]=YS8_xinter( Yi, YTemps);
	 Yting[16]=YS1_r_float( Yi, YTemps-1);
	 Yting[17]=YS2_r_float( Yi, YTemps-1);
	if (!strcmp(nmmod, "ur") || All)
	{
	 if (Yi>=yi && Yj>=yj && Yk>=yk) 
	 {	nbko += Ytestdf_ur(modop,KeKo,pdx,ptol,YNBI_ur,YNBS_ur,"ur",Yting,Yur[Yi]->Ystate[0],Yur[Yi]);
		 if (nbko>=koleft) return(nbko);
	 }
	}
	 Yur(Yi)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4], Yting[5], Yting[6], Yting[7], Yting[8], Yting[9], Yting[10], Yting[11], Yting[12], Yting[13], Yting[14], Yting[15], Yting[16], Yting[17] );
	 Yting[0]=YS1_r_float( Yi, YTemps-1);
	 Yting[1]=YS2_r_float( Yi, YTemps-1);
	 Yting[2]=YS1_ur( Yi, YTemps);
	 Yting[3]=YS2_ur( Yi, YTemps);
	if (!strcmp(nmmod, "r_float") || All)
	{
	 if (Yi>=yi && Yj>=yj && Yk>=yk) 
	 {	nbko += Ytestdf_r_float(modop,KeKo,pdx,ptol,YNBI_r_float,YNBS_r_float,"r_float",Yting,Yr_float[Yi]->Ystate[0],Yr_float[Yi]);
		 if (nbko>=koleft) return(nbko);
	 }
	}
	 Yr_float(Yi)->forward( Yting[0], Yting[1], Yting[2], Yting[3] );

}

	 return(nbko);
}

int Ydfward_space_maille_S_float(int modop, char *nmmod, int All, int KeKo,  float pdx, float ptol)
{
	 int nbko=0;
	 Yting[0]=YS1_uinter( Yi, YTemps);
	 Yting[1]=YS2_uinter( Yi, YTemps);
	 Yting[2]=YS3_uinter( Yi, YTemps);
	 Yting[3]=YS4_uinter( Yi, YTemps);
	 Yting[4]=YS5_uinter( Yi, YTemps);
	 Yting[5]=YS6_uinter( Yi, YTemps);
	 Yting[6]=YS7_uinter( Yi, YTemps);
	 Yting[7]=YS8_uinter( Yi, YTemps);
	 Yting[8]=YS1_xinter( Yi, YTemps);
	 Yting[9]=YS2_xinter( Yi, YTemps);
	 Yting[10]=YS3_xinter( Yi, YTemps);
	 Yting[11]=YS4_xinter( Yi, YTemps);
	 Yting[12]=YS5_xinter( Yi, YTemps);
	 Yting[13]=YS6_xinter( Yi, YTemps);
	 Yting[14]=YS7_xinter( Yi, YTemps);
	 Yting[15]=YS8_xinter( Yi, YTemps);
	 Yting[16]=YS1_r_float( Yi, YTemps-1);
	 Yting[17]=YS2_r_float( Yi, YTemps-1);
	if (!strcmp(nmmod, "ur") || All)
	{
	 nbko += Ytestdf_ur(modop,KeKo,pdx,ptol,YNBI_ur,YNBS_ur,"ur",Yting,Yur[Yi]->Ystate[0],Yur[Yi]);
	}
	 Yur(Yi)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4], Yting[5], Yting[6], Yting[7], Yting[8], Yting[9], Yting[10], Yting[11], Yting[12], Yting[13], Yting[14], Yting[15], Yting[16], Yting[17] );
	 Yting[0]=YS1_r_float( Yi, YTemps-1);
	 Yting[1]=YS2_r_float( Yi, YTemps-1);
	 Yting[2]=YS1_ur( Yi, YTemps);
	 Yting[3]=YS2_ur( Yi, YTemps);
	if (!strcmp(nmmod, "r_float") || All)
	{
	 nbko += Ytestdf_r_float(modop,KeKo,pdx,ptol,YNBI_r_float,YNBS_r_float,"r_float",Yting,Yr_float[Yi]->Ystate[0],Yr_float[Yi]);
	}
	 Yr_float(Yi)->forward( Yting[0], Yting[1], Yting[2], Yting[3] );
	 return(nbko);
}

 int Ybackward_space_S_float()
{
	 YA1=1; YA2=0; YA3=0;
Yi=-1; Yj=-1; Yk=-1; /* init des indices de maille: maj par la boucle si valide*/
for(Yi=YA1_S_float-1; Yi>=0; --Yi)
{
	 Yting[0]=YS1_r_float( Yi, YTemps-1);
	 Yting[1]=YS2_r_float( Yi, YTemps-1);
	 Yting[2]=YS1_ur( Yi, YTemps);
	 Yting[3]=YS2_ur( Yi, YTemps);
	 memset(Yjac, 0, 36*sizeof(double));
	 Yr_float(Yi)->backward( Yting[0], Yting[1], Yting[2], Yting[3] );
	 Yvsmatt (YNBS_r_float, YNBI_r_float, YMAX_JAC_NBI, &YG1Y_r_float, Yjac[0], Ytbeta);
	 YG1_r_float( Yi, YTemps-1) += Ytbeta[0];
	 YG2_r_float( Yi, YTemps-1) += Ytbeta[1];
	 YG1_ur( Yi, YTemps) += Ytbeta[2];
	 YG2_ur( Yi, YTemps) += Ytbeta[3];
	 if (Ytestad_module) {
			Yting[0]=YD1_r_float( Yi, YTemps-1);
			Yting[1]=YD2_r_float( Yi, YTemps-1);
			Yting[2]=YD1_ur( Yi, YTemps);
			Yting[3]=YD2_ur( Yi, YTemps);
			YLTRes = Yprosca(&YD1Y_r_float, &YG1Y_r_float, YNBS_r_float);
			YAdRes = Yprosca(Yting, Ytbeta, YNBI_r_float);
			if (!Ytesterrad_mod("r_float", YLTRes,YAdRes)) return(0);
	 }
	 Yting[0]=YS1_uinter( Yi, YTemps);
	 Yting[1]=YS2_uinter( Yi, YTemps);
	 Yting[2]=YS3_uinter( Yi, YTemps);
	 Yting[3]=YS4_uinter( Yi, YTemps);
	 Yting[4]=YS5_uinter( Yi, YTemps);
	 Yting[5]=YS6_uinter( Yi, YTemps);
	 Yting[6]=YS7_uinter( Yi, YTemps);
	 Yting[7]=YS8_uinter( Yi, YTemps);
	 Yting[8]=YS1_xinter( Yi, YTemps);
	 Yting[9]=YS2_xinter( Yi, YTemps);
	 Yting[10]=YS3_xinter( Yi, YTemps);
	 Yting[11]=YS4_xinter( Yi, YTemps);
	 Yting[12]=YS5_xinter( Yi, YTemps);
	 Yting[13]=YS6_xinter( Yi, YTemps);
	 Yting[14]=YS7_xinter( Yi, YTemps);
	 Yting[15]=YS8_xinter( Yi, YTemps);
	 Yting[16]=YS1_r_float( Yi, YTemps-1);
	 Yting[17]=YS2_r_float( Yi, YTemps-1);
	 memset(Yjac, 0, 36*sizeof(double));
	 Yur(Yi)->backward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4], Yting[5], Yting[6], Yting[7], Yting[8], Yting[9], Yting[10], Yting[11], Yting[12], Yting[13], Yting[14], Yting[15], Yting[16], Yting[17] );
	 Yvsmatt (YNBS_ur, YNBI_ur, YMAX_JAC_NBI, &YG1Y_ur, Yjac[0], Ytbeta);
	 YG1_uinter( Yi, YTemps) += Ytbeta[0];
	 YG2_uinter( Yi, YTemps) += Ytbeta[1];
	 YG3_uinter( Yi, YTemps) += Ytbeta[2];
	 YG4_uinter( Yi, YTemps) += Ytbeta[3];
	 YG5_uinter( Yi, YTemps) += Ytbeta[4];
	 YG6_uinter( Yi, YTemps) += Ytbeta[5];
	 YG7_uinter( Yi, YTemps) += Ytbeta[6];
	 YG8_uinter( Yi, YTemps) += Ytbeta[7];
	 YG1_xinter( Yi, YTemps) += Ytbeta[8];
	 YG2_xinter( Yi, YTemps) += Ytbeta[9];
	 YG3_xinter( Yi, YTemps) += Ytbeta[10];
	 YG4_xinter( Yi, YTemps) += Ytbeta[11];
	 YG5_xinter( Yi, YTemps) += Ytbeta[12];
	 YG6_xinter( Yi, YTemps) += Ytbeta[13];
	 YG7_xinter( Yi, YTemps) += Ytbeta[14];
	 YG8_xinter( Yi, YTemps) += Ytbeta[15];
	 YG1_r_float( Yi, YTemps-1) += Ytbeta[16];
	 YG2_r_float( Yi, YTemps-1) += Ytbeta[17];
	 if (Ytestad_module) {
			Yting[0]=YD1_uinter( Yi, YTemps);
			Yting[1]=YD2_uinter( Yi, YTemps);
			Yting[2]=YD3_uinter( Yi, YTemps);
			Yting[3]=YD4_uinter( Yi, YTemps);
			Yting[4]=YD5_uinter( Yi, YTemps);
			Yting[5]=YD6_uinter( Yi, YTemps);
			Yting[6]=YD7_uinter( Yi, YTemps);
			Yting[7]=YD8_uinter( Yi, YTemps);
			Yting[8]=YD1_xinter( Yi, YTemps);
			Yting[9]=YD2_xinter( Yi, YTemps);
			Yting[10]=YD3_xinter( Yi, YTemps);
			Yting[11]=YD4_xinter( Yi, YTemps);
			Yting[12]=YD5_xinter( Yi, YTemps);
			Yting[13]=YD6_xinter( Yi, YTemps);
			Yting[14]=YD7_xinter( Yi, YTemps);
			Yting[15]=YD8_xinter( Yi, YTemps);
			Yting[16]=YD1_r_float( Yi, YTemps-1);
			Yting[17]=YD2_r_float( Yi, YTemps-1);
			YLTRes = Yprosca(&YD1Y_ur, &YG1Y_ur, YNBS_ur);
			YAdRes = Yprosca(Yting, Ytbeta, YNBI_ur);
			if (!Ytesterrad_mod("ur", YLTRes,YAdRes)) return(0);
	 }

}

	 return(0);
}

 int Yforward_space_S_eulerlocate()
{
	 YA1=1; YA2=87; YA3=58;
Yi=-1; Yj=-1; Yk=-1; /* init des indices de maille: maj par la boucle si valide*/
for(Yi=0; Yi<YA1_S_eulerlocate; ++Yi)
for(Yj=0; Yj<YA2_S_eulerlocate; ++Yj)
for(Yk=0; Yk<YA3_S_eulerlocate; ++Yk)
{
	 Yting[0]=YS1_r_float( Yi, YTemps-1);
	 Yting[1]=YS2_r_float( Yi, YTemps-1);
	 Ylocate(Yi, Yj, Yk)->forward( Yting[0], Yting[1] );

}

	 return(0);
}

 int Ylinward_space_S_eulerlocate()
{
	 YA1=1; YA2=87; YA3=58;
Yi=-1; Yj=-1; Yk=-1; /* init des indices de maille: maj par la boucle si valide*/
for(Yi=0; Yi<YA1_S_eulerlocate; ++Yi)
for(Yj=0; Yj<YA2_S_eulerlocate; ++Yj)
for(Yk=0; Yk<YA3_S_eulerlocate; ++Yk)
{
	 Yting[0]=YS1_r_float( Yi, YTemps-1);
	 Yting[1]=YS2_r_float( Yi, YTemps-1);
	 Ytbeta[0]=YG1_r_float( Yi, YTemps-1);
	 Ytbeta[1]=YG2_r_float( Yi, YTemps-1);
	 memset(Yjac, 0, 18*sizeof(double));
	 Ylocate(Yi, Yj, Yk)->backward( Yting[0], Yting[1] );
	 Yvsmat (YNBS_locate, YNBI_locate, YMAX_JAC_NBI, &YG1Y_locate, Yjac[0], Ytbeta);
	 if (Ytestad_module)
	 {	memcpy(&YD1Y_locate, &YG1Y_locate, YNBS_locate*sizeof(YREAL));
	 }

}

	 return(0);
}

int Ydfward_space_S_eulerlocate(int modop, char *nmmod, int All, int KeKo, int koleft, float pdx, float ptol, int yi, int yj, int yk)
{
	 int nbko=0;
	 YA1=1; YA2=87; YA3=58;
Yi=-1; Yj=-1; Yk=-1; /* init des indices de maille: maj par la boucle si valide*/
for(Yi=0; Yi<YA1_S_eulerlocate; ++Yi)
for(Yj=0; Yj<YA2_S_eulerlocate; ++Yj)
for(Yk=0; Yk<YA3_S_eulerlocate; ++Yk)
{
	 Yting[0]=YS1_r_float( Yi, YTemps-1);
	 Yting[1]=YS2_r_float( Yi, YTemps-1);
	if (!strcmp(nmmod, "locate") || All)
	{
	 if (Yi>=yi && Yj>=yj && Yk>=yk) 
	 {	nbko += Ytestdf_locate(modop,KeKo,pdx,ptol,YNBI_locate,YNBS_locate,"locate",Yting,Ylocate[Yi][Yj][Yk]->Ystate[0],Ylocate[Yi][Yj][Yk]);
		 if (nbko>=koleft) return(nbko);
	 }
	}
	 Ylocate(Yi, Yj, Yk)->forward( Yting[0], Yting[1] );

}

	 return(nbko);
}

int Ydfward_space_maille_S_eulerlocate(int modop, char *nmmod, int All, int KeKo,  float pdx, float ptol)
{
	 int nbko=0;
	 Yting[0]=YS1_r_float( Yi, YTemps-1);
	 Yting[1]=YS2_r_float( Yi, YTemps-1);
	if (!strcmp(nmmod, "locate") || All)
	{
	 nbko += Ytestdf_locate(modop,KeKo,pdx,ptol,YNBI_locate,YNBS_locate,"locate",Yting,Ylocate[Yi][Yj][Yk]->Ystate[0],Ylocate[Yi][Yj][Yk]);
	}
	 Ylocate(Yi, Yj, Yk)->forward( Yting[0], Yting[1] );
	 return(nbko);
}

 int Ybackward_space_S_eulerlocate()
{
	 YA1=1; YA2=87; YA3=58;
Yi=-1; Yj=-1; Yk=-1; /* init des indices de maille: maj par la boucle si valide*/
for(Yi=YA1_S_eulerlocate-1; Yi>=0; --Yi)
for(Yj=YA2_S_eulerlocate-1; Yj>=0; --Yj)
for(Yk=YA3_S_eulerlocate-1; Yk>=0; --Yk)
{
	 Yting[0]=YS1_r_float( Yi, YTemps-1);
	 Yting[1]=YS2_r_float( Yi, YTemps-1);
	 memset(Yjac, 0, 18*sizeof(double));
	 Ylocate(Yi, Yj, Yk)->backward( Yting[0], Yting[1] );
	 Yvsmatt (YNBS_locate, YNBI_locate, YMAX_JAC_NBI, &YG1Y_locate, Yjac[0], Ytbeta);
	 YG1_r_float( Yi, YTemps-1) += Ytbeta[0];
	 YG2_r_float( Yi, YTemps-1) += Ytbeta[1];
	 if (Ytestad_module) {
			Yting[0]=YD1_r_float( Yi, YTemps-1);
			Yting[1]=YD2_r_float( Yi, YTemps-1);
			YLTRes = Yprosca(&YD1Y_locate, &YG1Y_locate, YNBS_locate);
			YAdRes = Yprosca(Yting, Ytbeta, YNBI_locate);
			if (!Ytesterrad_mod("locate", YLTRes,YAdRes)) return(0);
	 }

}

	 return(0);
}

int Yforward_traj_T_float(int nbp)
{
	Ycurward=FORWARD;
	 Yforward_space_S_eulerlocate();
	 Yforward_space_S_float();
	 return(0);
}

int Ylinward_traj_T_float(int nbp)
{
	Ycurward=LINWARD;
	 Ylinward_space_S_eulerlocate();
	 Ylinward_space_S_float();
	 return(0);
}

int Ydfward_traj_T_float(int modop, char *nmmod, int All, int KeKo, int koleft,float pdx, float ptol, int yi, int yj, int yk)
{
	 int nbko=0;
	 if (!Ydftesttt(0)) return(0);
	 nbko += Ydfward_space_S_eulerlocate(modop, nmmod, All, KeKo, koleft, pdx, ptol, yi, yj, yk);
	 nbko += Ydfward_space_S_float(modop, nmmod, All, KeKo, koleft, pdx, ptol, yi, yj, yk);
	 return(nbko);
}

int Ydfward_traj_maille_T_float(int modop, char *nmmod, int All, int KeKo,  float pdx, float ptol)
{
	 int nbko=0;
	 if (!Ydftesttt(0)) return(0);
	 nbko += Ydfward_space_maille_S_eulerlocate(modop, nmmod, All, KeKo, pdx, ptol);
	 nbko += Ydfward_space_maille_S_float(modop, nmmod, All, KeKo, pdx, ptol);
	 return(nbko);
}

int Ybackward_traj_T_float(int nbp)
{
	Ycurward=BACKWARD;
	 Ybackward_space_S_float();
	 Ybackward_space_S_eulerlocate();
	 return(0);
}

 int Yforward_space_S_float_d()
{
	 YA1=1; YA2=0; YA3=0;
Yi=-1; Yj=-1; Yk=-1; /* init des indices de maille: maj par la boucle si valide*/
for(Yi=0; Yi<YA1_S_float_d; ++Yi)
{
	 Yting[0]=YS1_uinter_d( Yi, YTemps);
	 Yting[1]=YS2_uinter_d( Yi, YTemps);
	 Yting[2]=YS3_uinter_d( Yi, YTemps);
	 Yting[3]=YS4_uinter_d( Yi, YTemps);
	 Yting[4]=YS5_uinter_d( Yi, YTemps);
	 Yting[5]=YS6_uinter_d( Yi, YTemps);
	 Yting[6]=YS7_uinter_d( Yi, YTemps);
	 Yting[7]=YS8_uinter_d( Yi, YTemps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[8]=YS1_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[9]=YS2_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[10]=YS3_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[11]=YS4_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[12]=YS5_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[13]=YS6_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[14]=YS7_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[15]=YS8_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps-1, YidTraj);
	 Yting[16]=YS1_r_float( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps-1, YidTraj);
	 Yting[17]=YS2_r_float( Yi, Ytps);
	 Yur_d(Yi)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4], Yting[5], Yting[6], Yting[7], Yting[8], Yting[9], Yting[10], Yting[11], Yting[12], Yting[13], Yting[14], Yting[15], Yting[16], Yting[17] );
	 Yting[0]=YS1_r_float_d( Yi, YTemps-1);
	 Yting[1]=YS2_r_float_d( Yi, YTemps-1);
	 Yting[2]=YS1_ur_d( Yi, YTemps);
	 Yting[3]=YS2_ur_d( Yi, YTemps);
	 Yr_float_d(Yi)->forward( Yting[0], Yting[1], Yting[2], Yting[3] );
	 Yting[0]=YS1_r_float_d( Yi, YTemps);
	 Yting[1]=YS2_r_float_d( Yi, YTemps);
	 Yr_cout_d(Yi)->forward( Yting[0], Yting[1] );

}

	 return(0);
}

 int Ylinward_space_S_float_d()
{
	 YA1=1; YA2=0; YA3=0;
Yi=-1; Yj=-1; Yk=-1; /* init des indices de maille: maj par la boucle si valide*/
for(Yi=0; Yi<YA1_S_float_d; ++Yi)
{
	 Yting[0]=YS1_uinter_d( Yi, YTemps);
	 Yting[1]=YS2_uinter_d( Yi, YTemps);
	 Yting[2]=YS3_uinter_d( Yi, YTemps);
	 Yting[3]=YS4_uinter_d( Yi, YTemps);
	 Yting[4]=YS5_uinter_d( Yi, YTemps);
	 Yting[5]=YS6_uinter_d( Yi, YTemps);
	 Yting[6]=YS7_uinter_d( Yi, YTemps);
	 Yting[7]=YS8_uinter_d( Yi, YTemps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[8]=YS1_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[9]=YS2_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[10]=YS3_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[11]=YS4_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[12]=YS5_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[13]=YS6_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[14]=YS7_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[15]=YS8_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps-1, YidTraj);
	 Yting[16]=YS1_r_float( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps-1, YidTraj);
	 Yting[17]=YS2_r_float( Yi, Ytps);
	 Ytbeta[0]=YG1_uinter_d( Yi, YTemps);
	 Ytbeta[1]=YG2_uinter_d( Yi, YTemps);
	 Ytbeta[2]=YG3_uinter_d( Yi, YTemps);
	 Ytbeta[3]=YG4_uinter_d( Yi, YTemps);
	 Ytbeta[4]=YG5_uinter_d( Yi, YTemps);
	 Ytbeta[5]=YG6_uinter_d( Yi, YTemps);
	 Ytbeta[6]=YG7_uinter_d( Yi, YTemps);
	 Ytbeta[7]=YG8_uinter_d( Yi, YTemps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Ytbeta[8]=YG1_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Ytbeta[9]=YG2_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Ytbeta[10]=YG3_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Ytbeta[11]=YG4_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Ytbeta[12]=YG5_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Ytbeta[13]=YG6_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Ytbeta[14]=YG7_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Ytbeta[15]=YG8_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps-1, YidTraj);
	 Ytbeta[16]=YG1_r_float( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps-1, YidTraj);
	 Ytbeta[17]=YG2_r_float( Yi, Ytps);
	 memset(Yjac, 0, 36*sizeof(double));
	 Yur_d(Yi)->backward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4], Yting[5], Yting[6], Yting[7], Yting[8], Yting[9], Yting[10], Yting[11], Yting[12], Yting[13], Yting[14], Yting[15], Yting[16], Yting[17] );
	 Yvsmat (YNBS_ur_d, YNBI_ur_d, YMAX_JAC_NBI, &YG1Y_ur_d, Yjac[0], Ytbeta);
	 if (Ytestad_module)
	 {	memcpy(&YD1Y_ur_d, &YG1Y_ur_d, YNBS_ur_d*sizeof(YREAL));
	 }
	 Yting[0]=YS1_r_float_d( Yi, YTemps-1);
	 Yting[1]=YS2_r_float_d( Yi, YTemps-1);
	 Yting[2]=YS1_ur_d( Yi, YTemps);
	 Yting[3]=YS2_ur_d( Yi, YTemps);
	 Ytbeta[0]=YG1_r_float_d( Yi, YTemps-1);
	 Ytbeta[1]=YG2_r_float_d( Yi, YTemps-1);
	 Ytbeta[2]=YG1_ur_d( Yi, YTemps);
	 Ytbeta[3]=YG2_ur_d( Yi, YTemps);
	 memset(Yjac, 0, 36*sizeof(double));
	 Yr_float_d(Yi)->backward( Yting[0], Yting[1], Yting[2], Yting[3] );
	 Yvsmat (YNBS_r_float_d, YNBI_r_float_d, YMAX_JAC_NBI, &YG1Y_r_float_d, Yjac[0], Ytbeta);
	 if (Ytestad_module)
	 {	memcpy(&YD1Y_r_float_d, &YG1Y_r_float_d, YNBS_r_float_d*sizeof(YREAL));
	 }
	 Yting[0]=YS1_r_float_d( Yi, YTemps);
	 Yting[1]=YS2_r_float_d( Yi, YTemps);
	 Ytbeta[0]=YG1_r_float_d( Yi, YTemps);
	 Ytbeta[1]=YG2_r_float_d( Yi, YTemps);
	 memset(Yjac, 0, 36*sizeof(double));
	 Yr_cout_d(Yi)->backward( Yting[0], Yting[1] );
	 Yvsmat (YNBS_r_cout_d, YNBI_r_cout_d, YMAX_JAC_NBI, &YG1Y_r_cout_d, Yjac[0], Ytbeta);
	 if (Ytestad_module)
	 {	memcpy(&YD1Y_r_cout_d, &YG1Y_r_cout_d, YNBS_r_cout_d*sizeof(YREAL));
	 }

}

	 return(0);
}

int Ydfward_space_S_float_d(int modop, char *nmmod, int All, int KeKo, int koleft, float pdx, float ptol, int yi, int yj, int yk)
{
	 int nbko=0;
	 YA1=1; YA2=0; YA3=0;
Yi=-1; Yj=-1; Yk=-1; /* init des indices de maille: maj par la boucle si valide*/
for(Yi=0; Yi<YA1_S_float_d; ++Yi)
{
	 Yting[0]=YS1_uinter_d( Yi, YTemps);
	 Yting[1]=YS2_uinter_d( Yi, YTemps);
	 Yting[2]=YS3_uinter_d( Yi, YTemps);
	 Yting[3]=YS4_uinter_d( Yi, YTemps);
	 Yting[4]=YS5_uinter_d( Yi, YTemps);
	 Yting[5]=YS6_uinter_d( Yi, YTemps);
	 Yting[6]=YS7_uinter_d( Yi, YTemps);
	 Yting[7]=YS8_uinter_d( Yi, YTemps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[8]=YS1_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[9]=YS2_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[10]=YS3_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[11]=YS4_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[12]=YS5_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[13]=YS6_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[14]=YS7_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[15]=YS8_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps-1, YidTraj);
	 Yting[16]=YS1_r_float( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps-1, YidTraj);
	 Yting[17]=YS2_r_float( Yi, Ytps);
	if (!strcmp(nmmod, "ur_d") || All)
	{
	 if (Yi>=yi && Yj>=yj && Yk>=yk) 
	 {	nbko += Ytestdf_ur_d(modop,KeKo,pdx,ptol,YNBI_ur_d,YNBS_ur_d,"ur_d",Yting,Yur_d[Yi]->Ystate[0],Yur_d[Yi]);
		 if (nbko>=koleft) return(nbko);
	 }
	}
	 Yur_d(Yi)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4], Yting[5], Yting[6], Yting[7], Yting[8], Yting[9], Yting[10], Yting[11], Yting[12], Yting[13], Yting[14], Yting[15], Yting[16], Yting[17] );
	 Yting[0]=YS1_r_float_d( Yi, YTemps-1);
	 Yting[1]=YS2_r_float_d( Yi, YTemps-1);
	 Yting[2]=YS1_ur_d( Yi, YTemps);
	 Yting[3]=YS2_ur_d( Yi, YTemps);
	if (!strcmp(nmmod, "r_float_d") || All)
	{
	 if (Yi>=yi && Yj>=yj && Yk>=yk) 
	 {	nbko += Ytestdf_r_float_d(modop,KeKo,pdx,ptol,YNBI_r_float_d,YNBS_r_float_d,"r_float_d",Yting,Yr_float_d[Yi]->Ystate[0],Yr_float_d[Yi]);
		 if (nbko>=koleft) return(nbko);
	 }
	}
	 Yr_float_d(Yi)->forward( Yting[0], Yting[1], Yting[2], Yting[3] );
	 Yting[0]=YS1_r_float_d( Yi, YTemps);
	 Yting[1]=YS2_r_float_d( Yi, YTemps);
	if (!strcmp(nmmod, "r_cout_d") || All)
	{
	 if (Yi>=yi && Yj>=yj && Yk>=yk) 
	 {	nbko += Ytestdf_r_cout_d(modop,KeKo,pdx,ptol,YNBI_r_cout_d,YNBS_r_cout_d,"r_cout_d",Yting,Yr_cout_d[Yi]->Ystate[0],Yr_cout_d[Yi]);
		 if (nbko>=koleft) return(nbko);
	 }
	}
	 Yr_cout_d(Yi)->forward( Yting[0], Yting[1] );

}

	 return(nbko);
}

int Ydfward_space_maille_S_float_d(int modop, char *nmmod, int All, int KeKo,  float pdx, float ptol)
{
	 int nbko=0;
	 Yting[0]=YS1_uinter_d( Yi, YTemps);
	 Yting[1]=YS2_uinter_d( Yi, YTemps);
	 Yting[2]=YS3_uinter_d( Yi, YTemps);
	 Yting[3]=YS4_uinter_d( Yi, YTemps);
	 Yting[4]=YS5_uinter_d( Yi, YTemps);
	 Yting[5]=YS6_uinter_d( Yi, YTemps);
	 Yting[6]=YS7_uinter_d( Yi, YTemps);
	 Yting[7]=YS8_uinter_d( Yi, YTemps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[8]=YS1_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[9]=YS2_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[10]=YS3_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[11]=YS4_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[12]=YS5_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[13]=YS6_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[14]=YS7_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[15]=YS8_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps-1, YidTraj);
	 Yting[16]=YS1_r_float( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps-1, YidTraj);
	 Yting[17]=YS2_r_float( Yi, Ytps);
	if (!strcmp(nmmod, "ur_d") || All)
	{
	 nbko += Ytestdf_ur_d(modop,KeKo,pdx,ptol,YNBI_ur_d,YNBS_ur_d,"ur_d",Yting,Yur_d[Yi]->Ystate[0],Yur_d[Yi]);
	}
	 Yur_d(Yi)->forward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4], Yting[5], Yting[6], Yting[7], Yting[8], Yting[9], Yting[10], Yting[11], Yting[12], Yting[13], Yting[14], Yting[15], Yting[16], Yting[17] );
	 Yting[0]=YS1_r_float_d( Yi, YTemps-1);
	 Yting[1]=YS2_r_float_d( Yi, YTemps-1);
	 Yting[2]=YS1_ur_d( Yi, YTemps);
	 Yting[3]=YS2_ur_d( Yi, YTemps);
	if (!strcmp(nmmod, "r_float_d") || All)
	{
	 nbko += Ytestdf_r_float_d(modop,KeKo,pdx,ptol,YNBI_r_float_d,YNBS_r_float_d,"r_float_d",Yting,Yr_float_d[Yi]->Ystate[0],Yr_float_d[Yi]);
	}
	 Yr_float_d(Yi)->forward( Yting[0], Yting[1], Yting[2], Yting[3] );
	 Yting[0]=YS1_r_float_d( Yi, YTemps);
	 Yting[1]=YS2_r_float_d( Yi, YTemps);
	if (!strcmp(nmmod, "r_cout_d") || All)
	{
	 nbko += Ytestdf_r_cout_d(modop,KeKo,pdx,ptol,YNBI_r_cout_d,YNBS_r_cout_d,"r_cout_d",Yting,Yr_cout_d[Yi]->Ystate[0],Yr_cout_d[Yi]);
	}
	 Yr_cout_d(Yi)->forward( Yting[0], Yting[1] );
	 return(nbko);
}

 int Ybackward_space_S_float_d()
{
	 YA1=1; YA2=0; YA3=0;
Yi=-1; Yj=-1; Yk=-1; /* init des indices de maille: maj par la boucle si valide*/
for(Yi=YA1_S_float_d-1; Yi>=0; --Yi)
{
	 Yting[0]=YS1_r_float_d( Yi, YTemps);
	 Yting[1]=YS2_r_float_d( Yi, YTemps);
	 memset(Yjac, 0, 36*sizeof(double));
	 Yr_cout_d(Yi)->backward( Yting[0], Yting[1] );
	 Yvsmatt (YNBS_r_cout_d, YNBI_r_cout_d, YMAX_JAC_NBI, &YG1Y_r_cout_d, Yjac[0], Ytbeta);
	 YG1_r_float_d( Yi, YTemps) += Ytbeta[0];
	 YG2_r_float_d( Yi, YTemps) += Ytbeta[1];
	 if (Ytestad_module) {
			Yting[0]=YD1_r_float_d( Yi, YTemps);
			Yting[1]=YD2_r_float_d( Yi, YTemps);
			YLTRes = Yprosca(&YD1Y_r_cout_d, &YG1Y_r_cout_d, YNBS_r_cout_d);
			YAdRes = Yprosca(Yting, Ytbeta, YNBI_r_cout_d);
			if (!Ytesterrad_mod("r_cout_d", YLTRes,YAdRes)) return(0);
	 }
	 Yting[0]=YS1_r_float_d( Yi, YTemps-1);
	 Yting[1]=YS2_r_float_d( Yi, YTemps-1);
	 Yting[2]=YS1_ur_d( Yi, YTemps);
	 Yting[3]=YS2_ur_d( Yi, YTemps);
	 memset(Yjac, 0, 36*sizeof(double));
	 Yr_float_d(Yi)->backward( Yting[0], Yting[1], Yting[2], Yting[3] );
	 Yvsmatt (YNBS_r_float_d, YNBI_r_float_d, YMAX_JAC_NBI, &YG1Y_r_float_d, Yjac[0], Ytbeta);
	 YG1_r_float_d( Yi, YTemps-1) += Ytbeta[0];
	 YG2_r_float_d( Yi, YTemps-1) += Ytbeta[1];
	 YG1_ur_d( Yi, YTemps) += Ytbeta[2];
	 YG2_ur_d( Yi, YTemps) += Ytbeta[3];
	 if (Ytestad_module) {
			Yting[0]=YD1_r_float_d( Yi, YTemps-1);
			Yting[1]=YD2_r_float_d( Yi, YTemps-1);
			Yting[2]=YD1_ur_d( Yi, YTemps);
			Yting[3]=YD2_ur_d( Yi, YTemps);
			YLTRes = Yprosca(&YD1Y_r_float_d, &YG1Y_r_float_d, YNBS_r_float_d);
			YAdRes = Yprosca(Yting, Ytbeta, YNBI_r_float_d);
			if (!Ytesterrad_mod("r_float_d", YLTRes,YAdRes)) return(0);
	 }
	 Yting[0]=YS1_uinter_d( Yi, YTemps);
	 Yting[1]=YS2_uinter_d( Yi, YTemps);
	 Yting[2]=YS3_uinter_d( Yi, YTemps);
	 Yting[3]=YS4_uinter_d( Yi, YTemps);
	 Yting[4]=YS5_uinter_d( Yi, YTemps);
	 Yting[5]=YS6_uinter_d( Yi, YTemps);
	 Yting[6]=YS7_uinter_d( Yi, YTemps);
	 Yting[7]=YS8_uinter_d( Yi, YTemps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[8]=YS1_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[9]=YS2_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[10]=YS3_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[11]=YS4_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[12]=YS5_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[13]=YS6_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[14]=YS7_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 Yting[15]=YS8_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps-1, YidTraj);
	 Yting[16]=YS1_r_float( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps-1, YidTraj);
	 Yting[17]=YS2_r_float( Yi, Ytps);
	 memset(Yjac, 0, 36*sizeof(double));
	 Yur_d(Yi)->backward( Yting[0], Yting[1], Yting[2], Yting[3], Yting[4], Yting[5], Yting[6], Yting[7], Yting[8], Yting[9], Yting[10], Yting[11], Yting[12], Yting[13], Yting[14], Yting[15], Yting[16], Yting[17] );
	 Yvsmatt (YNBS_ur_d, YNBI_ur_d, YMAX_JAC_NBI, &YG1Y_ur_d, Yjac[0], Ytbeta);
	 YG1_uinter_d( Yi, YTemps) += Ytbeta[0];
	 YG2_uinter_d( Yi, YTemps) += Ytbeta[1];
	 YG3_uinter_d( Yi, YTemps) += Ytbeta[2];
	 YG4_uinter_d( Yi, YTemps) += Ytbeta[3];
	 YG5_uinter_d( Yi, YTemps) += Ytbeta[4];
	 YG6_uinter_d( Yi, YTemps) += Ytbeta[5];
	 YG7_uinter_d( Yi, YTemps) += Ytbeta[6];
	 YG8_uinter_d( Yi, YTemps) += Ytbeta[7];
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 YG1_xinter( Yi, Ytps) += Ytbeta[8];
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 YG2_xinter( Yi, Ytps) += Ytbeta[9];
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 YG3_xinter( Yi, Ytps) += Ytbeta[10];
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 YG4_xinter( Yi, Ytps) += Ytbeta[11];
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 YG5_xinter( Yi, Ytps) += Ytbeta[12];
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 YG6_xinter( Yi, Ytps) += Ytbeta[13];
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 YG7_xinter( Yi, Ytps) += Ytbeta[14];
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
	 YG8_xinter( Yi, Ytps) += Ytbeta[15];
Ytps = Ytttt_pdt(0, YTemps-1, YidTraj);
	 YG1_r_float( Yi, Ytps) += Ytbeta[16];
Ytps = Ytttt_pdt(0, YTemps-1, YidTraj);
	 YG2_r_float( Yi, Ytps) += Ytbeta[17];
	 if (Ytestad_module) {
			Yting[0]=YD1_uinter_d( Yi, YTemps);
			Yting[1]=YD2_uinter_d( Yi, YTemps);
			Yting[2]=YD3_uinter_d( Yi, YTemps);
			Yting[3]=YD4_uinter_d( Yi, YTemps);
			Yting[4]=YD5_uinter_d( Yi, YTemps);
			Yting[5]=YD6_uinter_d( Yi, YTemps);
			Yting[6]=YD7_uinter_d( Yi, YTemps);
			Yting[7]=YD8_uinter_d( Yi, YTemps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
			Yting[8]=YD1_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
			Yting[9]=YD2_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
			Yting[10]=YD3_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
			Yting[11]=YD4_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
			Yting[12]=YD5_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
			Yting[13]=YD6_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
			Yting[14]=YD7_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps, YidTraj);
			Yting[15]=YD8_xinter( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps-1, YidTraj);
			Yting[16]=YD1_r_float( Yi, Ytps);
Ytps = Ytttt_pdt(0, YTemps-1, YidTraj);
			Yting[17]=YD2_r_float( Yi, Ytps);
			YLTRes = Yprosca(&YD1Y_ur_d, &YG1Y_ur_d, YNBS_ur_d);
			YAdRes = Yprosca(Yting, Ytbeta, YNBI_ur_d);
			if (!Ytesterrad_mod("ur_d", YLTRes,YAdRes)) return(0);
	 }

}

	 return(0);
}

 int Yforward_space_S_eulerlocate_d()
{
	 YA1=1; YA2=87; YA3=58;
Yi=-1; Yj=-1; Yk=-1; /* init des indices de maille: maj par la boucle si valide*/
for(Yi=0; Yi<YA1_S_eulerlocate_d; ++Yi)
for(Yj=0; Yj<YA2_S_eulerlocate_d; ++Yj)
for(Yk=0; Yk<YA3_S_eulerlocate_d; ++Yk)
{
	 Yting[0]=YS1_r_float_d( Yi, YTemps-1);
	 Yting[1]=YS2_r_float_d( Yi, YTemps-1);
	 Ylocate_d(Yi, Yj, Yk)->forward( Yting[0], Yting[1] );

}

	 return(0);
}

 int Ylinward_space_S_eulerlocate_d()
{
	 YA1=1; YA2=87; YA3=58;
Yi=-1; Yj=-1; Yk=-1; /* init des indices de maille: maj par la boucle si valide*/
for(Yi=0; Yi<YA1_S_eulerlocate_d; ++Yi)
for(Yj=0; Yj<YA2_S_eulerlocate_d; ++Yj)
for(Yk=0; Yk<YA3_S_eulerlocate_d; ++Yk)
{
	 Yting[0]=YS1_r_float_d( Yi, YTemps-1);
	 Yting[1]=YS2_r_float_d( Yi, YTemps-1);
	 Ytbeta[0]=YG1_r_float_d( Yi, YTemps-1);
	 Ytbeta[1]=YG2_r_float_d( Yi, YTemps-1);
	 memset(Yjac, 0, 18*sizeof(double));
	 Ylocate_d(Yi, Yj, Yk)->backward( Yting[0], Yting[1] );
	 Yvsmat (YNBS_locate_d, YNBI_locate_d, YMAX_JAC_NBI, &YG1Y_locate_d, Yjac[0], Ytbeta);
	 if (Ytestad_module)
	 {	memcpy(&YD1Y_locate_d, &YG1Y_locate_d, YNBS_locate_d*sizeof(YREAL));
	 }

}

	 return(0);
}

int Ydfward_space_S_eulerlocate_d(int modop, char *nmmod, int All, int KeKo, int koleft, float pdx, float ptol, int yi, int yj, int yk)
{
	 int nbko=0;
	 YA1=1; YA2=87; YA3=58;
Yi=-1; Yj=-1; Yk=-1; /* init des indices de maille: maj par la boucle si valide*/
for(Yi=0; Yi<YA1_S_eulerlocate_d; ++Yi)
for(Yj=0; Yj<YA2_S_eulerlocate_d; ++Yj)
for(Yk=0; Yk<YA3_S_eulerlocate_d; ++Yk)
{
	 Yting[0]=YS1_r_float_d( Yi, YTemps-1);
	 Yting[1]=YS2_r_float_d( Yi, YTemps-1);
	if (!strcmp(nmmod, "locate_d") || All)
	{
	 if (Yi>=yi && Yj>=yj && Yk>=yk) 
	 {	nbko += Ytestdf_locate_d(modop,KeKo,pdx,ptol,YNBI_locate_d,YNBS_locate_d,"locate_d",Yting,Ylocate_d[Yi][Yj][Yk]->Ystate[0],Ylocate_d[Yi][Yj][Yk]);
		 if (nbko>=koleft) return(nbko);
	 }
	}
	 Ylocate_d(Yi, Yj, Yk)->forward( Yting[0], Yting[1] );

}

	 return(nbko);
}

int Ydfward_space_maille_S_eulerlocate_d(int modop, char *nmmod, int All, int KeKo,  float pdx, float ptol)
{
	 int nbko=0;
	 Yting[0]=YS1_r_float_d( Yi, YTemps-1);
	 Yting[1]=YS2_r_float_d( Yi, YTemps-1);
	if (!strcmp(nmmod, "locate_d") || All)
	{
	 nbko += Ytestdf_locate_d(modop,KeKo,pdx,ptol,YNBI_locate_d,YNBS_locate_d,"locate_d",Yting,Ylocate_d[Yi][Yj][Yk]->Ystate[0],Ylocate_d[Yi][Yj][Yk]);
	}
	 Ylocate_d(Yi, Yj, Yk)->forward( Yting[0], Yting[1] );
	 return(nbko);
}

 int Ybackward_space_S_eulerlocate_d()
{
	 YA1=1; YA2=87; YA3=58;
Yi=-1; Yj=-1; Yk=-1; /* init des indices de maille: maj par la boucle si valide*/
for(Yi=YA1_S_eulerlocate_d-1; Yi>=0; --Yi)
for(Yj=YA2_S_eulerlocate_d-1; Yj>=0; --Yj)
for(Yk=YA3_S_eulerlocate_d-1; Yk>=0; --Yk)
{
	 Yting[0]=YS1_r_float_d( Yi, YTemps-1);
	 Yting[1]=YS2_r_float_d( Yi, YTemps-1);
	 memset(Yjac, 0, 18*sizeof(double));
	 Ylocate_d(Yi, Yj, Yk)->backward( Yting[0], Yting[1] );
	 Yvsmatt (YNBS_locate_d, YNBI_locate_d, YMAX_JAC_NBI, &YG1Y_locate_d, Yjac[0], Ytbeta);
	 YG1_r_float_d( Yi, YTemps-1) += Ytbeta[0];
	 YG2_r_float_d( Yi, YTemps-1) += Ytbeta[1];
	 if (Ytestad_module) {
			Yting[0]=YD1_r_float_d( Yi, YTemps-1);
			Yting[1]=YD2_r_float_d( Yi, YTemps-1);
			YLTRes = Yprosca(&YD1Y_locate_d, &YG1Y_locate_d, YNBS_locate_d);
			YAdRes = Yprosca(Yting, Ytbeta, YNBI_locate_d);
			if (!Ytesterrad_mod("locate_d", YLTRes,YAdRes)) return(0);
	 }

}

	 return(0);
}

int Yforward_traj_T_float_d(int nbp)
{
	Ycurward=FORWARD;
	 Yforward_space_S_eulerlocate_d();
	 Yforward_space_S_float_d();
	 return(0);
}

int Ylinward_traj_T_float_d(int nbp)
{
	Ycurward=LINWARD;
	 Ylinward_space_S_eulerlocate_d();
	 Ylinward_space_S_float_d();
	 return(0);
}

int Ydfward_traj_T_float_d(int modop, char *nmmod, int All, int KeKo, int koleft,float pdx, float ptol, int yi, int yj, int yk)
{
	 int nbko=0;
	 if (!Ydftesttt(2)) return(0);
	 nbko += Ydfward_space_S_eulerlocate_d(modop, nmmod, All, KeKo, koleft, pdx, ptol, yi, yj, yk);
	 nbko += Ydfward_space_S_float_d(modop, nmmod, All, KeKo, koleft, pdx, ptol, yi, yj, yk);
	 return(nbko);
}

int Ydfward_traj_maille_T_float_d(int modop, char *nmmod, int All, int KeKo,  float pdx, float ptol)
{
	 int nbko=0;
	 if (!Ydftesttt(2)) return(0);
	 nbko += Ydfward_space_maille_S_eulerlocate_d(modop, nmmod, All, KeKo, pdx, ptol);
	 nbko += Ydfward_space_maille_S_float_d(modop, nmmod, All, KeKo, pdx, ptol);
	 return(nbko);
}

int Ybackward_traj_T_float_d(int nbp)
{
	Ycurward=BACKWARD;
	 Ybackward_space_S_float_d();
	 Ybackward_space_S_eulerlocate_d();
	 return(0);
}

int Yforward_order()
{
	 Yforward_traj_T_float(0);
	 Yforward_traj_T_float_d(0);
	 return(0);
}

int Ylinward_order()
{
	 Ylinward_traj_T_float(0);
	 Ylinward_traj_T_float_d(0);
	 return(0);
}

int Ydfward_order(int modop, char *nmmod, int All, int KeKo, int koleft, float pdx, float ptol, int yi, int yj, int yk)
{
	 int nbko=0;
	 nbko += Ydfward_traj_T_float(modop, nmmod, All, KeKo, koleft, pdx, ptol, yi, yj, yk);
	 nbko += Ydfward_traj_T_float_d(modop, nmmod, All, KeKo, koleft, pdx, ptol, yi, yj, yk);
	 return(nbko);
}

int Ydfward_order_maille(int modop, char *nmmod, int All, int KeKo,  float pdx, float ptol)
{
	 int nbko=0;
	 nbko += Ydfward_traj_maille_T_float(modop, nmmod, All, KeKo, pdx, ptol);
	 nbko += Ydfward_traj_maille_T_float_d(modop, nmmod, All, KeKo, pdx, ptol);
	 return(nbko);
}

int Ybackward_order()
{
	 Ybackward_traj_T_float_d(0);
	 Ybackward_traj_T_float(0);
	 return(0);
}

void Yforward_operator (char type) 
{	
}
void Ylinward_operator (char type) 
{	
}
void Ybackward_operator (char type) 
{	
}

// ¤ ¤ ¤ ¤ ¤ ¤ ¤ ¤ LES FONCTIONS UTILISATEUR ... : 
int Yuser_call (int argc, char *argv[]) 
{
	 int codret=1;
	 if (strcmp(argv[0], "update_uv") == 0)
		 update_uv();
	 else if (strcmp(argv[0], "run_inc") == 0)
		 run_inc(argc, argv);
	 else codret=0;
	 return(codret);
}

//-------------- end Yao generation -----------------