Prof. Bangerth, Thank you for the prompt response! Yes please find attached the minimal program. The main body is pretty short (see below), modified from step-49. Please let me know if you cannot reproduce it!
```
int main ()
{
Triangulation<2, 3> triangulation;
GridIn<2, 3> gridin;
gridin.attach_triangulation(triangulation);
std::ifstream f("mesh.inp");
gridin.read_abaqus(f);
GridOut gridout;
std::ofstream fout("mesh.vtu");
gridout.write_vtu(triangulation, fout);
}
```
Shawn
On Tue, Sep 18, 2018 at 4:28 PM Wolfgang Bangerth <bange...@colostate.edu>
wrote:
> On 09/18/2018 03:18 PM, Yuxiang Wang wrote:
> >
> > Thank you for the pointer! I changed my code as you suggested and it
> finished
> > running.
>
> Great!
>
>
> > However, when I opened the output file (mesh.vtu), the 3D mesh was
> projected
> > to a 2D plane of z=0, and lost its 3D features. Was that expected
> behavior?
>
> No. Are you saying that if you read a mesh, then immediately output it
> again,
> that the z-coordinate is lost? If that's correct, then that's a bug and it
> would be nice to have a small program and input file that demonstrates
> this.
>
> Best
> W.
>
> --
> ------------------------------------------------------------------------
> Wolfgang Bangerth email: bange...@colostate.edu
> www: http://www.math.colostate.edu/~bangerth/
>
> --
> The deal.II project is located at http://www.dealii.org/
> For mailing list/forum options, see
> https://groups.google.com/d/forum/dealii?hl=en
> ---
> You received this message because you are subscribed to a topic in the
> Google Groups "deal.II User Group" group.
> To unsubscribe from this topic, visit
> https://groups.google.com/d/topic/dealii/Erk7pxwTEj4/unsubscribe.
> To unsubscribe from this group and all its topics, send an email to
> dealii+unsubscr...@googlegroups.com.
> For more options, visit https://groups.google.com/d/optout.
>
--
Yuxiang "Shawn" Wang, PhD
yw...@virginia.edu
+1 (434) 284-0836
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.
#include <deal.II/grid/tria.h>
#include <deal.II/grid/tria_accessor.h>
#include <deal.II/grid/tria_iterator.h>
#include <deal.II/grid/grid_generator.h>
#include <deal.II/grid/grid_tools.h>
#include <deal.II/grid/tria_boundary_lib.h>
#include <deal.II/grid/grid_out.h>
#include <deal.II/grid/grid_in.h>
#include <fstream>
#include <map>
using namespace dealii;
int main ()
{
Triangulation<2, 3> triangulation;
GridIn<2, 3> gridin;
gridin.attach_triangulation(triangulation);
std::ifstream f("mesh.inp");
gridin.read_abaqus(f);
GridOut gridout;
std::ofstream fout("mesh.vtu");
gridout.write_vtu(triangulation, fout);
}
mesh.inp
Description: Binary data
##
# CMake script for the step-1 tutorial program:
##
# Set the name of the project and target:
SET(TARGET "step-49")
# Declare all source files the target consists of:
SET(TARGET_SRC
${TARGET}.cc
# You can specify additional files here!
)
# Usually, you will not need to modify anything beyond this point...
CMAKE_MINIMUM_REQUIRED(VERSION 2.8.8)
FIND_PACKAGE(deal.II 8.0 QUIET
HINTS ${deal.II_DIR} ${DEAL_II_DIR} ../ ../../ $ENV{DEAL_II_DIR}
)
IF(NOT ${deal.II_FOUND})
MESSAGE(FATAL_ERROR "\n"
"*** Could not locate deal.II. ***\n\n"
"You may want to either pass a flag -DDEAL_II_DIR=/path/to/deal.II to
cmake\n"
"or set an environment variable \"DEAL_II_DIR\" that contains this path."
)
ENDIF()
DEAL_II_INITIALIZE_CACHED_VARIABLES()
PROJECT(${TARGET})
DEAL_II_INVOKE_AUTOPILOT()
mesh.vtu
Description: Binary data
