Hi Daniel, Thank you for your help! I just tried a code snippet (modified step-49) below and it does not work, either when I set the dim=2 or dim=3.
```
int main ()
{
Triangulation<2> triangulation;
GridIn<2> 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);
}
```
The complete code is attached at the end of this email. If you don't mind,
could you help provide any guesses in what potentially went wrong?
Shawn
On Mon, Sep 17, 2018 at 12:07 PM Daniel Arndt <
[email protected]> wrote:
> Shawn,
>
> If I understand correctly you want to read a 2D mesh embedded in a
> three-dimensional space.
> Which mesh format are you interested in?
> At least, GridIn::read_ucd and GridIn::read_msh should be able to do that.
>
> If you experience problems, we like to here about that as well.
>
> Best,
> Daniel
>
> --
> 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
> [email protected].
> For more options, visit https://groups.google.com/d/optout.
>
--
Yuxiang "Shawn" Wang, PhD
[email protected]
+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 [email protected].
For more options, visit https://groups.google.com/d/optout.
##
# 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()
#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> triangulation;
GridIn<2> 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
