Dear Mrs Marchandise,<br><br>Following our previous mail discussion, and after reading of the paper "High Quality Surface Remeshing Using Harmonic Maps", I have some additionnal questions and I would greatly appreciate any help for a better understanding of confusing subjects to me.<br>


<br>1) As you have been saying in a previous mail about the outputs of the surface remeshing process:<br>
>>  No the parametrization is not a kind of B-spline surface but a
 planar surface (the surface has been flattened during the 
parametrization process)<br>
If you remesh with the option -saveall, you will see a bunch of files that show you the mapping u(x))<br>
see for example UVX-idSurf.pos that shows the X coordinates of the initial surface on the parametric plane UV.<br>
<br>I would like to know what kind of surface is actually shown on the right of the pictures "1_surface1.png" & "2_plane1.png" enclosed here, and extracted from the paper "High Quality Surface Remeshing Using Harmonic Maps" (shapes displaying no mesh)? Are they in this case very fine meshes obtained after the surface division step into surfaces of genus = 0?<br>


<br>2) I have also a question on the surface division step. The surface shown on picture "1_surface1.png" presents if I understand correctly a genus g=0.<br>However, the subdivision returns finally 3 surfaces.<br>
I guess that the subdivision step also split the surface into smaller ones according sharp edges identification, am I wrong?<br>
Then, the outer surface of the shape displayed on picture 1 has a genus 1 and should be split. Is that right?<br>However, I seem to miss something because as highlighted on picture "3_plane2.png", the parametrization appears to be able to handle surfaces with holes.<br>
If till now, my understanding is correct, why couldn't it handle the outer surfaces of "1_surface1.png" as a single one instead of 2?<br>
<br>3) If there does be a sharp edge identification step, I am curious about its logic. Is there a reference toward a paper maybe presenting it?<br><br>4) In the paper, a benchmark is presented between different remeshing methods.<br>

The 2 first results present the performances of the surface remeshing based on the harmonic mapping algorithm implemented in GMSH, but using two different meshing algorithms: "meshadapt" and "del2d": is it possible to the user to chose which one to use with the current implementation?<br>

Currently which one is used?<br><br>5) Finally, I tried on a test case enclosed in teh 7zip archive: starting from "testUnconnectedTria.stl", and I would have liked to obtain a quad mesh similar to "testConnectedTriaQuad.bdf".<br>

Using GMSH 2.4.2, and the script "surfaceRemeshing.geo", I however obtain the error message below.<br>Do you know what is wrong?<br>~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~<br>Info    : Running '/opt/aster/outils/gmsh surfaceRemeshing.geo -2'<br>

Info    : Started on Thu Sep 16 22:10:08 2010<br>Info    : Reading 'surfaceRemeshing.geo'<br>Info    : Reading 'testUnconnectedTria.stl'<br>Info    : 29092 facets in solid 0<br>Info    : Read 'testUnconnectedTria.stl'<br>

Error   : 'surfaceRemeshing.geo', line 4 : syntax error (;)<br>Info    : Read 'surfaceRemeshing.geo'<br>Info    : Meshing 1D...<br>Info    : Mesh 1D complete (0 s)<br>Info    : Meshing 2D...<br>Info    : Mesh 2D complete (0 s)<br>

Info    : 15320 vertices 35204 elements<br>Info    : Writing 'surfaceRemeshing.msh'<br>Info    : Wrote 'surfaceRemeshing.msh'<br>Info    : Stopped on Thu Sep 16 22:11:33 2010<br>~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~<br>

<br>6) Also, do you know how I could ensure a quad mesh similar to the mesh: "testConnectedTriaQuad.bdf"?<br><br>I thank you in advance for all your help and advices which will be very appreciated.<br><br>Best regards,<br>

<br><br>Pierre<br><br><br><br><br><br><br><br><div class="gmail_quote">2010/8/17 Marchandise Emilie <span dir="ltr"><<a href="mailto:emilie.marchandise@uclouvain.be" target="_blank">emilie.marchandise@uclouvain.be</a>></span><br>

<blockquote class="gmail_quote" style="margin: 0pt 0pt 0pt 0.8ex; border-left: 1px solid rgb(204, 204, 204); padding-left: 1ex;">
Hi,<br>
<br>
As I said in a previous post, I encounter problem with surface remeshing<br>
considering the following *.geo and *.msh files.<br>
Could someone help me on this topic?<br>
I tried with a quad mesh saved in a *.msh format, and the message console<br>
tells me there is a syntax error at line 6 of my *.geo file containing the<br>
remeshing commands.<br>
The *.geo file is actually very similar to the one given as exemple in the<br>
wiki.<br>
<br>
<br>
>> I didn't find the .geo file in the mailing list (only the msh files). I can try to see what's wrong if you send it by e-mail.<br>
The reparametrization has been written for trinagular meshes. So, if you have a quad mesh, you should first transform it into a triangular mesh.<br>
<br>
Also, would someone give us an advice about the difference between Harmonic<br>
or Conformal mapping in terms maybe of quality and speed? Or in which kind<br>
of cases one is better suited than the other.<br>
<br>
>> Harmonic map parametrize the surface onto a unit disk while conformal map is an open boundray  parametrization (which improves the quality of the mesh near the boundary). It is however a bit slower and requires gmsh compiled with petsc and slepc.<br>



<br>
I tried with a quad mesh saved in a *.msh format, and the message console<br>
tells me there is a syntax error at line 6 of my *.geo file containing the<br>
remeshing commands.<br>
<br>
>> You should try with triangle mesh<br>
<br>
Finally, I have a last question: if I understand well, the logic for<br>
remeshing is to reparametrize a surface before recreating a mesh. Does it<br>
mean that the 'new' mesh is realized on a kind of 'B-spline surfaces' (which<br>
would be the reparametrized surface)?<br>
If yes, is it possible to save this re-created surface in brep or step?<br>
<br>
>>  No the parametrization is not a kind of B-spline surface but a planar surface (the surface has been flattened during the parametrization process)<br>
If you remesh with the option -saveall, you will see a bunch of files that show you the mapping u(x))<br>
see for example UVX-idSurf.pos that shows the X coordinates of the initial surface on the parametric plane UV.<br>
<br>
I thank you in advance for your advices and help.<br>
<br>
Bests,<br>
<br>
Pierre<br>
<br>
<br>
<br>
---------------------------------------------------<br>
Prof. Emilie Marchandise<br>
Université catholique de Louvain (UCL)<br>
Institute of Mechanics, Materials and Civil Engineering (iMMC)<br>
Center for Systems Engineering and Applied Mechanics (CESAME)<br>
Euler building, 4 avenue G. Lemaitre,<br>
B-1348 Louvain-la-neuve, BELGIUM<br>
Tel: 010/47.23.53<br>
e-mail: <a href="mailto:emilie.marchandise@uclouvain.be" target="_blank">emilie.marchandise@uclouvain.be</a><br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
</blockquote></div><br>