Hi Weston,

Thanks for the report and the test data. I'll see that the fix gets merged in 
for both HEX9 and HEX21.

        David

> On Feb 12, 2018, at 10:31 AM, Weston Ortiz <wes...@wortiz.com> wrote:
> 
> Using
> 
>   else if ((elemType.substr(0,3) == "HEX") && (binfo.BdsPerEntry[0] == 9))
>     { binfo.CellType=VTK_HEXAHEDRON;          binfo.PointsPerCell = 9; }
> 
> Appears to render correctly (see attached plot)
> 
> I changed the Points per cell to 9 as well, with it as 8 the output was the 
> same as before.
> 
> I'm guessing this is because exodus stores arrays by element (e.g. [<values 
> for element 1, nodes 1 to 9>, <values for element 2, nodes 1 to 9>, ... 
> etc.]) and PointsPerCell is used to index into these
> 
> If so your hex21 might also need to be changed, cubit/trelis doesn't support 
> hex21 so I couldn't check.
> 
> Thanks,
> 
> Weston
> 
> On Fri, Feb 9, 2018 at 4:16 PM, Weston Ortiz <wes...@wortiz.com> wrote:
> > Have you built ParaView yourself or do you install it from packages.
> 
> Usually just from the paraview tarballs, i can try adding the new change code 
> next week though.
> 
> > It would help to have a small example file for us to test against.
> 
> Attached is a small example file (lid driven cavity, ldc.exoII) and an 
> example output image of a slice with velocity vectors in ensight
> 
> Thanks,
> 
> Weston
> 
> On Fri, Feb 9, 2018 at 3:17 PM, David Thompson <david.thomp...@kitware.com> 
> wrote:
> Hi Weston,
> 
> > Is there any way I could render HEX9 exodusii data as as HEX8 data in 
> > paraview,  (ignoring the 9th node) or work with it directly?
> 
> Not without modifying the source. Have you built ParaView yourself or do you 
> install it from packages.
> 
> >
> > HEX9 numbering: 
> > https://github.com/gsjaardema/seacas/blob/master/docs/topology/hex09.png
> >
> > It appears that Paraview reads the file as HEX8 but is reading the center 
> > node as part of one of the hexahedron's edges (see attached image)
> 
> It would help to have a small example file for us to test against. If you 
> want to experiment, you can edit VTK/IO/Exodus/vtkExodusIIReader.cxx around 
> line 3237 (inside vtkExodusIIReaderPrivate::DetermineVtkCellType) and add
> 
>   else if ((elemType.substr(0,3) == "HEX") &&      (binfo.BdsPerEntry[0] == 
> 9))
>     { binfo.CellType=VTK_HEXAHEDRON;               binfo.PointsPerCell = 8; }
> 
> You can see that is what we do for 21-node hexes currently.
> 
>         David
> 
> 
> > <simple_cube.png>_______________________________________________
> > Powered by www.kitware.com
> >
> > Visit other Kitware open-source projects at 
> > http://www.kitware.com/opensource/opensource.html
> >
> > Please keep messages on-topic and check the ParaView Wiki at: 
> > http://paraview.org/Wiki/ParaView
> >
> > Search the list archives at: http://markmail.org/search/?q=ParaView
> >
> > Follow this link to subscribe/unsubscribe:
> > https://paraview.org/mailman/listinfo/paraview
> 
> 
> 
> <hex9plot.png>

_______________________________________________
Powered by www.kitware.com

Visit other Kitware open-source projects at 
http://www.kitware.com/opensource/opensource.html

Please keep messages on-topic and check the ParaView Wiki at: 
http://paraview.org/Wiki/ParaView

Search the list archives at: http://markmail.org/search/?q=ParaView

Follow this link to subscribe/unsubscribe:
https://public.kitware.com/mailman/listinfo/paraview

Reply via email to