Dear all,
I'm trying to use the class GridIn to import a mesh into the program. My
program compiled fine but had a runtime error as follows:
%./step-1
--------------------------------------------------------
An error occurred in line <1794> of file <source/grid/grid_in.cc> in
function
void dealii::GridIn<dim, spacedim>::read_tecplot(std::istream&) [with
int dim = 3, int spacedim = 3]
The violated condition was:
false
The name and call sequence of the exception was:
ExcNotImplemented()
Additional Information:
(none)
Stacktrace:
-----------
#0 ./step-1: dealii::GridIn<3, 3>::read_tecplot(std::istream&)
#1 ./step-1: main
--------------------------------------------------------
Aborted
My code is pasted below (without the includes, since everything compiled
fine. This is a smaller problem just to make sure I know how to import
things. My real problem is a 3D unstructured grid with about a million
nodes...):
using namespace dealii;
int main ()
{
Triangulation<3, 3> tria;
std::ifstream st("tet");
GridIn<3, 3> gridIn;
gridIn.attach_triangulation(tria);
gridIn.read_tecplot(st);
st.close();
std::ofstream out ("grid-1.eps");
GridOut grid_out;
grid_out.write_eps (tria, out);
return(0);
}
And the input file "tet" is just:
title = "tet!"
variables = "x", "y", "z"
zone n=4, e=1, et=TETRAHEDRON
1.0 1.0 1.0
1.0 -1.0 -1.0
-1.0 1.0 -1.0
-1.0 -1.0 1.0
2 4 1 3
(I also tried by replacing "et=TETRAHEDRON" with "DATAPACKING=POINT" and
"DATAPACKING=POINT, ZONTYPE=TETRAHEDRON", which all resulted in the same
error.)
Did I do something wrong with the GridIn class? I also noticed that the
GridIn class is able to read in many different types of data formats, but am
not sure the "default" format means or what format "read_ucd" takes
exactly.... I couldn't seem to find any example code that deals with
importation of meshes. If they actually exist on the site, could you guys
direct me to it? Otherwise, could you show me a simple example how this is
done?
Thank you very much!
Karen
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
