-- FreeFem++ v4.14 (mer. 06 mars 2024 16:59:04 CET - git v4.14-1-g2b2052ae)
   file : molding.edp
 Load: lg_fem lg_mesh lg_mesh3 eigenvalue 
    1 : //  molding.edp
    2 : //
    3 : //  Discussion:
    4 : //
    5 : //    Simulation of injection molding.
    6 : //
    7 : //  Licensing:
    8 : //
    9 : //    This code is distributed under the MIT license.
   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 : cout << "\n";
   28 : cout << "molding():\n";
   29 : cout << "  FreeFem++ version\n";
   30 : cout << "  Fill a hollow mold with a liquid.\n";
   31 : 
   32 : real[int] PA(2); PA = [0, 0];
   33 : real[int] PB(2); PB = [9, 0];
   34 : real[int] PC(2); PC = [9, 9];
   35 : real[int] PD(2); PD = [8, 9];
   36 : real[int] PE(2); PE = [0, 9];
   37 : real [int ] PF(2); PF = [0, 6];
   38 : real [int ] PG(2); PG = [0,4];
   39 : real [int ] PH(2); PH = [1, 1];
   40 : real [int ] PK(2); PK = [1,8];
   41 : real [int ] PL(2); PL = [8,8];
   42 : real [int ] PM(2); PM = [8,5];
   43 : real [int ] PO(2); PO = [5,5];
   44 : real [int ] PP(2); PP = [5,4];
   45 : real [int ] PQ(2); PQ = [8,4];
   46 : real [int ] PR(2); PR = [8,1];
   47 : 
   48 : border c1(t=0,1) {x =(1-t)*PA[0]+t*PB[0]; y =(1-t)*PA[1]+t*PB[1];}
   49 : border c2(t=0,1) {x =(1-t)*PB[0]+t*PC[0]; y =(1-t)*PB[1]+t*PC[1];}
   50 : border c3(t=0,1) {x =(1-t)*PC[0]+t*PD[0]; y =(1-t)*PC[1]+t*PD[1];}
   51 : border c4(t=0,1) {x =(1-t)*PD[0]+t*PE[0]; y =(1-t)*PD[1]+t*PE[1];}
   52 : border c5(t=0,1) {x =(1-t)*PE[0]+t*PF[0]; y =(1-t)*PE[1]+t*PF[1];}
   53 : border c6(t=0,1) {x =(1-t)*PF[0]+t*PG[0]; y =(1-t)*PF[1]+t*PG[1];}
   54 : border c7(t=0,1) {x =(1-t)*PG[0]+t*PA[0]; y =(1-t)*PG[1]+t*PA[1];}
   55 : //
   56 : border c8(t=0,1){x=(1-t)*PH[0]+t*PK[0]; y=(1-t)*PH[1]+t*PK[1];}
   57 : border c9(t=0,1){x=(1-t)*PK[0]+t*PL[0]; y=(1-t)*PK[1]+t*PL[1];}
   58 : border c10(t=0,1){x=(1-t)*PL[0]+t*PM[0]; y=(1-t)*PL[1]+t*PM[1];}
   59 : border c11(t=0,1){x=(1-t)*PM[0]+t*PO[0]; y=(1-t)*PM[1]+t*PO[1];}
   60 : border c12(t=0,1){x=(1-t)*PO[0]+t*PP[0]; y=(1-t)*PO[1]+t*PP[1];}
   61 : border c13(t=0,1){x=(1-t)*PP[0]+t*PQ[0]; y=(1-t)*PP[1]+t*PQ[1];}
   62 : border c14(t=0,1){x=(1-t)*PQ[0]+t*PR[0]; y=(1-t)*PQ[1]+t*PR[1];}
   63 : border c15(t=0,1){x=(1-t)*PR[0]+t*PH[0]; y=(1-t)*PR[1]+t*PH[1];}
   64 : //
   65 : //  Define the mesh, Th.
   66 : //
   67 : int m = 2;
   68 : 
   69 : mesh Th = buildmesh ( c1(80*m)+c2(50*m)+c3(10*m)+c4(40*m)
   70 :   + c5(30*m)+c6(10*m)+c7(30*m)
   71 :   + c8(80*m)+c9(40*m)+c10(25*m)+c11(20*m)
   72 :   + c12(5*m)+c13(20*m)+c14(25*m)+c15(30*m) );
   73 : //
   74 : //  Plot the mesh.
   75 : //
   76 : plot ( Th, ps = "molding_mesh.ps" );
   77 : //
   78 : real gravity = 9.81;
   79 : real rhom = 1.0;
   80 : real rhop = 20.0;
   81 : real num = 0.1;
   82 : real nup = 10;
   83 : real dt = 1.0;
   84 : //
   85 : //  Define the finite element spaces.
   86 : //
   87 : fespace Uh(Th, P1b);
   88 : fespace Vh(Th, P1);
   89 : fespace Wh(Th, P1b);
   90 : Uh u1, u2, u1old, u2old, u1h, u2h;
   91 : Vh p, ph;
   92 : Wh u1view, u2view;
   93 : Wh rho, rhoold, rho2, mu, muold;
   94 : //
   95 : // Initialization
   96 : //
   97 : rho = rhom;
   98 : mu = rhom * num;
   99 : rho = rho + (rhop-rhom)*((x<0.2)*(y>PG[1])*(y<PF[1]));
  100 : mu = mu + (rhop*nup-rhom*num)*((x<0.2)*(y>PG[1])*(y<PF[1]));
  101 : u1 = 0.0;
  102 : u1old = u1;
  103 : u2 = 0.0;
  104 : u2old = u2;
  105 : //
  106 : plot ( rho, nbiso = 60, fill = true );
  107 : //
  108 : //  Define the weak form of the Navier-Stokes equations.
  109 : //
  110 : problem NavierStokes ( [u1, u2, p], [u1h, u2h, ph] ) =
  111 :     int2d (Th) (rho*u1*u1h/dt)
  112 :   - int2d (Th) (rhoold*convect([u1old,u2old], -dt, u1old)*u1h/dt)
  113 :   + int2d (Th) (mu*dx(u1)*dx(u1h)+mu*dy(u1)*dy(u1h))
  114 :   + int2d (Th) (dx(p)*u1h)
  115 :   + int2d (Th) (rho*u2*u2h/dt)
  116 :   - int2d (Th) (rhoold*convect([u1old,u2old], -dt, u2old)*u2h/dt)
  117 :   + int2d (Th) (mu*dx(u2)*dx(u2h)+mu*dy(u2)*dy(u2h))
  118 :   + int2d (Th) (dy(p)*u2h)
  119 :   + int2d (Th) (rho*gravity*u2h)
  120 :   + int2d (Th) ((dx(u1)+dy(u2))*ph+0.000001*p*ph)
  121 :   + on (c6, u1=0.3, u2=0)
  122 :   + on (c1,c2,c4,c5,c7,c8,c9,c10,c11,c12,c13,c14,c15, u1=0, u2=0);
  123 : //
  124 : real[int] viso(1);
  125 : viso = [5, 10, 15];
  126 : 
  127 : for ( int it = 0; it < 40; it++ ) 
  128 : {
  129 :   for ( int sit = 0; sit < 5; sit++ )
  130 :   {
  131 : //
  132 : //  Solve the Navier-Stokes equations.
  133 : //
  134 :     NavierStokes;
  135 : //
  136 : //  Convect the density and dynamic viscosity.
  137 : //
  138 :     u1old = u1;
  139 :     u2old = u2;
  140 :     rhoold = rho;
  141 :     muold = mu;
  142 :     rho = convect ( [ u1old, u2old ], -dt, rho );
  143 :     mu = convect ( [ u1old, u2old ], -dt, mu );
  144 :   }
  145 : // Th = adaptmesh(Th, rho, u1, u2);
  146 : // rho = rho;
  147 : // mu = mu;
  148 : // u1=u1;
  149 : // u2=u2;
  150 : 
  151 :   rho2 = rhop - rho;
  152 : 
  153 :   plot ( rho, fill = true, nbiso = 3, viso = viso, value = true, 
  154 :     ps = "molding_rho_it"+it+".ps");
  155 : 
  156 : // plot ( rho2, nbiso = 60, fill = true, grey = true, bw = true,
  157 : //   ps = "molding_rho_it"+it+".ps" );
  158 : }
  159 : //
  160 : //  Terminate.
  161 : //
  162 : cout << "\n";
  163 : cout << "molding():\n";
  164 : cout << "  Normal end of execution.\n";
  165 : 
  166 : exit ( 0 );
  167 :  sizestack + 1024 =8424  ( 7400 )


