Changeset 98 for altifloat/src/floater.h

Timestamp:

06/05/13 16:00:10 (11 years ago)

Author:

jbrlod

Message:

add function run_inc

File:

• altifloat/src/floater.h (modified) (3 diffs)

altifloat/src/floater.h

@@ -7 +7 @@
 YREAL vr_dy=0;
 
+
+
+//Options du run incrŽmental
+int nb_extiter=6;
+short inc_save=1;
+char dirsave[50]="../obs_float/";
+char savename[50]="optim";
+extern float Y3ddf1;
 //Declaration 
 
@@ -15 +23 @@
 void after_it(int nit){
         
-        printf("u_d(9,16)=%f %f %f %f\n",YS1_u_d(9,16),YS1_u_d(10,16),YS1_u_d(10,17),YS1_u_d(9,17));
+        //printf("u_d(9,16)=%f %f %f %f\n",YS1_u_d(9,16),YS1_u_d(10,16),YS1_u_d(10,17),YS1_u_d(9,17));
         }
 void forward_before(int ctrp){
@@ -50 +58 @@
         
 //      }
+
+void run_inc(int argv, char *argc[]) {
+//Run the incremental optimisation
+//1 forward of the complete model
+//2 initialize du et dv to zero
+//3 minimize the incremental cost function
+//4 update u and v
+
+/*
+Options
+nb_extiter : number of extern loop (default=6)
+inc_save : save option
+        0  : no save
+        1 (default): save after each extern loop r_float,u and v
+dirsave : directory in case inc_save>0 (default="../obs_float/")
+savename : basename of save fils (default "optim")
+*/
+char filesave[100];
+int i;
+
+char sactiv[20]="activ";
+char sMD[20]="M";
+char sonly[20]="only";
+char *liste[3];
+liste[0]=sactiv;
+liste[1]=sMD;
+liste[2]=sonly;
+/*Check M1QN3 config*/
+ if (YNbItRun<=0)
+                {        printf("runm(2): number of run iteration not seted; use set_nbiter please\n");
+                         //return(0);
+                }
+        if (Y3ddf1<=0)
+                {       printf("runm(2): expected positive fcost decrease missed; use setm_ddf1\n");
+                         //return(0);
+                }
+        if (YioInsertObsCtr<0)
+                   printf("runm(2): warning : oh oh, run with no obs !!! \n");
+                
+        YTypeAdjust = ADJUST_M1QN3; //d'office avec M1QN3
+        
+//Savestate config
+ YioModulot = OFF; //On sauvegarde tous les pas de temps
+    YioWrite = ON;
+   // Yio_savestate(cdes[1], cdes[3], 0, cdes[7]);
+    YioState=0 ;//Save all states
+    YioBin=OFF; //ascii output
+    YioAscii=ON; //ascii output
+    YioAxes=ON; //save axe numbers
+
+for (i=0;i<nb_extiter;i++) {
+     strcpy(liste[1],"M");
+         Yactraj(3, liste);
+         Yset_modeltime(0);
+         before_it(1);
+         Yforward(-1, 0);
+         
+         if (inc_save==1 && i>0) {
+          //save rfloat
+          YioTime=ON;
+          sprintf(filesave,"%sr_float_%s_%2d.dat",dirsave,savename,i+1);
+          Yio_savestate("r_float","i",0,filesave);
+         }
+             strcpy(liste[1],"D");
+         Yactraj(3,liste);
+         Yset_modeltime(0);
+         
+        
+        Y3run ('0');
+    update_uv();
+    if (inc_save==1) {
+   
+    
+    //save u
+    YioTime=OFF;
+    sprintf(filesave,"%su%s%2d.dat",dirsave,savename,i+1);
+    Yio_savestate("u","ij",0,filesave);
+    
+    //save v
+    sprintf(filesave,"%sv%s%2d.dat",dirsave,savename,i+1);
+    Yio_savestate("v","ij",0,filesave);
+    
+    }
+    }
+    /*End of the optimization, we run another time the model to compute the final value of r_float
+    */
+     strcpy(liste[1],"M");
+    Yactraj(3, liste);
+         Yset_modeltime(0);
+         before_it(1);
+         Yforward(-1, 0);
+    
+}