[petsc-users] ex52_integrateElement.cu

[David Fuentes]

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>