molding():
  FreeFem++ version
  Fill a hollow mold with a liquid.
  --  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.70413
          min -1.8444  max 2.79445
          min 12801.7  max 13704.6
  -- Solve : 
          min -2.94273  max 4.03663
          min -3.35037  max 3.18752
          min 11718.4  max 15075.3
  -- Solve : 
          min -1.21132  max 1.21835
          min -1.25076  max 1.37688
          min 13336.1  max 14752.6
  -- Solve : 
          min -27.1512  max 18.1634
          min -34.8758  max 10.6112
          min 13840.3  max 15512.2
  -- Solve : 
          min -9.59556  max 2.62207
          min -1.79879  max 23.2575
          min 13994.4  max 15812.5
  -- Solve : 
          min -1.73474  max 2.605
          min -1.3351  max 0.820602
          min 14558.5  max 16263.7
  -- Solve : 
          min -0.429132  max 1.10034
          min -0.887061  max 0.852611
          min 14915.4  max 16624.6
  -- Solve : 
          min -2.07988  max 1.76994
          min -4.79061  max 1.06198
          min 14928.5  max 16792
  -- Solve : 
          min -80.4378  max 62.9699
          min -38.4417  max 72.787
          min 15117.6  max 17230.8
  -- Solve : 
          min -0.263572  max 0.858751
          min -0.459116  max 0.507618
          min 15362.9  max 17305.5
  -- Solve : 
          min -0.3953  max 0.406848
          min -0.425273  max 0.480729
          min 15512.9  max 17473.6
  -- Solve : 
          min -2.93913  max 2.1776
          min -3.46199  max 0.874523
          min 15639.4  max 17612
  -- Solve : 
          min -1.28869  max 2.05326
          min -1.65355  max 2.52113
          min 15713.7  max 17707.2
  -- Solve : 
          min -3.02497  max 1.65014
          min -1.75333  max 4.39648
          min 15750.6  max 17793.3
  -- Solve : 
          min -0.0838752  max 0.621589
          min -0.422743  max 0.321215
          min 15617.6  max 17891.2
  -- Solve : 
          min -0.0675227  max 0.664569
          min -0.425602  max 0.3191
          min 15449.7  max 17940.4
  -- Solve : 
          min -0.721728  max 1.40478
          min -0.711227  max 0.393685
          min 15463.8  max 18024.2
  -- Solve : 
          min -0.0662687  max 0.345013
          min -0.419653  max 0.325531
          min 15547.7  max 18118.1
  -- Solve : 
          min -0.0662535  max 0.344037
          min -0.419484  max 0.325911
          min 15455.7  max 18159.6
  -- Solve : 
          min -0.299363  max 0.487015
          min -0.419518  max 0.326088
          min 15444.4  max 18191.1
  -- Solve : 
          min -0.146119  max 0.692329
          min -0.421032  max 0.324812
          min 15456  max 18214.1
  -- Solve : 
          min -0.0660349  max 0.319734
          min -0.422808  max 0.323294
          min 15470.5  max 18237.1
  -- Solve : 
          min -0.0659453  max 0.315452
          min -0.423508  max 0.322849
          min 15470.5  max 18260.8
  -- Solve : 
          min -0.0658593  max 0.311331
          min -0.424655  max 0.321974
          min 15475.5  max 18280
  -- Solve : 
          min -0.0657695  max 0.308712
          min -0.425927  max 0.320967
          min 15483.7  max 18297.6
  -- Solve : 
          min -0.065681  max 0.308048
          min -0.42723  max 0.319918
          min 15497.5  max 18313.1
  -- Solve : 
          min -0.0655951  max 0.3089
          min -0.428292  max 0.319076
          min 15512.8  max 18329.7
  -- Solve : 
          min -0.0655116  max 0.307686
          min -0.428928  max 0.318607
          min 15491.2  max 18351
  -- Solve : 
          min -0.0654384  max 0.306744
          min -0.428961  max 0.318684
          min 15474.5  max 18377.4
  -- Solve : 
          min -0.0653752  max 0.305417
          min -0.428468  max 0.319219
          min 15456  max 18409.4
  -- Solve : 
          min -0.0653193  max 0.305151
          min -0.427997  max 0.319692
          min 15432.7  max 18440.5
  -- Solve : 
          min -0.0652636  max 0.30515
          min -0.427309  max 0.320341
          min 15435.4  max 18474.7
  -- Solve : 
          min -0.0652121  max 0.305182
          min -0.426198  max 0.321373
          min 15429.7  max 18513.3
  -- Solve : 
          min -0.065167  max 0.306231
          min -0.424917  max 0.322548
          min 15429  max 18552.5
  -- Solve : 
          min -0.065123  max 0.306149
          min -0.423626  max 0.32372
          min 15422.3  max 18591.1
  -- Solve : 
          min -0.0650792  max 0.305761
          min -0.422792  max 0.32445
          min 15405.8  max 18622.9
  -- Solve : 
          min -0.0650288  max 0.306049
          min -0.422101  max 0.325038
          min 15380.9  max 18653.2
  -- Solve : 
          min -0.0649768  max 0.306941
          min -0.42133  max 0.325701
          min 15377.4  max 18683.4
  -- Solve : 
          min -0.0649266  max 0.30805
          min -0.42053  max 0.326393
          min 15374.6  max 18712.5
  -- Solve : 
          min -0.0648741  max 0.309032
          min -0.419849  max 0.326969
          min 15373.9  max 18740.1
  -- Solve : 
          min -0.0648212  max 0.309383
          min -0.419243  max 0.327477
          min 15350.1  max 18765.7
  -- Solve : 
          min -0.0647677  max 0.309099
          min -0.418552  max 0.328066
          min 15329.6  max 18791.7
  -- Solve : 
          min -0.0647171  max 0.308514
          min -0.417887  max 0.328632
          min 15331.1  max 18816.3
  -- Solve : 
          min -0.0646689  max 0.307822
          min -0.417232  max 0.329191
          min 15329.4  max 18839.8
  -- Solve : 
          min -0.0646228  max 0.30713
          min -0.416577  max 0.32975
          min 15315.3  max 18862.8
  -- Solve : 
          min -0.0645773  max 0.306475
          min -0.415932  max 0.330298
          min 15296.6  max 18885.9
  -- Solve : 
          min -0.0645309  max 0.305681
          min -0.415128  max 0.330996
          min 15304.4  max 18911.3
  -- Solve : 
          min -0.0644859  max 0.305133
          min -0.414169  max 0.331841
          min 15297.6  max 18938.3
  -- Solve : 
          min -0.0644362  max 0.305132
          min -0.413161  max 0.332735
          min 15286  max 18963.7
  -- Solve : 
          min -0.0643637  max 0.305131
          min -0.412125  max 0.33367
          min 15266.7  max 18983.3
  -- Solve : 
          min -0.0642373  max 0.30513
          min -0.411102  max 0.334633
          min 15226.8  max 18987.8
  -- Solve : 
          min -0.0639986  max 0.305131
          min -0.410176  max 0.335586
          min 15164.5  max 18961.1
  -- Solve : 
          min -1.17175  max 0.478312
          min -0.489606  max 0.807469
          min 15063.5  max 18890.7
  -- Solve : 
          min -0.0629978  max 0.305135
          min -0.408657  max 0.337572
          min 14912.3  max 18775.8
  -- Solve : 
          min -0.062221  max 0.30514
          min -0.408086  max 0.338597
          min 14700.9  max 18616.6
  -- Solve : 
          min -0.0613561  max 0.305144
          min -0.407643  max 0.339553
          min 14478  max 18439.2
  -- Solve : 
          min -0.060609  max 0.305148
          min -0.407234  max 0.340371
          min 14313.6  max 18295.5
  -- Solve : 
          min -0.0601489  max 0.30515
          min -0.406779  max 0.341009
          min 14193.3  max 18223
  -- Solve : 
          min -0.0598189  max 0.305151
          min -0.406257  max 0.34163
          min 14119.5  max 18178.6
  -- Solve : 
          min -0.059202  max 0.305155
          min -0.405709  max 0.342537
          min 13965.3  max 18056.7
  -- Solve : 
          min -0.0583037  max 0.30516
          min -0.405113  max 0.343715
          min 13721.1  max 17868.9
  -- Solve : 
          min -0.0572123  max 0.305166
          min -0.404437  max 0.345117
          min 13438.8  max 17639.8
  -- Solve : 
          min -0.0563235  max 0.305171
          min -0.403692  max 0.346442
          min 13219.5  max 17459.7
  -- Solve : 
          min -0.0558704  max 0.305172
          min -0.402938  max 0.347396
          min 13105.9  max 17395.2
  -- Solve : 
          min -0.0556793  max 0.305172
          min -0.402151  max 0.348176
          min 13063.8  max 17395.1
  -- Solve : 
          min -0.0547938  max 0.305178
          min -0.401192  max 0.349908
          min 12747.3  max 17159.7
  -- Solve : 
          min -0.0565163  max 0.305172
          min -0.400393  max 0.34998
          min 12941.3  max 17375.3
  -- Solve : 
          min -0.101923  max 0.305138
          min -0.400166  max 0.346468
          min 14197.1  max 18525.1
  -- Solve : 
          min -0.124992  max 0.305131
          min -0.399495  max 0.346417
          min 14423.9  max 18759.1
  -- Solve : 
          min -0.0778824  max 0.355477
          min -0.394836  max 0.369838
          min 7952.29  max 13064.4
  -- Solve : 
          min -0.155531  max 0.305087
          min -0.398371  max 0.342481
          min 16113.2  max 20286.9
  -- Solve : 
          min -2.12343  max 3.98876
          min -1.85041  max 2.17132
          min 14682.6  max 19033.2
  -- Solve : 
          min -0.066708  max 0.305113
          min -0.396336  max 0.347743
          min 14996.6  max 19331.9
  -- Solve : 
          min -0.0786872  max 0.315076
          min -0.394099  max 0.355068
          min 13191.8  max 17746.5
  -- Solve : 
          min -0.164834  max 0.312988
          min -0.826461  max 0.384985
          min 9357.56  max 14380.1
  -- Solve : 
          min -0.0655094  max 0.342721
          min -0.394931  max 0.347083
          min 15461.2  max 19911.7
  -- Solve : 
          min -0.073376  max 0.305047
          min -0.396269  max 0.340553
          min 17348.5  max 21560.8
  -- Solve : 
          min -0.351421  max 0.654425
          min -0.399778  max 2.21214
          min -20080.5  max -9831
  -- Solve : 
          min -0.329114  max 0.305122
          min -0.562526  max 0.539668
          min 14292.3  max 19109.2
  -- Solve : 
          min -0.62373  max 0.99047
          min -0.392188  max 0.352065
          min 14377.2  max 19269.8
  -- Solve : 
          min -0.0802145  max 0.305109
          min -0.392167  max 0.351928
          min 14426.2  max 19330.7
  -- Solve : 
          min -0.0839487  max 0.305108
          min -0.392077  max 0.351764
          min 14496.9  max 19403.4
  -- Solve : 
          min -0.0872885  max 0.305106
          min -0.391908  max 0.351607
          min 14577.4  max 19482.4
  -- Solve : 
          min -0.0900597  max 0.305104
          min -0.391613  max 0.35167
          min 14630.1  max 19540
  -- Solve : 
          min -0.0927312  max 0.305103
          min -0.391143  max 0.352006
          min 14650.8  max 19569.1
  -- Solve : 
          min -0.0962253  max 0.305102
          min -0.390617  max 0.352475
          min 14650.1  max 19583.4
  -- Solve : 
          min -0.235699  max 0.305102
          min -0.390076  max 0.352984
          min 14643.5  max 19594.5
  -- Solve : 
          min -0.103509  max 0.305101
          min -0.38955  max 0.353479
          min 14637.2  max 19605.5
  -- Solve : 
          min -0.106487  max 0.305101
          min -0.389024  max 0.353976
          min 14632.7  max 19614.9
  -- Solve : 
          min -6.87081  max 2.58548
          min -2.94093  max 2.64612
          min 13423.7  max 19632.8
  -- Solve : 
          min -0.12171  max 0.305099
          min -0.388118  max 0.50907
          min 14634.4  max 19637.1
  -- Solve : 
          min -0.150669  max 0.305099
          min -0.387646  max 0.651545
          min 14629.8  max 19638.9
  -- Solve : 
          min -0.285724  max 0.305099
          min -0.387268  max 0.756385
          min 14629.5  max 19645
  -- Solve : 
          min -0.241326  max 0.305098
          min -0.386861  max 0.935213
          min 14630.5  max 19654.2
  -- Solve : 
          min -0.28493  max 0.305098
          min -0.3864  max 0.796095
          min 14631.3  max 19664.7
  -- Solve : 
          min -0.13365  max 0.305097
          min -0.385854  max 0.369243
          min 14630.5  max 19675.6
  -- Solve : 
          min -0.95125  max 0.628799
          min -1.37505  max 0.863752
          min 14628.6  max 19685.8
  -- Solve : 
          min -0.111597  max 0.305096
          min -0.384731  max 0.357995
          min 14625.6  max 19694.5
  -- Solve : 
          min -0.108955  max 0.305095
          min -0.384217  max 0.358484
          min 14622.4  max 19701.5
  -- Solve : 
          min -0.107553  max 0.305095
          min -0.383741  max 0.358937
          min 14619.6  max 19708.1
  -- Solve : 
          min -0.112931  max 0.305094
          min -0.383309  max 0.359349
          min 14617  max 19714.2
  -- Solve : 
          min -0.1161  max 0.305094
          min -0.382915  max 0.359725
          min 14614.7  max 19720.1
  -- Solve : 
          min -0.16469  max 0.305094
          min -0.382585  max 0.360035
          min 14613.5  max 19726.2
  -- Solve : 
          min -0.222423  max 0.305093
          min -0.382278  max 0.360321
          min 14613.4  max 19733.5
  -- Solve : 
          min -0.66632  max 0.37365
          min -1.02288  max 0.449002
          min 14614.1  max 19742.2
  -- Solve : 
          min -0.135874  max 0.305092
          min -0.381682  max 0.360864
          min 14615  max 19750.1
  -- Solve : 
          min -0.11184  max 0.305092
          min -0.381449  max 0.36107
          min 14617  max 19759
  -- Solve : 
          min -0.103711  max 0.305092
          min -0.381216  max 0.361285
          min 14617  max 19764.6
  -- Solve : 
          min -0.0989206  max 0.305091
          min -0.381041  max 0.361432
          min 14619.9  max 19773.7
  -- Solve : 
          min -0.09568  max 0.305091
          min -0.380948  max 0.361482
          min 14628  max 19787.4
  -- Solve : 
          min -0.0937069  max 0.305091
          min -0.380863  max 0.361533
          min 14634.7  max 19798.5
  -- Solve : 
          min -0.276215  max 0.30509
          min -0.38076  max 0.361614
          min 14638.1  max 19805.7
  -- Solve : 
          min -0.0950338  max 0.30509
          min -0.380639  max 0.361718
          min 14640.1  max 19811.1
  -- Solve : 
          min -0.235583  max 0.542031
          min -0.380509  max 0.683664
          min 14641.6  max 19815.8
  -- Solve : 
          min -0.409516  max 0.33465
          min -0.91947  max 0.727779
          min 14642.7  max 19819.4
  -- Solve : 
          min -0.440312  max 0.30509
          min -0.571345  max 0.561903
          min 14643.4  max 19822.3
  -- Solve : 
          min -1.14465  max 0.68436
          min -0.706184  max 0.600582
          min 14643.5  max 19824.7
  -- Solve : 
          min -0.136389  max 0.305089
          min -0.380032  max 0.362276
          min 14643.1  max 19826.5
  -- Solve : 
          min -0.0955048  max 0.305089
          min -0.379921  max 0.362381
          min 14642.6  max 19828
  -- Solve : 
          min -0.0815921  max 0.305089
          min -0.379818  max 0.362481
          min 14642  max 19829.4
  -- Solve : 
          min -0.080933  max 0.305089
          min -0.379718  max 0.362577
          min 14641.3  max 19830.6
  -- Solve : 
          min -0.0805685  max 0.305089
          min -0.379622  max 0.36267
          min 14640.6  max 19831.9
  -- Solve : 
          min -0.0804055  max 0.305089
          min -0.379523  max 0.362764
          min 14640  max 19833.2
  -- Solve : 
          min -0.0803689  max 0.305089
          min -0.379421  max 0.362862
          min 14639.5  max 19834.7
  -- Solve : 
          min -0.0801226  max 0.305089
          min -0.379334  max 0.362945
          min 14639.4  max 19836.1
  -- Solve : 
          min -0.0801644  max 0.305089
          min -0.379266  max 0.363009
          min 14639.4  max 19837.2
  -- Solve : 
          min -0.0814747  max 0.305088
          min -0.379201  max 0.363072
          min 14639.2  max 19838.1
  -- Solve : 
          min -0.232266  max 1.20472
          min -1.03522  max 1.38841
          min 14637.9  max 19839.7
  -- Solve : 
          min -0.0770606  max 0.339929
          min -0.379113  max 0.363157
          min 14638.6  max 19839.2
  -- Solve : 
          min -0.0758372  max 0.354313
          min -0.379068  max 0.363199
          min 14638.2  max 19839.8
  -- Solve : 
          min -0.0739734  max 0.305088
          min -0.379016  max 0.36325
          min 14638.1  max 19840.8
  -- Solve : 
          min -5.39151  max 2.14397
          min -5.05331  max 5.41112
          min 13648.7  max 19830.4
  -- Solve : 
          min -0.0731944  max 0.305088
          min -0.379106  max 0.363167
          min 14639.8  max 19841.4
  -- Solve : 
          min -0.0705817  max 0.305088
          min -0.379031  max 0.363233
          min 14640.3  max 19842.8
  -- Solve : 
          min -0.0691296  max 0.305088
          min -0.378981  max 0.363274
          min 14640  max 19843.1
  -- Solve : 
          min -0.067972  max 0.305088
          min -0.378946  max 0.363309
          min 14639.6  max 19843.4
  -- Solve : 
          min -0.0668861  max 0.305088
          min -0.378916  max 0.363338
          min 14639.2  max 19843.7
  -- Solve : 
          min -0.0823872  max 0.305088
          min -0.378889  max 0.363365
          min 14638.9  max 19844.2
  -- Solve : 
          min -0.0647642  max 0.305088
          min -0.378863  max 0.36339
          min 14638.6  max 19844.7
  -- Solve : 
          min -0.0638874  max 0.305088
          min -0.378837  max 0.363414
          min 14638.5  max 19845.3
  -- Solve : 
          min -0.0632437  max 0.305088
          min -0.378814  max 0.363436
          min 14638.5  max 19845.9
  -- Solve : 
          min -0.0629365  max 0.305088
          min -0.378797  max 0.363451
          min 14638.6  max 19846.5
  -- Solve : 
          min -0.0626942  max 0.305088
          min -0.378783  max 0.363464
          min 14638.6  max 19846.9
  -- Solve : 
          min -0.0624963  max 0.305088
          min -0.378767  max 0.363479
          min 14638.5  max 19847.3
  -- Solve : 
          min -0.0623307  max 0.305088
          min -0.378749  max 0.363496
          min 14638.4  max 19847.7
  -- Solve : 
          min -0.0622091  max 0.305088
          min -0.378732  max 0.363511
          min 14638.4  max 19848.2
  -- Solve : 
          min -0.559556  max 0.570677
          min -0.818119  max 0.363525
          min 14638.6  max 19848.8
  -- Solve : 
          min -0.0620952  max 0.305088
          min -0.378698  max 0.363541
          min 14638.9  max 19849.3
  -- Solve : 
          min -0.19977  max 0.305088
          min -0.378682  max 0.363556
          min 14639  max 19849.7
  -- Solve : 
          min -0.0736012  max 0.305088
          min -0.378661  max 0.363577
          min 14639  max 19850.1
  -- Solve : 
          min -0.0619097  max 0.305088
          min -0.378639  max 0.363597
          min 14638.8  max 19850.3
  -- Solve : 
          min -0.0618103  max 0.305088
          min -0.378618  max 0.363618
          min 14638.6  max 19850.5
  -- Solve : 
          min -0.0617058  max 0.305088
          min -0.378597  max 0.363639
          min 14638.3  max 19850.7
  -- Solve : 
          min -0.447866  max 0.416373
          min -0.378581  max 0.363654
          min 14638.1  max 19850.9
  -- Solve : 
          min -0.0614611  max 0.305088
          min -0.378564  max 0.36367
          min 14637.9  max 19851.2
  -- Solve : 
          min -0.0613672  max 0.305088
          min -0.378555  max 0.363679
          min 14637.9  max 19851.4
  -- Solve : 
          min -0.0612897  max 0.305088
          min -0.378552  max 0.363681
          min 14637.9  max 19851.6
  -- Solve : 
          min -0.193191  max 0.653632
          min -0.378544  max 0.363689
          min 14638  max 19851.9
  -- Solve : 
          min -0.0612695  max 0.305088
          min -0.378568  max 0.363665
          min 14638.1  max 19851.7
  -- Solve : 
          min -0.0612457  max 0.305088
          min -0.378581  max 0.363652
          min 14638.3  max 19851.7
  -- Solve : 
          min -0.0612882  max 0.305088
          min -0.378594  max 0.363639
          min 14638.5  max 19851.8
  -- Solve : 
          min -0.0613678  max 0.305088
          min -0.378602  max 0.36363
          min 14638.6  max 19852
  -- Solve : 
          min -0.0614209  max 0.305088
          min -0.378605  max 0.363627
          min 14638.8  max 19852.1
  -- Solve : 
          min -0.0614528  max 0.305088
          min -0.378603  max 0.363628
          min 14638.8  max 19852.3
  -- Solve : 
          min -0.0614718  max 0.305088
          min -0.378598  max 0.363632
          min 14638.9  max 19852.6
  -- Solve : 
          min -0.0614789  max 0.305088
          min -0.378589  max 0.36364
          min 14638.9  max 19852.8
  -- Solve : 
          min -0.0614836  max 0.305088
          min -0.378578  max 0.36365
          min 14638.9  max 19853.1
  -- Solve : 
          min -0.0614819  max 0.305088
          min -0.378564  max 0.363663
          min 14638.9  max 19853.3
  -- Solve : 
          min -0.0614706  max 0.305088
          min -0.378548  max 0.363678
          min 14638.8  max 19853.6
  -- Solve : 
          min -0.0614528  max 0.305088
          min -0.378531  max 0.363695
          min 14638.8  max 19853.9
  -- Solve : 
          min -0.0614302  max 0.305088
          min -0.378512  max 0.363712
          min 14638.7  max 19854.2
  -- Solve : 
          min -0.0614039  max 0.305088
          min -0.378493  max 0.36373
          min 14638.7  max 19854.6
  -- Solve : 
          min -0.061375  max 0.305088
          min -0.378473  max 0.363749
          min 14638.6  max 19854.9
  -- Solve : 
          min -0.0613443  max 0.305088
          min -0.378453  max 0.363768
          min 14638.5  max 19855.2
  -- Solve : 
          min -0.0613126  max 0.305088
          min -0.378433  max 0.363787
          min 14638.5  max 19855.5
  -- Solve : 
          min -0.0612807  max 0.305087
          min -0.378413  max 0.363806
          min 14638.4  max 19855.8
  -- Solve : 
          min -0.0612491  max 0.305087
          min -0.378394  max 0.363825
          min 14638.3  max 19856.1
  -- Solve : 
          min -0.061218  max 0.305087
          min -0.378375  max 0.363843
          min 14638.2  max 19856.4
  -- Solve : 
          min -0.0611877  max 0.305087
          min -0.378356  max 0.363861
          min 14638.1  max 19856.6
  -- Solve : 
          min -0.0611581  max 0.305087
          min -0.378338  max 0.363878
          min 14638  max 19856.9
  -- Solve : 
          min -0.0611294  max 0.305087
          min -0.378321  max 0.363895
          min 14638  max 19857.1
  -- Solve : 
          min -0.0611016  max 0.305087
          min -0.378305  max 0.36391
          min 14637.9  max 19857.4
  -- Solve : 
          min -0.0610746  max 0.305087
          min -0.378289  max 0.363925
          min 14637.9  max 19857.7
  -- Solve : 
          min -0.0610485  max 0.305087
          min -0.378274  max 0.363939
          min 14637.8  max 19857.9
  -- Solve : 
          min -0.0610232  max 0.305087
          min -0.37826  max 0.363953
          min 14637.8  max 19858.1
  -- Solve : 
          min -0.060999  max 0.305087
          min -0.378246  max 0.363966
          min 14637.8  max 19858.3
  -- Solve : 
          min -0.0609757  max 0.305087
          min -0.378233  max 0.363979
          min 14637.7  max 19858.5
  -- Solve : 
          min -0.0609538  max 0.305087
          min -0.378221  max 0.363991
          min 14637.6  max 19858.7
  -- Solve : 
          min -0.0609333  max 0.305087
          min -0.378209  max 0.364002
          min 14637.6  max 19858.8
  -- Solve : 
          min -0.0609144  max 0.305087
          min -0.378198  max 0.364013
          min 14637.5  max 19859
  -- Solve : 
          min -0.0608974  max 0.305087
          min -0.378187  max 0.364023
          min 14637.5  max 19859.1
  -- Solve : 
          min -0.0608822  max 0.305087
          min -0.378177  max 0.364033
          min 14637.4  max 19859.3
  -- Solve : 
          min -0.0608688  max 0.305087
          min -0.378167  max 0.364042
          min 14637.4  max 19859.4
  -- Solve : 
          min -0.0608573  max 0.305087
          min -0.378158  max 0.36405
          min 14637.3  max 19859.6
  -- Solve : 
          min -0.0608474  max 0.305087
          min -0.37815  max 0.364059
          min 14637.3  max 19859.7
  -- Solve : 
          min -0.0608391  max 0.305087
          min -0.378142  max 0.364066
          min 14637.3  max 19859.8
  -- Solve : 
          min -0.0608322  max 0.305087
          min -0.378134  max 0.364074
          min 14637.2  max 19860

molding():
  Normal end of execution.
  current line = 166
exit(0)
 err code 0 ,  mpirank 0
 CodeAlloc : nb ptr  5207,  size :573216 mpirank: 0
Ok: Normal End