Generate error message if attempting to use VTK writer with OutputSolnPoints.

Author: baagaard-usgs

libsrc/pylith/meshio/OutputSolnPoints.cc

@@ -13,6 +13,7 @@
 #include "pylith/meshio/OutputSolnPoints.hh" // implementation of class methods
 
 #include "pylith/meshio/DataWriter.hh" // USES DataWriter
+#include "pylith/meshio/DataWriterVTK.hh" // USES DataWriterVTK
 #include "pylith/meshio/MeshBuilder.hh" // USES MeshBuilder
 
 #include "pylith/topology/Mesh.hh" // USES Mesh
@@ -106,6 +107,10 @@ pylith::meshio::OutputSolnPoints::_writeSolnStep(const PylithReal t,
     PYLITH_METHOD_BEGIN;
     PYLITH_COMPONENT_DEBUG("_writeSolnStep(t="<<t<<", tindex="<<tindex<<", solution="<<solution.getLabel()<<")");
 
+    if (dynamic_cast<pylith::meshio::DataWriterVTK*>(_writer)) {
+        PYLITH_JOURNAL_LOGICERROR("PETSc VTK writer using the VTU format does not support output at points. Use the default DataWriterHDF5 writer.");
+    } // if
+
     if (!_interpolator) {
         _setupInterpolator(solution);
     } // if

pylith/meshio/OutputSolnPoints.py

@@ -83,6 +83,18 @@ def preinitialize(self, problem):
 
     # PRIVATE METHODS ////////////////////////////////////////////////////
 
+    def _validate(self, context):
+        from .DataWriterVTK import DataWriterVTK
+        if isinstance(self.inventory.writer, DataWriterVTK):
+            trait = self.inventory.getTrait("writer")
+            self._validationError(context, trait, "PETSc VTK writer using the VTU format does not support output at points. Use the default DataWriterHDF5 writer.")
+
+    def _validationError(self, context, trait, msg):
+        from pythia.pyre.inventory.Item import Item
+        error = ValueError(msg)
+        descriptor = self.getTraitDescriptor(trait.name)
+        context.error(error, items=[Item(trait, descriptor)])
+
     def _createModuleObj(self):
         """Create handle to C++ object.
         """