-- FreeFem++ v4.6 (Thu Apr 2 15:47:38 CEST 2020 - git v4.6) Load: lg_fem lg_mesh lg_mesh3 eigenvalue 1 : // Discussion: 2 : // 3 : // Read the battery_o.msh and battery_g.msh files created by BAMG, 4 : // and save a standard FreeFem++ mesh as battery.msh. 5 : // 6 : // Location: 7 : // 8 : // http://people.sc.fsu.edu/~jburkardt/freefem_src/mitchell_05/battery_convert.edp 9 : // 10 : // Licensing: 11 : // 12 : // This code is distributed under the GNU LGPL license. 13 : // 14 : // Modified: 15 : // 16 : // 30 December 2014 17 : // 18 : // Author: 19 : // 20 : // John Burkardt 21 : // 22 : // Reference: 23 : // 24 : // Frederic Hecht, 25 : // Freefem++, 26 : // Third Edition, version 3.22 27 : // 28 : // William Mitchell, 29 : // A collection of 2D elliptic problems for testing adaptive 30 : // grid refinement algorithms, 31 : // Applied Mathematics and Computation, 32 : // Volume 220, 1 September 2013, pages 350-364. 33 : // 34 : cout << "\n"; 35 : cout << "battery_convert\n"; 36 : cout << " FreeFem++ version\n"; 37 : cout << " Read BAMG battery mesh, create a FreeFem++ mesh ... : .\n"; 38 : // 39 : // Th: the triangulation. 40 : // 41 : mesh Th = readmesh ( "battery_o.msh" ); 42 : // 43 : // Plot the mesh. 44 : // 45 : plot ( Th, wait = 1, ps = "battery.eps" ); 46 : // 47 : // Save the mesh file. 48 : // 49 : savemesh ( Th, "battery.msh" ); 50 : // 51 : // Terminate. 52 : // 53 : cout << "\n"; 54 : cout << "battery_convert\n"; 55 : cout << " Normal end of execution.\n"; 56 : 57 : sizestack + 1024 =1256 ( 232 ) battery_convert FreeFem++ version Read BAMG battery mesh, create a FreeFem++ mesh. -- mesh: Nb of Triangles = 1799, Nb of Vertices 966 number of required edges : 0 battery_convert Normal end of execution. times: compile 0.004264s, execution 0.015546s, mpirank:0 CodeAlloc : nb ptr 3443, size :467736 mpirank: 0 Ok: Normal End -- FreeFem++ v4.6 (Thu Apr 2 15:47:38 CEST 2020 - git v4.6) Load: lg_fem lg_mesh lg_mesh3 eigenvalue 1 : // Discussion: 2 : // 3 : // -d/dx p dudx - d/dy q dudy = f on [0,8.4]x[0,24] 4 : // Mixed boundary conditions 5 : // 6 : // f = given 7 : // g unknown 8 : // 9 : // P, Q and F are piecewise constant functions: 10 : // 11 : // Domain P Q F 12 : // ------ ---- ---- ---- 13 : // 1 25.0 25.0 0.0 14 : // 2 7.0 0.8 1.0 15 : // 3 5.0 0.0001 1.0 16 : // 4 0.2 0.2 0.0 17 : // 5 0.05 0.05 0.0 18 : // 19 : // Y8 +------------------+ 20 : // | K = 1 | 21 : // Y7 +---------------+ | 22 : // | K = 5 | | 23 : // Y6 +--------+---+ | | 24 : // | K = 2 | | | | 25 : // Y5 +--------+ K | | | 26 : // | K = 3 | = | | | 27 : // Y4 +--------+ 4 | | | 28 : // | K = 2 | | | | 29 : // Y3 +--------+ | | | 30 : // | K = 5 | | | | 31 : // Y2 +--------+---+--+ | 32 : // | | 33 : // Y1 +------------------+ 34 : // X1 X2 X3 X4 X5 35 : // 36 : // The rectangular subdomains are characterized by 5 X values: 37 : // 0.0, 6.1, 6.5, 8.0, 8.4 38 : // and 8 Y values: 39 : // 0.0, 0.8, 1.6, 3.6, 18.8, 21.2, 23.2, 24.0 40 : // 41 : // The boundary conditions are of the form: 42 : // p dudx dyds - q dudy dxds + c * u = g 43 : // 44 : // using the following values: 45 : // 46 : // side C gN Boundary condition 47 : // ------- --- --- ------------------ 48 : // bottom 3.0 1.0: - q dudy + 3 u = 1 49 : // right 2.0 2.0: p dudx + 2 u = 2 50 : // top 1.0 3.0: q dudy + u = 3 51 : // left 0.0 0.0: - p dudx = 0 52 : // 53 : // The mesh generated by BAMG will assign the following boundary labels: 54 : // 55 : // side label 56 : // ------- ----- 57 : // bottom 1 58 : // right 2 59 : // top 3 60 : // left 4 61 : // 62 : // Location: 63 : // 64 : // http://people.sc.fsu.edu/~jburkardt/freefem_src/mitchell_05/mitchell_05.edp 65 : // 66 : // Licensing: 67 : // 68 : // This code is distributed under the GNU LGPL license. 69 : // 70 : // Modified: 71 : // 72 : // 10 February 2015 73 : // 74 : // Author: 75 : // 76 : // John Burkardt 77 : // 78 : // Reference: 79 : // 80 : // Frederic Hecht, 81 : // Freefem++, 82 : // Third Edition, version 3.22 83 : // 84 : // William Mitchell, 85 : // A collection of 2D elliptic problems for testing adaptive 86 : // grid refinement algorithms, 87 : // Applied Mathematics and Computation, 88 : // Volume 220, 1 September 2013, pages 350-364. 89 : // 90 : cout << "\n"; 91 : cout << "mitchell_05\n"; 92 : cout << " FreeFem++ version\n"; 93 : cout << " Mitchel test problem 5, the battery.\n"; 94 : // 95 : // Set parameters. 96 : // 97 : real x1 = 0.0; 98 : real x2 = 6.1; 99 : real x3 = 6.5; 100 : real x4 = 8.0; 101 : real x5 = 8.4; 102 : real y1 = 0.0; *** Warning The identifier y1 hide a Global identifier 103 : real y2 = 0.8; 104 : real y3 = 1.6; 105 : real y4 = 3.6; 106 : real y5 = 18.8; 107 : real y6 = 21.2; 108 : real y7 = 23.2; 109 : real y8 = 24.0; 110 : // 111 : // Read the mesh. 112 : // 113 : mesh Th = readmesh ( "battery.msh" ); 114 : // 115 : // Define Vh, the finite element space defined over Th, using P1 basis functions. 116 : // 117 : fespace Vh ( Th, P1 ); 118 : // 119 : // Define U, V, and F, piecewise continuous functions over Th. 120 : // 121 : Vh u; 122 : Vh v; 123 : // 124 : // Define the function F. 125 : // 126 : func real f ( real xx, real yy ) 127 : { 128 : real value; 129 : 130 : value = 0.0; 131 : // 132 : // Region 1: 133 : // 134 : if ( x4 <= xx || y7 <= yy || yy <= y2 ) 135 : { 136 : value = 0.0; 137 : } 138 : // 139 : // Region 2: 140 : // 141 : else if ( y3 <= yy && yy <= y4 && xx <= x2 ) 142 : { 143 : value = 1.0; 144 : } 145 : else if ( y5 <= yy && yy <= y6 && xx <= x2 ) 146 : { 147 : value = 1.0; 148 : } 149 : // 150 : // Region 3; 151 : // 152 : else if ( y4 <= yy && yy <= y5 && xx <= x2 ) 153 : { 154 : value = 1.0; 155 : } 156 : // 157 : // Region 4: 158 : // 159 : else if ( y2 <= yy && yy <= y6 && x2 <= xx && xx <= x3 ) 160 : { 161 : value = 0.0; 162 : } 163 : // 164 : // Region 5: 165 : // 166 : else if ( y6 <= yy && yy <= y7 && xx <= x4 ) 167 : { 168 : value = 0.0; 169 : } 170 : else if ( y2 <= yy && yy <= y7 && x3 <= xx && xx <= x4 ) 171 : { 172 : value = 0.0; 173 : } 174 : else 175 : { 176 : cout << "F: Stray values X = " << x << " Y = " << y << "\n"; 177 : } 178 : return value; 179 : } 180 : // 181 : // Define the boundary condition functions G1, G2, G3, G4: 182 : // 183 : func g1 = 1.0; 184 : func g2 = 2.0; 185 : func g3 = 3.0; 186 : func g4 = 0.0; 187 : // 188 : // Define the coefficient function P. 189 : // 190 : func real p ( real xx, real yy ) 191 : { 192 : real value; 193 : 194 : value = 0.0; 195 : // 196 : // Region 1: 197 : // 198 : if ( x4 <= xx || y7 <= yy || yy <= y2 ) 199 : { 200 : value = 25.0; 201 : } 202 : // 203 : // Region 2: 204 : // 205 : else if ( y3 <= yy && yy <= y4 && xx <= x2 ) 206 : { 207 : value = 7.0; 208 : } 209 : else if ( y5 <= yy && yy <= y6 && xx <= x2 ) 210 : { 211 : value = 7.0; 212 : } 213 : // 214 : // Region 3; 215 : // 216 : else if ( y4 <= yy && yy <= y5 && xx <= x2 ) 217 : { 218 : value = 5.0; 219 : } 220 : // 221 : // Region 4: 222 : // 223 : else if ( y2 <= yy && yy <= y6 && x2 <= xx && xx <= x3 ) 224 : { 225 : value = 0.2; 226 : } 227 : // 228 : // Region 5: 229 : // 230 : else if ( y6 <= yy && yy <= y7 && xx <= x4 ) 231 : { 232 : value = 0.05; 233 : } 234 : else if ( y2 <= yy && yy <= y7 && x3 <= xx && xx <= x4 ) 235 : { 236 : value = 0.05; 237 : } 238 : else if ( y2 <= yy && yy <= y3 && x1 <= xx && xx <= x2 ) 239 : { 240 : value = 0.05; 241 : } 242 : else 243 : { 244 : cout << "P: Stray values: X = " << x << " Y = " << y << "\n"; 245 : } 246 : return value; 247 : } 248 : // 249 : // Define the coefficient function Q. 250 : // 251 : func real q ( real xx, real yy ) 252 : { 253 : real value; 254 : 255 : value = 0.0; 256 : // 257 : // Region 1: 258 : // 259 : if ( x4 <= xx || y7 <= yy || yy <= y2 ) 260 : { 261 : value = 25.0; 262 : } 263 : // 264 : // Region 2: 265 : // 266 : else if ( y3 <= yy && yy <= y4 && xx <= x2 ) 267 : { 268 : value = 0.8; 269 : } 270 : else if ( y5 <= yy && yy <= y6 && xx <= x2 ) 271 : { 272 : value = 0.8; 273 : } 274 : // 275 : // Region 3; 276 : // 277 : else if ( y4 <= yy && yy <= y5 && xx <= x2 ) 278 : { 279 : value = 0.0001; 280 : } 281 : // 282 : // Region 4: 283 : // 284 : else if ( y2 <= yy && yy <= y6 && x2 <= xx && xx <= x3 ) 285 : { 286 : value = 0.2; 287 : } 288 : // 289 : // Region 5: 290 : // Two regions, one of which is described by two rectangles. 291 : // 292 : else if ( y6 <= yy && yy <= y7 && xx <= x4 ) 293 : { 294 : value = 0.05; 295 : } 296 : else if ( y2 <= yy && yy <= y7 && x3 <= xx && xx <= x4 ) 297 : { 298 : value = 0.05; 299 : } 300 : else if ( y2 <= yy && yy <= y3 && x1 <= xx && xx <= x2 ) 301 : { 302 : value = 0.05; 303 : } 304 : else 305 : { 306 : cout << "Q: Stray values: X = " << x << " Y = " << y << "\n"; 307 : } 308 : return value; 309 : } 310 : // 311 : // Solve the variational problem with Robin boundary conditions. 312 : // 313 : solve Laplace ( u, v ) 314 : = int2d ( Th ) ( p(x,y) * dx(u)*dx(v) + q(x,y) * dy(u)*dy(v) ) 315 : - int2d ( Th ) ( f(x,y) * v ) 316 : - int1d ( Th, 1 ) ( 3.0 * u * v ) 317 : - int1d ( Th, 1 ) ( - 1.0 * v ) 318 : - int1d ( Th, 2 ) ( 2.0 * u * v ) 319 : - int1d ( Th, 2 ) ( - 2.0 * v ) 320 : - int1d ( Th, 3 ) ( u * v ) 321 : - int1d ( Th, 3 ) ( - 3.0 * v ); 322 : // 323 : // Plot the solution. 324 : // 325 : plot ( u, wait = 1, fill = true, ps = "mitchell_05_u.eps" ); 326 : // 327 : // Plot the mesh. 328 : // 329 : plot ( Th, wait = 1, ps = "mitchell_05_mesh.eps" ); 330 : // 331 : // Save the mesh file. 332 : // 333 : savemesh ( Th, "mitchell_05.msh" ); 334 : // 335 : // Terminate. 336 : // 337 : cout << "\n"; 338 : cout << "mitchell_05:\n"; 339 : cout << " Normal end of execution.\n"; 340 : 341 : sizestack + 1024 =3248 ( 2224 ) mitchell_05 FreeFem++ version Mitchel test problem 5, the battery. -- mesh: Nb of Triangles = 1799, Nb of Vertices 966 F: Stray values X = 0.338889 Y = 1.06667 F: Stray values X = 0.456846 Y = 1.15949 F: Stray values X = 0.456846 Y = 0.881029 F: Stray values X = 0.102975 Y = 1.15949 F: Stray values X = 0.269344 Y = 1.01194 F: Stray values X = 0.269344 Y = 1.17611 F: Stray values X = 0.477978 Y = 1.01194 F: Stray values X = 0.169444 Y = 0.933333 F: Stray values X = 0.0514873 Y = 1.11897 F: Stray values X = 0.405359 Y = 0.840515 F: Stray values X = 0.0514873 Y = 0.840515 F: Stray values X = 0.238989 Y = 0.823886 F: Stray values X = 0.0303556 Y = 0.988057 F: Stray values X = 0.238989 Y = 0.988057 F: Stray values X = 5.93056 Y = 0.933333 F: Stray values X = 5.69464 Y = 0.840515 F: Stray values X = 6.04851 Y = 1.11897 F: Stray values X = 6.04851 Y = 0.840515 F: Stray values X = 6.06964 Y = 0.988057 F: Stray values X = 5.86101 Y = 0.823886 F: Stray values X = 5.86101 Y = 0.988057 F: Stray values X = 5.93056 Y = 1.46667 F: Stray values X = 6.04851 Y = 1.55949 F: Stray values X = 6.04851 Y = 1.28103 F: Stray values X = 5.69464 Y = 1.55949 F: Stray values X = 5.86101 Y = 1.41194 F: Stray values X = 5.86101 Y = 1.57611 F: Stray values X = 6.06964 Y = 1.41194 F: Stray values X = 5.76111 Y = 1.33333 F: Stray values X = 5.64315 Y = 1.51897 F: Stray values X = 5.99703 Y = 1.24051 F: Stray values X = 5.64315 Y = 1.24051 F: Stray values X = 5.83066 Y = 1.22389 F: Stray values X = 5.62202 Y = 1.38806 F: Stray values X = 5.83066 Y = 1.38806 F: Stray values X = 5.25278 Y = 1.33333 F: Stray values X = 5.48869 Y = 1.24051 F: Stray values X = 5.13482 Y = 1.24051 F: Stray values X = 5.13482 Y = 1.51897 F: Stray values X = 5.11369 Y = 1.38806 F: Stray values X = 5.32232 Y = 1.38806 F: Stray values X = 5.32232 Y = 1.22389 F: Stray values X = 5.76111 Y = 1.06667 F: Stray values X = 5.99703 Y = 1.15949 F: Stray values X = 5.64315 Y = 0.881029 F: Stray values X = 5.64315 Y = 1.15949 F: Stray values X = 5.62202 Y = 1.01194 F: Stray values X = 5.83066 Y = 1.17611 F: Stray values X = 5.83066 Y = 1.01194 F: Stray values X = 5.42222 Y = 1.46667 F: Stray values X = 5.54018 Y = 1.55949 F: Stray values X = 5.54018 Y = 1.28103 F: Stray values X = 5.18631 Y = 1.55949 F: Stray values X = 5.35268 Y = 1.41194 F: Stray values X = 5.35268 Y = 1.57611 F: Stray values X = 5.56131 Y = 1.41194 F: Stray values X = 5.42222 Y = 0.933333 F: Stray values X = 5.18631 Y = 0.840515 F: Stray values X = 5.54018 Y = 1.11897 F: Stray values X = 5.54018 Y = 0.840515 F: Stray values X = 5.56131 Y = 0.988057 F: Stray values X = 5.35268 Y = 0.823886 F: Stray values X = 5.35268 Y = 0.988057 F: Stray values X = 4.74444 Y = 1.46667 F: Stray values X = 4.98036 Y = 1.55949 F: Stray values X = 4.62649 Y = 1.28103 F: Stray values X = 4.62649 Y = 1.55949 F: Stray values X = 4.60536 Y = 1.41194 F: Stray values X = 4.81399 Y = 1.57611 F: Stray values X = 4.81399 Y = 1.41194 F: Stray values X = 4.91389 Y = 1.06667 F: Stray values X = 5.03185 Y = 1.15949 F: Stray values X = 5.03185 Y = 0.881029 F: Stray values X = 4.67797 Y = 1.15949 F: Stray values X = 4.84434 Y = 1.01194 F: Stray values X = 4.84434 Y = 1.17611 F: Stray values X = 5.05298 Y = 1.01194 F: Stray values X = 5.25278 Y = 1.06667 F: Stray values X = 5.48869 Y = 1.15949 F: Stray values X = 5.13482 Y = 0.881029 F: Stray values X = 5.13482 Y = 1.15949 F: Stray values X = 5.11369 Y = 1.01194 F: Stray values X = 5.32232 Y = 1.17611 F: Stray values X = 5.32232 Y = 1.01194 F: Stray values X = 4.91389 Y = 1.33333 F: Stray values X = 5.03185 Y = 1.51897 F: Stray values X = 5.03185 Y = 1.24051 F: Stray values X = 4.67797 Y = 1.24051 F: Stray values X = 4.84434 Y = 1.22389 F: Stray values X = 4.84434 Y = 1.38806 F: Stray values X = 5.05298 Y = 1.38806 F: Stray values X = 4.40556 Y = 1.06667 F: Stray values X = 4.52351 Y = 1.15949 F: Stray values X = 4.52351 Y = 0.881029 F: Stray values X = 4.16964 Y = 1.15949 F: Stray values X = 4.33601 Y = 1.01194 F: Stray values X = 4.33601 Y = 1.17611 F: Stray values X = 4.54464 Y = 1.01194 F: Stray values X = 4.74444 Y = 0.933333 F: Stray values X = 4.98036 Y = 0.840515 F: Stray values X = 4.62649 Y = 0.840515 F: Stray values X = 4.62649 Y = 1.11897 F: Stray values X = 4.60536 Y = 0.988057 F: Stray values X = 4.81399 Y = 0.988057 F: Stray values X = 4.81399 Y = 0.823886 F: Stray values X = 4.40556 Y = 1.33333 F: Stray values X = 4.52351 Y = 1.51897 F: Stray values X = 4.52351 Y = 1.24051 F: Stray values X = 4.16964 Y = 1.24051 F: Stray values X = 4.33601 Y = 1.22389 F: Stray values X = 4.33601 Y = 1.38806 F: Stray values X = 4.54464 Y = 1.38806 F: Stray values X = 4.23611 Y = 0.933333 F: Stray values X = 4.11815 Y = 1.11897 F: Stray values X = 4.47203 Y = 0.840515 F: Stray values X = 4.11815 Y = 0.840515 F: Stray values X = 4.30566 Y = 0.823886 F: Stray values X = 4.09702 Y = 0.988057 F: Stray values X = 4.30566 Y = 0.988057 F: Stray values X = 4.23611 Y = 1.46667 F: Stray values X = 4.11815 Y = 1.55949 F: Stray values X = 4.47203 Y = 1.55949 F: Stray values X = 4.11815 Y = 1.28103 F: Stray values X = 4.30566 Y = 1.41194 F: Stray values X = 4.09702 Y = 1.41194 F: Stray values X = 4.30566 Y = 1.57611 F: Stray values X = 3.77014 Y = 1.46667 F: Stray values X = 3.97656 Y = 1.55949 F: Stray values X = 3.71116 Y = 1.28103 F: Stray values X = 3.62269 Y = 1.55949 F: Stray values X = 3.64844 Y = 1.41194 F: Stray values X = 3.80491 Y = 1.57611 F: Stray values X = 3.85707 Y = 1.41194 F: Stray values X = 3.51597 Y = 1.33333 F: Stray values X = 3.63393 Y = 1.24051 F: Stray values X = 3.36853 Y = 1.24051 F: Stray values X = 3.54546 Y = 1.51897 F: Stray values X = 3.44643 Y = 1.38806 F: Stray values X = 3.6029 Y = 1.38806 F: Stray values X = 3.49859 Y = 1.22389 F: Stray values X = 3.93958 Y = 1.06667 F: Stray values X = 3.76265 Y = 1.15949 F: Stray values X = 4.02805 Y = 1.15949 F: Stray values X = 4.02805 Y = 0.881029 F: Stray values X = 4.0439 Y = 1.01194 F: Stray values X = 3.88742 Y = 1.01194 F: Stray values X = 3.88742 Y = 1.17611 F: Stray values X = 3.30417 Y = 1.46667 F: Stray values X = 3.12723 Y = 1.55949 F: Stray values X = 3.4811 Y = 1.55949 F: Stray values X = 3.30417 Y = 1.28103 F: Stray values X = 3.40848 Y = 1.41194 F: Stray values X = 3.19985 Y = 1.41194 F: Stray values X = 3.30417 Y = 1.57611 F: Stray values X = 3.05 Y = 1.33333 F: Stray values X = 3.22694 Y = 1.24051 F: Stray values X = 2.87306 Y = 1.24051 F: Stray values X = 3.05 Y = 1.51897 F: Stray values X = 2.94568 Y = 1.38806 F: Stray values X = 3.15432 Y = 1.38806 F: Stray values X = 3.05 Y = 1.22389 F: Stray values X = 3.30417 Y = 0.933333 F: Stray values X = 3.4811 Y = 0.840515 F: Stray values X = 3.12723 Y = 0.840515 F: Stray values X = 3.30417 Y = 1.11897 F: Stray values X = 3.19985 Y = 0.988057 F: Stray values X = 3.40848 Y = 0.988057 F: Stray values X = 3.30417 Y = 0.823886 F: Stray values X = 3.93958 Y = 1.33333 F: Stray values X = 3.76265 Y = 1.24051 F: Stray values X = 4.02805 Y = 1.51897 F: Stray values X = 4.02805 Y = 1.24051 F: Stray values X = 4.0439 Y = 1.38806 F: Stray values X = 3.88742 Y = 1.22389 F: Stray values X = 3.88742 Y = 1.38806 F: Stray values X = 3.51597 Y = 1.06667 F: Stray values X = 3.54546 Y = 0.881029 F: Stray values X = 3.36853 Y = 1.15949 F: Stray values X = 3.63393 Y = 1.15949 F: Stray values X = 3.49859 Y = 1.17611 F: Stray values X = 3.6029 Y = 1.01194 F: Stray values X = 3.44643 Y = 1.01194 F: Stray values X = 3.77014 Y = 0.933333 F: Stray values X = 3.62269 Y = 0.840515 F: Stray values X = 3.71116 Y = 1.11897 F: Stray values X = 3.97656 Y = 0.840515 F: Stray values X = 3.85707 Y = 0.988057 F: Stray values X = 3.80491 Y = 0.823886 F: Stray values X = 3.64844 Y = 0.988057 F: Stray values X = 3.05 Y = 1.06667 F: Stray values X = 3.05 Y = 0.881029 F: Stray values X = 2.87306 Y = 1.15949 F: Stray values X = 3.22694 Y = 1.15949 F: Stray values X = 3.05 Y = 1.17611 F: Stray values X = 3.15432 Y = 1.01194 F: Stray values X = 2.94568 Y = 1.01194 F: Stray values X = 2.79583 Y = 1.46667 F: Stray values X = 2.79583 Y = 1.28103 F: Stray values X = 2.6189 Y = 1.55949 F: Stray values X = 2.97277 Y = 1.55949 F: Stray values X = 2.79583 Y = 1.57611 F: Stray values X = 2.90015 Y = 1.41194 F: Stray values X = 2.69152 Y = 1.41194 F: Stray values X = 2.54167 Y = 1.33333 F: Stray values X = 2.36473 Y = 1.24051 F: Stray values X = 2.54167 Y = 1.51897 F: Stray values X = 2.7186 Y = 1.24051 F: Stray values X = 2.64598 Y = 1.38806 F: Stray values X = 2.54167 Y = 1.22389 F: Stray values X = 2.43735 Y = 1.38806 F: Stray values X = 2.79583 Y = 0.933333 F: Stray values X = 2.6189 Y = 0.840515 F: Stray values X = 2.79583 Y = 1.11897 F: Stray values X = 2.97277 Y = 0.840515 F: Stray values X = 2.90015 Y = 0.988057 F: Stray values X = 2.79583 Y = 0.823886 F: Stray values X = 2.69152 Y = 0.988057 F: Stray values X = 2.2875 Y = 1.46667 F: Stray values X = 2.2875 Y = 1.28103 F: Stray values X = 2.11056 Y = 1.55949 F: Stray values X = 2.46444 Y = 1.55949 F: Stray values X = 2.2875 Y = 1.57611 F: Stray values X = 2.39182 Y = 1.41194 F: Stray values X = 2.18318 Y = 1.41194 F: Stray values X = 2.54167 Y = 1.06667 F: Stray values X = 2.36473 Y = 1.15949 F: Stray values X = 2.7186 Y = 1.15949 F: Stray values X = 2.54167 Y = 0.881029 F: Stray values X = 2.64598 Y = 1.01194 F: Stray values X = 2.43735 Y = 1.01194 F: Stray values X = 2.54167 Y = 1.17611 F: Stray values X = 2.03333 Y = 1.33333 F: Stray values X = 1.8564 Y = 1.24051 F: Stray values X = 2.03333 Y = 1.51897 F: Stray values X = 2.21027 Y = 1.24051 F: Stray values X = 2.13765 Y = 1.38806 F: Stray values X = 2.03333 Y = 1.22389 F: Stray values X = 1.92902 Y = 1.38806 F: Stray values X = 2.2875 Y = 0.933333 F: Stray values X = 2.11056 Y = 0.840515 F: Stray values X = 2.2875 Y = 1.11897 F: Stray values X = 2.46444 Y = 0.840515 F: Stray values X = 2.39182 Y = 0.988057 F: Stray values X = 2.2875 Y = 0.823886 F: Stray values X = 2.18318 Y = 0.988057 F: Stray values X = 1.77917 Y = 1.46667 F: Stray values X = 1.77917 Y = 1.28103 F: Stray values X = 1.60223 Y = 1.55949 F: Stray values X = 1.9561 Y = 1.55949 F: Stray values X = 1.77917 Y = 1.57611 F: Stray values X = 1.88348 Y = 1.41194 F: Stray values X = 1.67485 Y = 1.41194 F: Stray values X = 1.31319 Y = 1.46667 F: Stray values X = 1.37217 Y = 1.28103 F: Stray values X = 1.10677 Y = 1.55949 F: Stray values X = 1.46064 Y = 1.55949 F: Stray values X = 1.27842 Y = 1.57611 F: Stray values X = 1.4349 Y = 1.41194 F: Stray values X = 1.22626 Y = 1.41194 F: Stray values X = 1.77917 Y = 0.933333 F: Stray values X = 1.9561 Y = 0.840515 F: Stray values X = 1.60223 Y = 0.840515 F: Stray values X = 1.77917 Y = 1.11897 F: Stray values X = 1.67485 Y = 0.988057 F: Stray values X = 1.88348 Y = 0.988057 F: Stray values X = 1.77917 Y = 0.823886 F: Stray values X = 2.03333 Y = 1.06667 F: Stray values X = 2.03333 Y = 0.881029 F: Stray values X = 1.8564 Y = 1.15949 F: Stray values X = 2.21027 Y = 1.15949 F: Stray values X = 2.03333 Y = 1.17611 F: Stray values X = 2.13765 Y = 1.01194 F: Stray values X = 1.92902 Y = 1.01194 F: Stray values X = 1.56736 Y = 1.33333 F: Stray values X = 1.4494 Y = 1.24051 F: Stray values X = 1.53787 Y = 1.51897 F: Stray values X = 1.71481 Y = 1.24051 F: Stray values X = 1.63691 Y = 1.38806 F: Stray values X = 1.58475 Y = 1.22389 F: Stray values X = 1.48043 Y = 1.38806 F: Stray values X = 1.14375 Y = 1.06667 F: Stray values X = 1.05528 Y = 0.881029 F: Stray values X = 1.05528 Y = 1.15949 F: Stray values X = 1.32069 Y = 1.15949 F: Stray values X = 1.19591 Y = 1.17611 F: Stray values X = 1.19591 Y = 1.01194 F: Stray values X = 1.03943 Y = 1.01194 F: Stray values X = 1.56736 Y = 1.06667 F: Stray values X = 1.53787 Y = 0.881029 F: Stray values X = 1.4494 Y = 1.15949 F: Stray values X = 1.71481 Y = 1.15949 F: Stray values X = 1.58475 Y = 1.17611 F: Stray values X = 1.63691 Y = 1.01194 F: Stray values X = 1.48043 Y = 1.01194 F: Stray values X = 1.14375 Y = 1.33333 F: Stray values X = 1.05528 Y = 1.24051 F: Stray values X = 1.05528 Y = 1.51897 F: Stray values X = 1.32069 Y = 1.24051 F: Stray values X = 1.19591 Y = 1.38806 F: Stray values X = 1.19591 Y = 1.22389 F: Stray values X = 1.03943 Y = 1.38806 F: Stray values X = 0.847222 Y = 1.46667 F: Stray values X = 0.611308 Y = 1.55949 F: Stray values X = 0.965179 Y = 1.55949 F: Stray values X = 0.965179 Y = 1.28103 F: Stray values X = 0.986311 Y = 1.41194 F: Stray values X = 0.777678 Y = 1.41194 F: Stray values X = 0.777678 Y = 1.57611 F: Stray values X = 1.31319 Y = 0.933333 F: Stray values X = 1.10677 Y = 0.840515 F: Stray values X = 1.37217 Y = 1.11897 F: Stray values X = 1.46064 Y = 0.840515 F: Stray values X = 1.4349 Y = 0.988057 F: Stray values X = 1.27842 Y = 0.823886 F: Stray values X = 1.22626 Y = 0.988057 F: Stray values X = 0.677778 Y = 1.06667 F: Stray values X = 0.559821 Y = 1.15949 F: Stray values X = 0.913692 Y = 1.15949 F: Stray values X = 0.559821 Y = 0.881029 F: Stray values X = 0.747322 Y = 1.01194 F: Stray values X = 0.538689 Y = 1.01194 F: Stray values X = 0.747322 Y = 1.17611 F: Stray values X = 0.677778 Y = 1.33333 F: Stray values X = 0.559821 Y = 1.51897 F: Stray values X = 0.913692 Y = 1.24051 F: Stray values X = 0.559821 Y = 1.24051 F: Stray values X = 0.747322 Y = 1.22389 F: Stray values X = 0.538689 Y = 1.38806 F: Stray values X = 0.747322 Y = 1.38806 F: Stray values X = 0.338889 Y = 1.33333 F: Stray values X = 0.456846 Y = 1.51897 F: Stray values X = 0.456846 Y = 1.24051 F: Stray values X = 0.102975 Y = 1.24051 F: Stray values X = 0.269344 Y = 1.22389 F: Stray values X = 0.269344 Y = 1.38806 F: Stray values X = 0.477978 Y = 1.38806 F: Stray values X = 0.847222 Y = 0.933333 F: Stray values X = 0.965179 Y = 0.840515 F: Stray values X = 0.611308 Y = 0.840515 F: Stray values X = 0.965179 Y = 1.11897 F: Stray values X = 0.777678 Y = 0.988057 F: Stray values X = 0.986311 Y = 0.988057 F: Stray values X = 0.777678 Y = 0.823886 F: Stray values X = 0.169444 Y = 1.46667 F: Stray values X = 0.0514873 Y = 1.28103 F: Stray values X = 0.0514873 Y = 1.55949 F: Stray values X = 0.405359 Y = 1.55949 F: Stray values X = 0.238989 Y = 1.57611 F: Stray values X = 0.238989 Y = 1.41194 F: Stray values X = 0.0303556 Y = 1.41194 -- Solve : min -8.37751 max 192.294 number of required edges : 0 mitchell_05: Normal end of execution. times: compile 0.004716s, execution 0.795058s, mpirank:0 CodeAlloc : nb ptr 4062, size :493760 mpirank: 0 Ok: Normal End