<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Dear Guillaume,<div class=""><br class=""></div><div class="">Thanks for the quick answer, it is then an actual bug. IN my understanding, the current orientation should provide the orientation of the cut. I have already implemented the proposed solution, it is not so ideal but it leads to a temporary fix.</div><div class=""><br class=""></div><div class="">Best,</div><div class=""><br class=""></div><div class="">Frederic</div><div class=""><br class=""></div><div class=""><br class=""><div><blockquote type="cite" class=""><div class="">On Feb 13, 2017, at 01:40, DILASSER Guillaume <<a href="mailto:Guillaume.DILASSER@cea.fr" class="">Guillaume.DILASSER@cea.fr</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div class="WordSection1" style="page: WordSection1; font-family: Helvetica; font-size: 14px; font-style: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;"><div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;" class="">Dear Frederic,<o:p class=""></o:p></div><div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;" class=""><o:p class=""> </o:p></div><div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;" class=""><span lang="EN-US" class="">Your issue is caused by inconsistent orientation of cohomology cuts. If you mesh and generate cohomology cuts for your model within the Gmsh GUI, you will see that 3 of your cuts are oriented in one direction and the one corresponding to the first turn is in the opposite direction. Therefore the current flows in the opposite direction in the first turn. A quick fix is to set up your current constraint in the pro file on a per turn basis (put +I on each turn with the correct orientation and -I on the other) but this is not practical when the number of turns grows. Another solution is to post-process the msh file and reverse the orientation of each line element within the cuts that have bad orientation (this is presently what I do). There may be a way to handle the issue directly in Gmsh but unfortunately I was not able to find it.<o:p class=""></o:p></span></div><div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;" class=""><span lang="EN-US" class=""><o:p class=""> </o:p></span></div><div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;" class=""><span lang="EN-US" class="">Sincerely Yours,<o:p class=""></o:p></span></div><div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;" class=""><o:p class=""> </o:p></div><div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;" class=""><span style="font-family: Tahoma, sans-serif; color: rgb(0, 0, 153);" class="">Guillaume DILASSER</span><span style="font-size: 10pt; font-family: Tahoma, sans-serif; color: rgb(0, 0, 153);" class=""><o:p class=""></o:p></span></div><div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;" class=""><span style="font-size: 7.5pt; font-family: Arial, sans-serif; color: rgb(0, 0, 153);" class="">Doctorant SACM / LEAS<o:p class=""></o:p></span></div><div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;" class=""><span style="font-size: 7.5pt; font-family: Arial, sans-serif; color: rgb(0, 0, 153);" class="">CEA - Centre de Saclay - Bât.123 - PC 319c</span><span style="color: rgb(0, 0, 153);" class=""><o:p class=""></o:p></span></div><div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;" class=""><span style="font-size: 7.5pt; font-family: Arial, sans-serif; color: rgb(0, 0, 153);" class="">91191 Gif sur Yvette Cedex - France -</span><span style="color: rgb(0, 0, 153);" class=""><o:p class=""></o:p></span></div><div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;" class=""><span style="font-size: 7.5pt; font-family: Arial, sans-serif; color: rgb(0, 0, 153);" class=""><o:p class=""> </o:p></span></div><div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;" class=""><span style="font-size: 7.5pt; font-family: Arial, sans-serif; color: rgb(0, 0, 153);" class=""><a href="mailto:guillaume.dilasser@cea.fr" style="color: rgb(149, 79, 114); text-decoration: underline;" class=""><span style="color: blue;" class="">guillaume.dilasser@cea.fr</span></a><o:p class=""></o:p></span></div><div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;" class=""><o:p class=""> </o:p></div><div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;" class=""><span lang="EN-US" class=""><o:p class=""> </o:p></span></div><div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;" class=""><span class="">-----Message d'origine-----<br class="">De : getdp [<a href="mailto:getdp-bounces@ace20.montefiore.ulg.ac.be" class="">mailto:getdp-bounces@ace20.montefiore.ulg.ac.be</a>] De la part de Frederic Trillaud Pighi<br class="">Envoyé : samedi 11 février 2017 02:25<br class="">À : <a href="mailto:getdp@geuz.org" class="">getdp@geuz.org</a><br class="">Cc : Cesar Simón Lopez Monsalvo <<a href="mailto:cslopezmo@conacyt.mx" class="">cslopezmo@conacyt.mx</a>><br class="">Objet : [Getdp] Issue with cohomology solver</span></div><div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;" class=""><o:p class=""> </o:p></div><div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;" class="">Dear all,<o:p class=""></o:p></div><div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;" class=""><o:p class=""> </o:p></div><div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;" class="">I do not know if I am doing wrong here. I have an axisymmetric model of superconducting tapes in parallel carrying the same impressed current.<o:p class=""></o:p></div><div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;" class=""><o:p class=""> </o:p></div><div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;" class="">In the case of triangular mesh, all the tapes carry the current in the same direction leading to the expected total current integrated over all the tapes. However, when I use quadrangles (exact same code), the first tape does not have the same direction of current as the remaining tape. It seems that there is an issue with the association of the current with the cut. I may think that it is some kind of a bug resulting from the choice of quadrangle meshing.<o:p class=""></o:p></div><div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;" class=""><o:p class=""> </o:p></div><div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;" class="">I have joined the code based and resulting results.<o:p class=""></o:p></div><div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;" class=""><o:p class=""> </o:p></div><div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;" class="">Best,<o:p class=""></o:p></div><div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;" class=""><o:p class=""> </o:p></div><div style="margin: 0cm 0cm 0.0001pt; font-size: 11pt; font-family: Calibri, sans-serif;" class="">Frederic</div></div></div></blockquote></div><br class=""></div></body></html>