Thanks! It works now.
I was only using Mesh.CharacteristicLengthFromPoints = 0; but I was missing that other option. In addition to that I had to delete the lines Mesh.CharacteristicLengthExtendFromBoundary = 1; Mesh.CharacteristicLengthFromPoints = 1; from the .opt file because gmsh reads first the .geo and then the .opt Sincerely, Héctor Espinoza ________________________________ De: Christophe Geuzaine <[email protected]> Enviado: viernes, 22 de junio de 2018 1:39:06 Para: Hector Gabriel Espinoza Román Cc: [email protected] Asunto: Re: [Gmsh] MathEval field in a circle with Open Cascade On 22 Jun 2018, at 00:55, Hector Gabriel Espinoza Román <[email protected]<mailto:[email protected]>> wrote: Hello: This simple script should draw a circle (2D surface) and when meshed it should be refined on the boundary, but I am not getting that. I have done the same without Open Cascade and it works. cx = 0.0 ; cy = 0.0 ; cz = 0.0 ; D = 1.0 ; R = D/2; SetFactory("OpenCASCADE"); Disk(1) = {cx,cy,cz,R} ; // creates a circle surface Field[1] = MathEval ; //Field[1].F = "0.15 - 0.13*((x^2)/(0.5^2))" ; // works 2D, 617 elements: 0.02 to 0.15 Field[1].F = "0.15 - 0.13*(x^2 + y^2)/(0.5^2)" ; // doesnt work 2D Background Field = 1 ; What I am doing wrong? Add Mesh.CharacteristicLengthExtendFromBoundary = 0; Christophe Sincerely, Héctor _______________________________________________ gmsh mailing list [email protected]<mailto:[email protected]> http://onelab.info/mailman/listinfo/gmsh — Prof. Christophe Geuzaine University of Liege, Electrical Engineering and Computer Science http://www.montefiore.ulg.ac.be/~geuzaine Free software: http://gmsh.info | http://getdp.info | http://onelab.info
_______________________________________________ gmsh mailing list [email protected] http://onelab.info/mailman/listinfo/gmsh
