-- FreeFem++ v4.6 (Thu Apr  2 15:47:38 CEST 2020 - git v4.6)
 Load: lg_fem lg_mesh lg_mesh3 eigenvalue 
    1 : //  molding.edp
    2 : //
    3 : //  Discussion:
    4 : //
    5 : //    Simulation of injection molding.
    6 : //
    7 : //  Location:
    8 : //
    9 : //    http://people.sc.fsu.edu/~jburkardt/freefem++/molding/molding.edp
   10 : //
   11 : //  Modified:
   12 : //
   13 : //    20 May 2015
   14 : //
   15 : //  Author:
   16 : //
   17 : //    Florian De Vuyst
   18 : //
   19 : //  Reference:
   20 : //
   21 : //    Florian De Vuyst,
   22 : //    Numerical modeling of transport problems using freefem++ software -
   23 : //    with examples in biology, CFD, traffic flow and energy transfer,
   24 : //    HAL id: cel-00842234
   25 : //    https://cel.archives-ouvertes.fr/cel-00842234
   26 : //
   27 : real[int] PA(2); PA = [0, 0];
   28 : real[int] PB(2); PB = [9, 0];
   29 : real[int] PC(2); PC = [9, 9];
   30 : real[int] PD(2); PD = [8, 9];
   31 : real[int] PE(2); PE = [0, 9];
   32 : real [int ] PF(2); PF = [0, 6];
   33 : real [int ] PG(2); PG = [0,4];
   34 : real [int ] PH(2); PH = [1, 1];
   35 : real [int ] PK(2); PK = [1,8];
   36 : real [int ] PL(2); PL = [8,8];
   37 : real [int ] PM(2); PM = [8,5];
   38 : real [int ] PO(2); PO = [5,5];
   39 : real [int ] PP(2); PP = [5,4];
   40 : real [int ] PQ(2); PQ = [8,4];
   41 : real [int ] PR(2); PR = [8,1];
   42 : 
   43 : border c1(t=0,1) {x =(1-t)*PA[0]+t*PB[0]; y =(1-t)*PA[1]+t*PB[1];}
   44 : border c2(t=0,1) {x =(1-t)*PB[0]+t*PC[0]; y =(1-t)*PB[1]+t*PC[1];}
   45 : border c3(t=0,1) {x =(1-t)*PC[0]+t*PD[0]; y =(1-t)*PC[1]+t*PD[1];}
   46 : border c4(t=0,1) {x =(1-t)*PD[0]+t*PE[0]; y =(1-t)*PD[1]+t*PE[1];}
   47 : border c5(t=0,1) {x =(1-t)*PE[0]+t*PF[0]; y =(1-t)*PE[1]+t*PF[1];}
   48 : border c6(t=0,1) {x =(1-t)*PF[0]+t*PG[0]; y =(1-t)*PF[1]+t*PG[1];}
   49 : border c7(t=0,1) {x =(1-t)*PG[0]+t*PA[0]; y =(1-t)*PG[1]+t*PA[1];}
   50 : //
   51 : border c8(t=0,1){x=(1-t)*PH[0]+t*PK[0]; y=(1-t)*PH[1]+t*PK[1];}
   52 : border c9(t=0,1){x=(1-t)*PK[0]+t*PL[0]; y=(1-t)*PK[1]+t*PL[1];}
   53 : border c10(t=0,1){x=(1-t)*PL[0]+t*PM[0]; y=(1-t)*PL[1]+t*PM[1];}
   54 : border c11(t=0,1){x=(1-t)*PM[0]+t*PO[0]; y=(1-t)*PM[1]+t*PO[1];}
   55 : border c12(t=0,1){x=(1-t)*PO[0]+t*PP[0]; y=(1-t)*PO[1]+t*PP[1];}
   56 : border c13(t=0,1){x=(1-t)*PP[0]+t*PQ[0]; y=(1-t)*PP[1]+t*PQ[1];}
   57 : border c14(t=0,1){x=(1-t)*PQ[0]+t*PR[0]; y=(1-t)*PQ[1]+t*PR[1];}
   58 : border c15(t=0,1){x=(1-t)*PR[0]+t*PH[0]; y=(1-t)*PR[1]+t*PH[1];}
   59 : //
   60 : //  Define the mesh, Th.
   61 : //
   62 : int m = 2;
   63 : 
   64 : mesh Th = buildmesh ( c1(80*m)+c2(50*m)+c3(10*m)+c4(40*m)
   65 :   + c5(30*m)+c6(10*m)+c7(30*m)
   66 :   + c8(80*m)+c9(40*m)+c10(25*m)+c11(20*m)
   67 :   + c12(5*m)+c13(20*m)+c14(25*m)+c15(30*m) );
   68 : //
   69 : //  Plot the mesh.
   70 : //
   71 : plot ( Th, ps = "molding_mesh.ps" );
   72 : //
   73 : real gravity = 9.81;
   74 : real rhom = 1.0;
   75 : real rhop = 20.0;
   76 : real num = 0.1;
   77 : real nup = 10;
   78 : real dt = 1.0;
   79 : //
   80 : //  Define the finite element spaces.
   81 : //
   82 : fespace Uh(Th, P1b);
   83 : fespace Vh(Th, P1);
   84 : fespace Wh(Th, P1b);
   85 : Uh u1, u2, u1old, u2old, u1h, u2h;
   86 : Vh p, ph;
   87 : Wh u1view, u2view;
   88 : Wh rho, rhoold, rho2, mu, muold;
   89 : //
   90 : // Initialization
   91 : //
   92 : rho = rhom;
   93 : mu = rhom * num;
   94 : rho = rho + (rhop-rhom)*((x<0.2)*(y>PG[1])*(y<PF[1]));
   95 : mu = mu + (rhop*nup-rhom*num)*((x<0.2)*(y>PG[1])*(y<PF[1]));
   96 : u1 = 0.0;
   97 : u1old = u1;
   98 : u2 = 0.0;
   99 : u2old = u2;
  100 : //
  101 : plot ( rho, nbiso = 60, fill = 1 );
  102 : //
  103 : //  Define the weak form of the Navier-Stokes equations.
  104 : //
  105 : problem NavierStokes ( [u1, u2, p], [u1h, u2h, ph] ) =
  106 :     int2d (Th) (rho*u1*u1h/dt)
  107 :   - int2d (Th) (rhoold*convect([u1old,u2old], -dt, u1old)*u1h/dt)
  108 :   + int2d (Th) (mu*dx(u1)*dx(u1h)+mu*dy(u1)*dy(u1h))
  109 :   + int2d (Th) (dx(p)*u1h)
  110 :   + int2d (Th) (rho*u2*u2h/dt)
  111 :   - int2d (Th) (rhoold*convect([u1old,u2old], -dt, u2old)*u2h/dt)
  112 :   + int2d (Th) (mu*dx(u2)*dx(u2h)+mu*dy(u2)*dy(u2h))
  113 :   + int2d (Th) (dy(p)*u2h)
  114 :   + int2d (Th) (rho*gravity*u2h)
  115 :   + int2d (Th) ((dx(u1)+dy(u2))*ph+0.000001*p*ph)
  116 :   + on (c6, u1=0.3, u2=0)
  117 :   + on (c1,c2,c4,c5,c7,c8,c9,c10,c11,c12,c13,c14,c15, u1=0, u2=0);
  118 : //
  119 : real[int] viso(1);
  120 : viso = [5, 10, 15];
  121 : 
  122 : for ( int it = 0; it < 40; it++ ) 
  123 : {
  124 :   for ( int sit = 0; sit < 5; sit++ )
  125 :   {
  126 : //
  127 : //  Solve the Navier-Stokes equations.
  128 : //
  129 :     NavierStokes;
  130 : //
  131 : //  Convect the density and dynamic viscosity.
  132 : //
  133 :     u1old = u1;
  134 :     u2old = u2;
  135 :     rhoold = rho;
  136 :     muold = mu;
  137 :     rho = convect ( [ u1old, u2old ], -dt, rho );
  138 :     mu = convect ( [ u1old, u2old ], -dt, mu );
  139 :   }
  140 : // Th = adaptmesh(Th, rho, u1, u2); rho = rho; mu = mu; u1=u1; u2=u2;
  141 :   rho2 = rhop - rho;
  142 :   plot ( rho, fill=1, nbiso=3, viso=viso, value=1, ps = "molding_rho_it"+it+".ps");
  143 : // plot (rho2, nbiso=60, fill =1, grey=1, bw=1, ps="molding_rho_it"+it+".ps");
  144 : }
  145 : //
  146 : //  Terminate.
  147 : //
  148 : cout << "\n";
  149 : cout << "MOLDING:\n";
  150 : cout << "  Normal end of execution.\n";
  151 : 
  152 :  sizestack + 1024 =8424  ( 7400 )

  --  mesh:  Nb of Triangles =  14910, Nb of Vertices 7950
  -- Solve : 
          min -0.161264  max 0.313166
          min -0.504528  max 0.520987
          min 16059.3  max 16776.5
  -- Solve : 
          min -0.347373  max 0.339209
          min -0.769008  max 0.906583
          min 15713  max 16430.1
  -- Solve : 
          min -3.68807  max 2.02522
          min -2.28693  max 3.8414
          min 14511.3  max 15266.4
  -- Solve : 
          min -2.27344  max 1.70414
          min -1.8444  max 2.79446
          min 12801.7  max 13704.6
  -- Solve : 
          min -2.943  max 4.03701
          min -3.35063  max 3.18784
          min 11718.2  max 15075.5
  -- Solve : 
          min -1.0594  max 1.03367
          min -1.25083  max 1.37666
          min 13336  max 14752.6
  -- Solve : 
          min -14.9091  max 9.76833
          min -19.0609  max 5.97461
          min 13842.4  max 15513.4
  -- Solve : 
          min -3.35847  max 4.04991
          min -3.4706  max 5.58283
          min 12394.2  max 16222.3
  -- Solve : 
          min -1.98456  max 2.53922
          min -1.654  max 3.26465
          min 14522.5  max 16229.5
  -- Solve : 
          min -0.328563  max 1.02709
          min -0.704563  max 0.54227
          min 14896.2  max 16606.6
  -- Solve : 
          min -5.72496  max 2.43119
          min -5.46129  max 2.68883
          min 15013.2  max 16800.4
  -- Solve : 
          min -2.50387  max 1.40508
          min -1.9394  max 3.51929
          min 15134.5  max 17032.9
  -- Solve : 
          min -1.90309  max 4.32104
          min -2.46897  max 5.25321
          min 15367.8  max 17307.7
  -- Solve : 
          min -0.896164  max 1.32129
          min -0.654846  max 0.667041
          min 15489.2  max 17442.9
  -- Solve : 
          min -18.8162  max 30.0045
          min -13.3253  max 43.529
          min 15569.6  max 17525.3
  -- Solve : 
          min -0.144398  max 0.325555
          min -0.426401  max 0.348072
          min 15690.3  max 17658.4
  -- Solve : 
          min -0.06653  max 0.451563
          min -0.426167  max 0.317412
          min 15736.7  max 17757.5
  -- Solve : 
          min -0.101846  max 0.60597
          min -0.425791  max 0.318608
          min 15619.8  max 17848.5
  -- Solve : 
          min -0.0663272  max 0.702977
          min -0.429729  max 0.315428
          min 15422.4  max 17883
  -- Solve : 
          min -0.0677045  max 0.392297
          min -0.429016  max 0.316505
          min 15423.4  max 17956.3
  -- Solve : 
          min -0.415588  max 0.341378
          min -0.423741  max 0.32184
          min 15559.1  max 18063.7
  -- Solve : 
          min -0.173023  max 0.317964
          min -0.421978  max 0.323698
          min 15433.3  max 18128.4
  -- Solve : 
          min -0.0966456  max 0.35263
          min -0.421488  max 0.324342
          min 15426.5  max 18165.7
  -- Solve : 
          min -0.111762  max 0.362415
          min -0.422215  max 0.323791
          min 15443.3  max 18196.2
  -- Solve : 
          min -0.0660164  max 0.309003
          min -0.423803  max 0.322434
          min 15452.5  max 18217.6
  -- Solve : 
          min -0.0659204  max 0.309393
          min -0.425614  max 0.320889
          min 15447.6  max 18232.4
  -- Solve : 
          min -0.0658212  max 0.31118
          min -0.427337  max 0.319459
          min 15451.9  max 18245.9
  -- Solve : 
          min -0.0657229  max 0.313166
          min -0.429024  max 0.318073
          min 15462  max 18258.2
  -- Solve : 
          min -0.0656284  max 0.313056
          min -0.430293  max 0.317078
          min 15479.2  max 18274.5
  -- Solve : 
          min -0.0655419  max 0.312489
          min -0.431084  max 0.316516
          min 15497.3  max 18295.2
  -- Solve : 
          min -0.0654609  max 0.310091
          min -0.431574  max 0.316203
          min 15472.6  max 18318.4
  -- Solve : 
          min -0.0653908  max 0.309012
          min -0.431263  max 0.316616
          min 15458.1  max 18348.6
  -- Solve : 
          min -0.0653315  max 0.308028
          min -0.430649  max 0.317264
          min 15440.1  max 18381.5
  -- Solve : 
          min -0.0652791  max 0.306199
          min -0.430054  max 0.317847
          min 15419.3  max 18413.1
  -- Solve : 
          min -0.065226  max 0.305151
          min -0.429212  max 0.318631
          min 15411.7  max 18448.4
  -- Solve : 
          min -0.0651765  max 0.305212
          min -0.427784  max 0.319953
          min 15413.7  max 18491
  -- Solve : 
          min -0.0651339  max 0.306508
          min -0.426026  max 0.321569
          min 15418.2  max 18537.1
  -- Solve : 
          min -0.0650949  max 0.306841
          min -0.424617  max 0.322835
          min 15412.8  max 18577.2
  -- Solve : 
          min -0.0650544  max 0.306221
          min -0.423426  max 0.323888
          min 15395  max 18613.8
  -- Solve : 
          min -0.0650086  max 0.3057
          min -0.422455  max 0.32473
          min 15373.4  max 18647.5
  -- Solve : 
          min -0.0649605  max 0.306251
          min -0.421648  max 0.325418
          min 15372.5  max 18678
  -- Solve : 
          min -0.0649127  max 0.307256
          min -0.420738  max 0.326209
          min 15362.9  max 18708.4
  -- Solve : 
          min -0.0648625  max 0.308347
          min -0.419936  max 0.326894
          min 15363.7  max 18737.6
  -- Solve : 
          min -0.0648107  max 0.308946
          min -0.419176  max 0.327541
          min 15348  max 18765.8
  -- Solve : 
          min -0.06476  max 0.308903
          min -0.418463  max 0.328146
          min 15330.7  max 18792.1
  -- Solve : 
          min -0.0647097  max 0.308507
          min -0.417854  max 0.328654
          min 15320.4  max 18816.4
  -- Solve : 
          min -0.0646616  max 0.307889
          min -0.417231  max 0.32918
          min 15327.1  max 18839.5
  -- Solve : 
          min -0.0646141  max 0.307128
          min -0.416529  max 0.329784
          min 15314.8  max 18863.2
  -- Solve : 
          min -0.0645683  max 0.306333
          min -0.415789  max 0.330423
          min 15294.4  max 18887.2
  -- Solve : 
          min -0.0645225  max 0.305489
          min -0.414971  max 0.331133
          min 15294.2  max 18912.6
  -- Solve : 
          min -0.0644782  max 0.305133
          min -0.413985  max 0.332003
          min 15296.8  max 18939.9
  -- Solve : 
          min -0.0644325  max 0.305131
          min -0.412889  max 0.332977
          min 15286.2  max 18967.5
  -- Solve : 
          min -0.105288  max 0.30513
          min -0.411768  max 0.33398
          min 15270.4  max 18991.5
  -- Solve : 
          min -0.0642785  max 0.30513
          min -0.410718  max 0.334937
          min 15240.1  max 19004.1
  -- Solve : 
          min -0.0641012  max 0.30513
          min -0.409778  max 0.335851
          min 15197.9  max 18992.8
  -- Solve : 
          min -0.063776  max 0.30513
          min -0.408876  max 0.336845
          min 15110.6  max 18943.8
  -- Solve : 
          min -0.063267  max 0.305133
          min -0.408016  max 0.337939
          min 14986.2  max 18850.8
  -- Solve : 
          min -0.0625811  max 0.305136
          min -0.407248  max 0.339078
          min 14791.9  max 18715.8
  -- Solve : 
          min -0.0617583  max 0.305141
          min -0.406626  max 0.340175
          min 14573  max 18548.9
  -- Solve : 
          min -0.0609506  max 0.305145
          min -0.40608  max 0.341172
          min 14394.1  max 18390.2
  -- Solve : 
          min -0.0603549  max 0.305148
          min -0.405509  max 0.342022
          min 14247.7  max 18286.3
  -- Solve : 
          min -0.0600022  max 0.305149
          min -0.404904  max 0.342732
          min 14173.2  max 18240.5
  -- Solve : 
          min -0.0595462  max 0.305151
          min -0.4043  max 0.343562
          min 14058.2  max 18162
  -- Solve : 
          min -0.0587303  max 0.305155
          min -0.403661  max 0.344729
          min 13842.1  max 17993
  -- Solve : 
          min -0.0576429  max 0.305162
          min -0.402984  max 0.346136
          min 13552.8  max 17763
  -- Solve : 
          min -0.0564937  max 0.305168
          min -0.402307  max 0.347589
          min 13276.7  max 17523.6
  -- Solve : 
          min -0.0557826  max 0.305172
          min -0.401656  max 0.348683
          min 13078.9  max 17386.2
  -- Solve : 
          min -0.055534  max 0.305172
          min -0.400985  max 0.349373
          min 13033.9  max 17378.6
  -- Solve : 
          min -0.0548695  max 0.305176
          min -0.400159  max 0.350737
          min 12798.3  max 17215.5
  -- Solve : 
          min -0.0574268  max 0.305175
          min -0.39938  max 0.351459
          min 12765.7  max 17231
  -- Solve : 
          min -0.0854401  max 0.305151
          min -0.399106  max 0.348998
          min 13666.4  max 18065.5
  -- Solve : 
          min -0.441727  max 0.305127
          min -0.556792  max 0.346598
          min 14539.8  max 18868.4
  -- Solve : 
          min -0.090169  max 0.305159
          min -0.397265  max 0.351907
          min 13240.8  max 17749.9
  -- Solve : 
          min -0.064127  max 0.305167
          min -0.396046  max 0.353977
          min 12922.3  max 17476.7
  -- Solve : 
          min -0.0730837  max 0.305183
          min -0.394564  max 0.357335
          min 12263  max 16898.4
  -- Solve : 
          min -0.0801192  max 0.308654
          min -0.391775  max 0.367197
          min 9827.68  max 14754.6
  -- Solve : 
          min -0.0823927  max 0.30502
          min -0.398186  max 0.335176
          min 18572.4  max 22602.9
  -- Solve : 
          min -0.0798528  max 0.30517
          min -0.392606  max 0.358419
          min 12568.4  max 17228.2
  -- Solve : 
          min -0.118042  max 0.30777
          min -0.389726  max 0.368267
          min 10105.8  max 15087.7
  -- Solve : 
          min -0.395248  max 0.607406
          min -2.33252  max 0.290258
          min 30128.8  max 33707.6
  -- Solve : 
          min -0.77264  max 0.305092
          min -0.717081  max 0.555764
          min 15364.9  max 19826
  -- Solve : 
          min -0.507763  max 0.305108
          min -0.393227  max 0.350929
          min 14923.1  max 19459.3
  -- Solve : 
          min -0.290081  max 0.469549
          min -0.369496  max 1.08392
          min -6383.95  max 819.376
  -- Solve : 
          min -0.17654  max 0.344142
          min -0.392332  max 0.354896
          min 14955.2  max 19764.5
  -- Solve : 
          min -0.0922636  max 0.305094
          min -0.392195  max 0.350444
          min 14983.3  max 19831.9
  -- Solve : 
          min -0.097501  max 0.305096
          min -0.39219  max 0.350546
          min 14920.1  max 19789.9
  -- Solve : 
          min -0.105194  max 0.305097
          min -0.392069  max 0.350636
          min 14903.2  max 19779
  -- Solve : 
          min -0.120358  max 0.305097
          min -0.391835  max 0.3508
          min 14879.4  max 19760.6
  -- Solve : 
          min -0.677653  max 0.305098
          min -0.406952  max 1.95961
          min 14867.1  max 19749.9
  -- Solve : 
          min -0.213241  max 0.305098
          min -0.391156  max 0.351429
          min 14860.5  max 19744.7
  -- Solve : 
          min -0.329545  max 0.305098
          min -0.407089  max 0.959201
          min 14862.2  max 19744
  -- Solve : 
          min -0.511676  max 1.08975
          min -1.67121  max 1.0085
          min 14870.9  max 19754.4
  -- Solve : 
          min -0.204419  max 0.305097
          min -0.390155  max 0.40127
          min 14883.5  max 19766.6
  -- Solve : 
          min -0.170224  max 0.305096
          min -0.389876  max 0.443185
          min 14894.4  max 19777.6
  -- Solve : 
          min -0.131129  max 0.305096
          min -0.389592  max 0.36865
          min 14903.3  max 19787.4
  -- Solve : 
          min -0.34569  max 0.313958
          min -0.509609  max 0.353108
          min 14908.7  max 19795.1
  -- Solve : 
          min -0.222606  max 0.305095
          min -0.389045  max 0.357792
          min 14910.7  max 19800
  -- Solve : 
          min -0.109931  max 0.305095
          min -0.388739  max 0.427165
          min 14909.3  max 19803
  -- Solve : 
          min -0.407919  max 0.305095
          min -0.388381  max 0.668571
          min 14906.4  max 19803.2
  -- Solve : 
          min -0.113887  max 0.305094
          min -0.388033  max 0.516581
          min 14902.2  max 19809.3
  -- Solve : 
          min -0.131679  max 0.305094
          min -0.387609  max 0.557846
          min 14900.7  max 19815.1
  -- Solve : 
          min -0.75724  max 0.34316
          min -0.596442  max 1.11633
          min 14897.5  max 19815.1
  -- Solve : 
          min -0.564314  max 0.961327
          min -0.882159  max 0.776507
          min 14891.6  max 19825.2
  -- Solve : 
          min -0.49592  max 0.469916
          min -0.385799  max 1.16428
          min 14889.3  max 19835.3
  -- Solve : 
          min -1.08734  max 0.366685
          min -1.57126  max 1.4238
          min 14882.8  max 19845.2
  -- Solve : 
          min -0.516212  max 0.566476
          min -0.384599  max 0.607389
          min 14884.2  max 19856.2
  -- Solve : 
          min -0.26605  max 0.327933
          min -0.525804  max 0.578816
          min 14878.3  max 19860.7
  -- Solve : 
          min -0.153187  max 0.30509
          min -0.891301  max 0.688997
          min 14874.1  max 19865.4
  -- Solve : 
          min -0.914707  max 0.78676
          min -0.382756  max 0.730206
          min 14869  max 19868.7
  -- Solve : 
          min -0.674243  max 0.462636
          min -0.56014  max 1.02647
          min 14863.8  max 19874.6
  -- Solve : 
          min -2.99706  max 5.48414
          min -5.64528  max 5.42827
          min 14870.5  max 19921.8
  -- Solve : 
          min -0.116546  max 0.305088
          min -0.381244  max 0.540865
          min 14869.2  max 19899.9
  -- Solve : 
          min -0.155826  max 0.305088
          min -0.380671  max 0.361392
          min 14846.2  max 19896.4
  -- Solve : 
          min -0.131624  max 0.305087
          min -0.380282  max 0.361735
          min 14848.1  max 19910.2
  -- Solve : 
          min -0.088436  max 0.305087
          min -0.379881  max 0.362119
          min 14844.1  max 19915.7
  -- Solve : 
          min -0.0866425  max 0.305086
          min -0.379484  max 0.362512
          min 14839.9  max 19917.6
  -- Solve : 
          min -0.0996582  max 0.305086
          min -0.379084  max 0.362915
          min 14835.8  max 19917
  -- Solve : 
          min -0.0877012  max 0.305086
          min -0.378689  max 0.363316
          min 14831.2  max 19915.8
  -- Solve : 
          min -0.116501  max 0.305086
          min -0.378358  max 0.363643
          min 14826.3  max 19917.6
  -- Solve : 
          min -2.11236  max 2.77025
          min -3.04915  max 3.53411
          min 14820.6  max 19924.4
  -- Solve : 
          min -0.209353  max 0.305085
          min -0.37789  max 0.79578
          min 14822.6  max 19927.2
  -- Solve : 
          min -0.258974  max 0.305085
          min -0.377591  max 0.364387
          min 14818.2  max 19927
  -- Solve : 
          min -0.144605  max 0.322924
          min -0.377359  max 0.36461
          min 14814.8  max 19930.4
  -- Solve : 
          min -0.0760774  max 0.305085
          min -0.377163  max 0.364797
          min 14812.8  max 19933.8
  -- Solve : 
          min -0.0719598  max 0.305084
          min -0.39356  max 0.364958
          min 14810.8  max 19937.6
  -- Solve : 
          min -0.0707694  max 0.305084
          min -0.376843  max 0.365096
          min 14809.7  max 19941.4
  -- Solve : 
          min -0.0698574  max 0.305084
          min -0.37671  max 0.365219
          min 14808.2  max 19945.1
  -- Solve : 
          min -0.0690887  max 0.305084
          min -0.376588  max 0.365331
          min 14807.1  max 19948.6
  -- Solve : 
          min -0.068478  max 0.305084
          min -0.376476  max 0.365433
          min 14806  max 19952
  -- Solve : 
          min -0.0679249  max 0.305083
          min -0.37639  max 0.365507
          min 14805.2  max 19956.3
  -- Solve : 
          min -0.0673369  max 0.305083
          min -0.376361  max 0.365514
          min 14806.9  max 19963.4
  -- Solve : 
          min -0.0667212  max 0.305083
          min -0.376336  max 0.36552
          min 14809.1  max 19969.4
  -- Solve : 
          min -0.0661347  max 0.305083
          min -0.376296  max 0.36555
          min 14810.2  max 19973.1
  -- Solve : 
          min -0.0654892  max 0.305083
          min -0.376256  max 0.365581
          min 14810.9  max 19975.9
  -- Solve : 
          min -0.0647609  max 0.305083
          min -0.376221  max 0.365609
          min 14811.5  max 19978.3
  -- Solve : 
          min -0.0640003  max 0.305083
          min -0.376197  max 0.365626
          min 14812.3  max 19980.3
  -- Solve : 
          min -0.0632921  max 0.305083
          min -0.376186  max 0.365632
          min 14813.3  max 19981.9
  -- Solve : 
          min -0.0626713  max 0.305082
          min -0.376183  max 0.365632
          min 14814.4  max 19983.2
  -- Solve : 
          min -0.0621608  max 0.305082
          min -0.376185  max 0.365627
          min 14815.5  max 19984.3
  -- Solve : 
          min -0.0617688  max 0.305082
          min -0.37619  max 0.365619
          min 14816.5  max 19985.2
  -- Solve : 
          min -0.0614814  max 0.305082
          min -0.376198  max 0.365609
          min 14817.4  max 19985.8
  -- Solve : 
          min -0.0612099  max 0.305082
          min -0.376213  max 0.365593
          min 14818.1  max 19986.2
  -- Solve : 
          min -0.0609605  max 0.305082
          min -0.376232  max 0.365573
          min 14818.8  max 19986.5
  -- Solve : 
          min -0.0608556  max 0.305082
          min -0.376255  max 0.36555
          min 14819.3  max 19986.6
  -- Solve : 
          min -0.0607957  max 0.305082
          min -0.376278  max 0.365526
          min 14819.8  max 19986.8
  -- Solve : 
          min -0.0607531  max 0.305082
          min -0.376302  max 0.365502
          min 14820.3  max 19987
  -- Solve : 
          min -0.0607111  max 0.305082
          min -0.376334  max 0.365469
          min 14821  max 19987.2
  -- Solve : 
          min -0.0606933  max 0.305082
          min -0.376383  max 0.36542
          min 14822  max 19987.3
  -- Solve : 
          min -0.0606814  max 0.305083
          min -0.376439  max 0.365363
          min 14823.1  max 19987.3
  -- Solve : 
          min -0.0606713  max 0.305083
          min -0.376499  max 0.365304
          min 14824.1  max 19987.3
  -- Solve : 
          min -0.221763  max 0.400429
          min -0.37655  max 0.365253
          min 14824.9  max 19987.3
  -- Solve : 
          min -0.0606962  max 0.305083
          min -0.376594  max 0.365208
          min 14825.5  max 19987.3
  -- Solve : 
          min -0.0607282  max 0.305083
          min -0.376632  max 0.36517
          min 14826.2  max 19987.3
  -- Solve : 
          min -0.0607609  max 0.305083
          min -0.376671  max 0.365131
          min 14826.9  max 19987.3
  -- Solve : 
          min -0.0607949  max 0.305083
          min -0.376715  max 0.365086
          min 14827.6  max 19987.3
  -- Solve : 
          min -0.0608273  max 0.305083
          min -0.376766  max 0.365036
          min 14828.4  max 19987.2
  -- Solve : 
          min -0.060854  max 0.305083
          min -0.376818  max 0.364983
          min 14829.1  max 19987.1
  -- Solve : 
          min -0.0608751  max 0.305083
          min -0.376869  max 0.364932
          min 14829.9  max 19986.9
  -- Solve : 
          min -0.0608916  max 0.305083
          min -0.376916  max 0.364886
          min 14830.6  max 19986.9
  -- Solve : 
          min -0.0609018  max 0.305083
          min -0.376955  max 0.364846
          min 14831.2  max 19986.8
  -- Solve : 
          min -0.0609039  max 0.305083
          min -0.376987  max 0.364814
          min 14831.7  max 19986.8
  -- Solve : 
          min -0.0608969  max 0.305083
          min -0.377015  max 0.364786
          min 14832.2  max 19986.8
  -- Solve : 
          min -0.0608812  max 0.305083
          min -0.377043  max 0.364758
          min 14832.6  max 19986.7
  -- Solve : 
          min -0.0608583  max 0.305083
          min -0.377071  max 0.364729
          min 14833  max 19986.7
  -- Solve : 
          min -0.0608334  max 0.305083
          min -0.377097  max 0.364703
          min 14833.4  max 19986.6
  -- Solve : 
          min -0.0608125  max 0.305083
          min -0.377122  max 0.364678
          min 14833.7  max 19986.6
  -- Solve : 
          min -0.060796  max 0.305083
          min -0.377147  max 0.364653
          min 14834.1  max 19986.6
  -- Solve : 
          min -0.0607774  max 0.305083
          min -0.377171  max 0.364629
          min 14834.4  max 19986.6
  -- Solve : 
          min -0.0607529  max 0.305083
          min -0.377192  max 0.364608
          min 14834.7  max 19986.6
  -- Solve : 
          min -0.0607205  max 0.305083
          min -0.377209  max 0.36459
          min 14835  max 19986.7
  -- Solve : 
          min -0.0606813  max 0.305083
          min -0.377224  max 0.364575
          min 14835.2  max 19986.7
  -- Solve : 
          min -0.0606649  max 0.305083
          min -0.377238  max 0.36456
          min 14835.4  max 19986.7
  -- Solve : 
          min -0.0606675  max 0.305083
          min -0.377254  max 0.364544
          min 14835.7  max 19986.8
  -- Solve : 
          min -0.0606874  max 0.305083
          min -0.377272  max 0.364526
          min 14835.9  max 19986.8
  -- Solve : 
          min -0.0607222  max 0.305083
          min -0.377295  max 0.364503
          min 14836.3  max 19986.8
  -- Solve : 
          min -0.0607631  max 0.305083
          min -0.377329  max 0.364469
          min 14836.8  max 19986.7
  -- Solve : 
          min -0.0607931  max 0.305083
          min -0.377374  max 0.364425
          min 14837.5  max 19986.6
  -- Solve : 
          min -0.060858  max 0.305083
          min -0.377426  max 0.364373
          min 14838.2  max 19986.4
  -- Solve : 
          min -0.0609739  max 0.305083
          min -0.37748  max 0.364319
          min 14839  max 19986.2
  -- Solve : 
          min -0.0611075  max 0.305083
          min -0.377535  max 0.364266
          min 14839.8  max 19986
  -- Solve : 
          min -0.0612475  max 0.305083
          min -0.377585  max 0.364216
          min 14840.5  max 19985.9
  -- Solve : 
          min -0.0613828  max 0.305083
          min -0.377625  max 0.364176
          min 14841.1  max 19985.7
  -- Solve : 
          min -0.0615074  max 0.305083
          min -0.377654  max 0.364147
          min 14841.5  max 19985.6
  -- Solve : 
          min -0.0616231  max 0.305083
          min -0.377673  max 0.364128
          min 14841.8  max 19985.6
  -- Solve : 
          min -0.0617309  max 0.305083
          min -0.377685  max 0.364115
          min 14842.1  max 19985.7
  -- Solve : 
          min -0.0618345  max 0.305083
          min -0.377693  max 0.364106
          min 14842.3  max 19985.7
  -- Solve : 
          min -0.062009  max 0.305083
          min -0.377698  max 0.364101
          min 14842.4  max 19985.8
  -- Solve : 
          min -0.0621821  max 0.305083
          min -0.377701  max 0.364097
          min 14842.6  max 19985.9
  -- Solve : 
          min -0.0623456  max 0.305083
          min -0.377702  max 0.364096
          min 14842.7  max 19986
  -- Solve : 
          min -0.0624965  max 0.305083
          min -0.3777  max 0.364097
          min 14842.8  max 19986.1
  -- Solve : 
          min -0.0626328  max 0.305083
          min -0.377697  max 0.3641
          min 14842.9  max 19986.2
  -- Solve : 
          min -0.0627529  max 0.305083
          min -0.377693  max 0.364104
          min 14842.9  max 19986.2
  -- Solve : 
          min -0.0628564  max 0.305083
          min -0.377687  max 0.364109
          min 14842.9  max 19986.3
  -- Solve : 
          min -0.062944  max 0.305083
          min -0.377681  max 0.364115
          min 14842.8  max 19986.4
  -- Solve : 
          min -0.0630167  max 0.305083
          min -0.377674  max 0.364122
          min 14842.8  max 19986.4
  -- Solve : 
          min -0.0630762  max 0.305083
          min -0.377667  max 0.364128
          min 14842.7  max 19986.5
  -- Solve : 
          min -0.0631244  max 0.305083
          min -0.37766  max 0.364136
          min 14842.7  max 19986.6
  -- Solve : 
          min -0.063163  max 0.305083
          min -0.377652  max 0.364143
          min 14842.6  max 19986.7
  -- Solve : 
          min -0.0631936  max 0.305083
          min -0.377644  max 0.36415
          min 14842.6  max 19986.8
  -- Solve : 
          min -0.0632175  max 0.305083
          min -0.377637  max 0.364158
          min 14842.5  max 19987

MOLDING:
  Normal end of execution.
times: compile 0.007117s, execution 949.446s,  mpirank:0
 CodeAlloc : nb ptr  4829,  size :524872 mpirank: 0
Ok: Normal End