<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div class="">Dear Petsc Developer team,</div><div class=""><br class=""></div><div class="">I solved a linear elastic problem in 3D using a DMPLEX. My system is converging, then I would like to write out my solution vector to a vtk file where I use unstructured mesh. Currently, I tried two algorithms and I have the same result.</div><div class=""><br class=""></div><div class="">1) Algorithm 1 </div><div class="">err = SNESSolve(_snes, bc_vec_test, solution);</div><div class="">CHKERRABORT(FOX::Parallel::COMM_WORLD,err);</div><div class="">PetscViewer vtk; </div><div class="">PetscViewerVTKOpen(FOX::Parallel::COMM_WORLD,"solution.vtk",FILE_MODE_WRITE,&vtk); </div><div class="">VecView(solution,vtk);</div><div class="">PetscViewerDestroy(&vtk);</div><div class=""><br class=""></div><div class=""><br class=""></div><div class="">2) Algorithm 2</div><div class="">err = SNESSolve(_snes, bc_vec_test, solution);</div><div class="">CHKERRABORT(FOX::Parallel::COMM_WORLD,err);</div><div class="">PetscViewer vtk; </div><div class="">PetscViewerCreate(FOX::Parallel::COMM_WORLD, &vtk); </div><div class="">PetscViewerSetType(vtk, PETSCVIEWERVTK); </div><div class="">PetscViewerFileSetName(vtk, "sol.vtk"); </div><div class="">VecView(solution, vtk); </div><div class="">PetscViewerDestroy(&vtk);</div><div class=""> </div><div class="">The result seems correct except for the rotation order of the nodes (<b class="">see the red lines on gmsh and vtk file respectively</b>). Then, I visualized my vtk file with paraview, and I remarked that my geometry is not correct and not conserved when comparing it with my gmsh file. So, I didn’t understand why the rotation order of nodes is not conserved when saving my result to a vtk file?</div><div class=""> </div><div class="">Other information used:</div><div class="">- gmsh format 2.2 </div><div class="">- Vtk version: 7.1.1</div><div class="">- Petsc version: 3.13/opt </div><div class=""><br class=""></div><div class="">Below my two files gmsh and vtk:</div><div class=""><br class=""></div><div class=""><b class="">Gmsh file:</b></div><div class="">$MeshFormat</div><div class="">2.2 0 8</div><div class="">$EndMeshFormat</div><div class="">$Nodes</div><div class="">12</div><div class="">1 0.0 10.0 10.0</div><div class="">2 0.0 0.0 10.0</div><div class="">3 0.0 0.0 0.0</div><div class="">4 0.0 10.0 0.0</div><div class="">5 10.0 10.0 10.0</div><div class="">6 10.0 0.0 10.0</div><div class="">7 10.0 0.0 0.0</div><div class="">8 10.0 10.0 0.0</div><div class="">9 20.0 10.0 10.0</div><div class="">10 20.0 0.0 10.0</div><div class="">11 20.0 0.0 0.0</div><div class="">12 20.0 10.0 0.0</div><div class="">$EndNodes</div><div class="">$Elements</div><div class="">2</div><div class="">1 5 2 68 60 <font color="#ff2600" class="">1 2 3 4 5 6 7 8</font></div><div class="">2 5 2 68 60 <font color="#ff2600" class="">5 6 7 8 9 10 11 12</font></div><div class="">$EndElements</div><div class=""><br class=""></div><div class=""><b class="">Vtk file :</b></div><div class=""># vtk DataFile Version 2.0</div><div class="">Simplicial Mesh Example</div><div class="">ASCII</div><div class="">DATASET UNSTRUCTURED_GRID</div><div class="">POINTS 12 double</div><div class="">0.000000e+00 1.000000e+01 1.000000e+01</div><div class="">0.000000e+00 0.000000e+00 1.000000e+01</div><div class="">0.000000e+00 0.000000e+00 0.000000e+00</div><div class="">0.000000e+00 1.000000e+01 0.000000e+00</div><div class="">1.000000e+01 1.000000e+01 1.000000e+01</div><div class="">1.000000e+01 0.000000e+00 1.000000e+01</div><div class="">1.000000e+01 0.000000e+00 0.000000e+00</div><div class="">1.000000e+01 1.000000e+01 0.000000e+00</div><div class="">2.000000e+01 1.000000e+01 1.000000e+01</div><div class="">2.000000e+01 0.000000e+00 1.000000e+01</div><div class="">2.000000e+01 0.000000e+00 0.000000e+00</div><div class="">2.000000e+01 1.000000e+01 0.000000e+00</div><div class="">CELLS 2 18</div><div class="">8 <font color="#ff2600" class="">0 3 2 1 4 5 6 7</font></div><div class="">8 <font color="#ff2600" class="">4 7 6 5 8 9 10 11</font></div><div class="">CELL_TYPES 2</div><div class="">12</div><div class="">12</div><div class="">POINT_DATA 12</div><div class="">VECTORS dU_x double</div><div class="">2.754808e-10 -8.653846e-11 -8.653846e-11</div><div class="">2.754808e-10 8.653846e-11 -8.653846e-11</div><div class="">2.754808e-10 8.653846e-11 8.653846e-11</div><div class="">2.754808e-10 -8.653846e-11 8.653846e-11</div><div class="">4.678571e-01 -9.107143e-02 -9.107143e-02</div><div class="">4.678571e-01 9.107143e-02 -9.107143e-02</div><div class="">4.678571e-01 9.107143e-02 9.107143e-02</div><div class="">4.678571e-01 -9.107143e-02 9.107143e-02</div><div class="">1.000000e+00 -7.500000e-02 -7.500000e-02</div><div class="">1.000000e+00 7.500000e-02 -7.500000e-02</div><div class="">1.000000e+00 7.500000e-02 7.500000e-02</div><div class="">1.000000e+00 -7.500000e-02 7.500000e-02</div><div class=""><br class=""></div><div class="">Thank you in advance and have a good day !</div><div class=""><br class=""></div><div class="">Sami,</div><div class=""><br class=""></div><div class="">
<meta charset="UTF-8" class=""><div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">--</div><div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Dr. Sami BEN ELHAJ SALAH<br class="">Ingénieur de Recherche (CNRS)<br class="">Institut Pprime - ISAE - ENSMA<br class="">Mobile: 06.62.51.26.74<br class=""><a href="mailto:sami.ben-elhaj-salah@ensma.fr" class="">Email: sami.ben-elhaj-salah@ensma.fr</a></div><div dir="auto" style="caret-color: rgb(0, 0, 0); color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><a href="https://samiben91.github.io/samibenelhajsalah/index.html" class="">www.samibenelhajsalah.com</a><br class=""><br class=""><br class=""></div></div></div></div>
</div>
<br class=""></body></html>