Skip to content

support visualizing higher order fields with paraview  #45

@cwsmith

Description

@cwsmith

See

// FIXME support 2d and 3d and fields with order>1
template <typename Field> void writeVtk(Field &field) {
using FieldDataType = typename decltype(field.vtxField)::BaseType;
// HACK assumes there is a vertex field.. in the Field Mixin object
auto field_view = field.vtxField.serialize();
Omega_h::Write<FieldDataType> field_write(field_view);
mesh.add_tag(0, "field", 1, Omega_h::read(field_write));
Omega_h::vtk::write_parallel("foo.vtk", &mesh, mesh.dim());
}

Metadata

Metadata

Assignees

No one assigned

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions