thanks! when generating the header file ex52.h what should the input arguments be to the PetscGenerateFEMQuadrature.py file ?
SCRGP2$ $PETSC_DIR/bin/pythonscripts/PetscGenerateFEMQuadrature.py 3 1 8 1 laplacian ex52.h [{(-1.0, -1.0, -1.0): [(1.0, ())]}, {(1.0, -1.0, -1.0): [(1.0, ())]}, {(-1.0, 1.0, -1.0): [(1.0, ())]}, {(-1.0, -1.0, 1.0): [(1.0, ())]}] {0: {0: [0], 1: [1], 2: [2], 3: [3]}, 1: {0: [], 1: [], 2: [], 3: [], 4: [], 5: []}, 2: {0: [], 1: [], 2: [], 3: []}, 3: {0: []}} Perm: [0, 1, 2, 3] Creating /home/fuentes/snestutorials/ex52.h Traceback (most recent call last): File "/opt/apps/PETSC/petsc-dev/bin/pythonscripts/PetscGenerateFEMQuadrature.py", line 33, in <module> generator.run(elements, numBlocks, operator, filename) File "/opt/apps/PETSC/petsc-dev/config/PETSc/FEM.py", line 967, in run self.outputElementSource(self.getElementSource(elements, numBlocks, operator, sourceType = 'GPU'), os.path.splitext(filename)[0]+'_gpu'+os.path.splitext(filename)[1]) File "/opt/apps/PETSC/petsc-dev/config/PETSc/FEM.py", line 926, in getElementSource defns.extend(self.getComputationTypes(element, n)) File "/opt/apps/PETSC/petsc-dev/config/PETSc/FEM.py", line 440, in getComputationTypes vecType.type = self.Cxx.typeMap['float'+str(dim)] KeyError: 'float3' On Tue, Mar 27, 2012 at 2:10 PM, Blaise Bourdin <bourdin at lsu.edu> wrote: > > On Mar 27, 2012, at 1:23 PM, Matthew Knepley wrote: > > On Tue, Mar 27, 2012 at 12:58 PM, David Fuentes <fuentesdt at gmail.com>wrote: > >> Hi, >> >> I had a question about the status of example 52. >> >> >> http://petsc.cs.iit.edu/petsc/petsc-dev/file/a8e2f2c19319/src/snes/examples/tutorials/ex52.c >> >> http://petsc.cs.iit.edu/petsc/petsc-dev/file/a8e2f2c19319/src/snes/examples/tutorials/ex52_integrateElement.cu >> >> >> Can this example be used with a DM object created from an unstructured >> exodusII mesh, DMMeshCreateExodus, And the FEM assembly done on GPU ? >> > > 1) I have pushed many more tests for it now. They can be run using the > Python build system > > ./config/builder2.py check src/snes/examples/tutorials/ex52.c > > in fact, you can build any set of files this way. > > 2) The Exodus creation has to be converted to DMComplex from DMMesh. That > should not take me very long. Blaise maintains that > so maybe there will be help :) You will just replace > DMComplexCreateBoxMesh() with DMComplexCreateExodus(). If you request > it, I will bump it up the list. > > > DMMeshCreateExodusNG is much more flexible than DMMeshCreateExodus in that > it can read meshes with multiple element types and should have a much lower > memory footprint. The code should be fairly easy to read. you can email me > directly if you have specific questions. I had looked at creating a > DMComplex and it did not look too difficult, as long as interpolation is > not needed. I have plans to write DMComplexCreateExodus, but haven't had > time too so far. Updating the Vec viewers and readers may be a bit more > involved. In perfect world, one would write an EXODUS viewer following the > lines of the VTK and HDF5 ones. > > Blaise > > > > Let me know if you can run the tests. > > Thanks > > Matt > > >> Thanks, >> David >> > -- > What most experimenters take for granted before they begin their > experiments is infinitely more interesting than any results to which their > experiments lead. > -- Norbert Wiener > > > -- > Department of Mathematics and Center for Computation & Technology > Louisiana State University, Baton Rouge, LA 70803, USA > Tel. +1 (225) 578 1612, Fax +1 (225) 578 4276 > http://www.math.lsu.edu/~bourdin > > > > > > > > -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20120327/105d5e4b/attachment.htm>