Hello all,
for a task, one model contains shells (surface elements) and beams (line elements). The goal was to get a simple working model. Please see the picture:
If this is loaded with GridapGmsh an error message occurs - does anyone know how to use Gridap with two different dimensional entities, i.e. line and surface elements?
using Gridap
using GridapGmsh
model = GmshDiscreteModel("test.msh")
error message:
Info : Reading ‘Code/Meshes/Data/0_MSH/beam_7.msh’…
Info : 13 entities
Info : 422 nodes
Info : 814 elements
Info : Done reading ‘Code/Meshes/Data/0_MSH/beam_7.msh’
ERROR: BoundsError: attempt to access 1129-element Vector{Int64} at index [0]
Stacktrace:
[1] setindex!
@ ./array.jl:966 [inlined]
[2] setindex!
@ ./multidimensional.jl:674 [inlined]
[3] _apply_offset_for_faces!(face_to_label::Vector{Int64}, gface_to_entity::Vector{Int64}, gface_to_face::Vector{Int32}, offset::Int64)
@ GridapGmsh ~/.julia/packages/GridapGmsh/Q8ZwW/src/GmshDiscreteModels.jl:604
[4] _fill_dim_to_face_to_label!(node_to_vertex::UnitRange{Int64}, dim_to_face_to_label::Vector{Vector{Int64}}, grid_topology::Gridap.Geometry.UnstructuredGridTopology{2, 3, Float64, Gridap.Geometry.NonOriented}, dim_to_gface_to_nodes::Vector{Any}, dim_gface_to_entity::Vector{Any}, dim_to_offset::Vector{Int64})
@ GridapGmsh ~/.julia/packages/GridapGmsh/Q8ZwW/src/GmshDiscreteModels.jl:576
[5] _setup_face_labels(vertex_to_node::UnitRange{Int64}, node_to_vertex::UnitRange{Int64}, grid_topology::Gridap.Geometry.UnstructuredGridTopology{2, 3, Float64, Gridap.Geometry.NonOriented}, dim_to_gface_to_nodes::Vector{Any}, dim_gface_to_entity::Vector{Any}, dim_to_offset::Vector{Int64}, dim_to_group_to_entities::Vector{Vector{Vector{Int64}}}, dim_to_group_to_name::Vector{Vector{String}}, node_to_label::Vector{Int64})
@ GridapGmsh ~/.julia/packages/GridapGmsh/Q8ZwW/src/GmshDiscreteModels.jl:531
[6] _setup_labeling(gmsh::Module, grid::Gridap.Geometry.UnstructuredGrid{2, 3, Float64, Gridap.Geometry.NonOriented, Gridap.Arrays.LazyArray{FillArrays.Fill{typeof(Gridap.Fields.constant_field), 1, Tuple{Base.OneTo{Int64}}}, Gridap.Fields.ConstantField{VectorValue{3, Float64}}, 1, Tuple{Vector{VectorValue{3, Float64}}}}}, grid_topology::Gridap.Geometry.UnstructuredGridTopology{2, 3, Float64, Gridap.Geometry.NonOriented}, cell_to_entity::Vector{Int64}, vertex_to_node::UnitRange{Int64}, node_to_vertex::UnitRange{Int64})
@ GridapGmsh ~/.julia/packages/GridapGmsh/Q8ZwW/src/GmshDiscreteModels.jl:414
[7] GmshDiscreteModel(mshfile::String; renumber::Bool)
@ GridapGmsh ~/.julia/packages/GridapGmsh/Q8ZwW/src/GmshDiscreteModels.jl:24
[8] GmshDiscreteModel(mshfile::String)
@ GridapGmsh ~/.julia/packages/GridapGmsh/Q8ZwW/src/GmshDiscreteModels.jl:6
[9] top-level scope
@ ~/Schreibtisch/Javed/1_Projekte/1_K2H2/6_Tasks/1_FEM/6_Meshing/2_Meshes_Jl/Code/First_Steps/Gm_1st.jl: