Thanks for the tip > On Jul 13, 2018, at 10:33 PM, Drew Davidson <davidson...@gmail.com> wrote: > > Hello, > > It seems possible to generate a parameterized 1D mesh using gmsh within a > python script, and then generate the corresponding 1D FiPy mesh in that same > python script: > > 1) Use pygmsh in a python script to construct a parameterized 1D (gmsh) mesh. > > 2) From that 1D mesh, construct a list of dx's > > 3) Follow https://www.mail-archive.com/fipy@nist.gov/msg01650.html and create > a 1D FiPy mesh from a list of dx's. > > It was pure luck to have found the mailing list message in (3). > > I did not see how to get the Physical Groups data from the pygmsh mesh into > the FiPy mesh. > > Thanks > > > > On Thu, Jul 12, 2018 at 3:26 PM Guyer, Jonathan E. Dr. (Fed) > <jonathan.gu...@nist.gov> wrote: > No, FiPy doesn't presently have any ability to read 1D meshes. > > I do have some experimental code that supports this, but it needs a lot of > work. I have no prognosis for when this will be made public. > > - Jon > > > On Jul 12, 2018, at 1:21 PM, Drew Davidson <davidson...@gmail.com> wrote: > > > > Hello, > > > > Can Gmsh2D read a 1D gmsh mesh? > > > > I have tried a similar script to the FiPy diffusion example circle.py: > > > > SCRIPT BEGINS > > import fipy as fp > > from IPython import embed > > import sys > > > > #I created this gmsh .geo file in gmsh GUI, then copied and pasted here; I > > cannot see a syntax error if it exists??? > > > > geo=""" > > L_Sample=%(L_Sample)g; > > L_Substrate=%(L_Substrate)g; > > cellSize=%(cellSize)g; > > Point(1) = {0, 0, 0, cellSize}; > > Point(2) = {L_Substrate, 0, 0, cellSize}; > > Point(3) = {L_Substrate+L_Sample, 0, 0, cellSize}; > > Line(1) = {1, 2}; > > Line(2) = {2, 3}; > > Physical Point("SubstrateEnd") = {1}; > > Physical Point("SampleEnd") = {3}; > > Physical Point("ThermalContact") = {2}; > > Physical Line("Substrate") = {1}; > > Physical Line("Sample") = {2}; > > """ % locals() > > > > #I am not having success with: > > mesh = fp.Gmsh2D(geo,coordDimensions=1) #TODO throwing error about there > > being no cells; is it a problem that mesh is 1D and Gmsh2D is called 2D? > > mesh = fp.Gmsh2D(geo) #TODO throwing error about there being no cells > > > > #then I went into gmsh GUI and saved to .msh, but still no luck: > > mesh = fp.Gmsh2D('OneDimensionalSampleAndSubstrate.msh') #TODO throwing > > error about there being no cells > > mesh = fp.Gmsh2D('OneDimensionalSampleAndSubstrate.msh',coordDimensions=1) > > #TODO throwing error about there being no cells > > SCRIPT ENDS > > > > > > If coordDimensions=1 is unsupported or meaningless, it would be great if > > FiPy printed a message saying so. > > > > I believe my FiPy installation is OK and was downloaded from github > > 6/30/2018 (fipy-develop). I can run circle.py, although I have never gotten > > it to plot the mesh. I can run anisotropy.py from diffusivity examples. > > Nearly all fipy tests pass (setup.py). Can you please check if Gmsh2D is > > capable of reading a 1D gmsh mesh? Or is there another FiPy class which can > > do it, and how? > > > > gmsh –version returns 3.0.6. Installation is system-wide via download > > method in Ubuntu 16.04. > > > > If I cannot read in a parameterized 1D gmsh mesh using FiPy functionality, > > I am at a loss as to how to have a single Python script that allows me to > > set cellSize, L_Sample, and L_Substrate, produce a mesh, and efficiently > > conduct a parameter study with FiPy. Gmsh is desirable because I want to go > > to a more complicated nonuniform/graded mesh, and would rather not write > > code to make a mesh myself. > > > > Thanks. > > > > > > _______________________________________________ > > fipy mailing list > > fipy@nist.gov > > http://www.ctcms.nist.gov/fipy > > [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ] > > > _______________________________________________ > fipy mailing list > fipy@nist.gov > http://www.ctcms.nist.gov/fipy > [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ] > _______________________________________________ > fipy mailing list > fipy@nist.gov > http://www.ctcms.nist.gov/fipy > [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]
_______________________________________________ fipy mailing list fipy@nist.gov http://www.ctcms.nist.gov/fipy [ NIST internal ONLY: https://email.nist.gov/mailman/listinfo/fipy ]