gridap / gridap/GridapEmbedded.jl

Bounds Error when no cut cells are present

Open
#37 2 comments 0 reactions 0 assignees View on GitHub
Dominant language
Julia
Stars
50
Forks
16
PR merge metrics
No merged PRs in 30d

Description

@fverdugo @santiagobadia
If we consider the following code
```
using Gridap
using GridapEmbedded
using Gridap.ReferenceFEs

u(x) = x[1] + x[2]
f(x) = 0.0

bgdomain = (0,1,0,1)
n = 10
partition = (n,n)
bgmodel = CartesianDiscreteModel(bgdomain,partition)

geo = square(L=0.8,x0=Point(0.5,0.5))
cutgeo = cut(bgmodel,geo)
model = DiscreteModel(cutgeo)

strategy = AggregateCutCellsByThreshold(1.0)
aggregates = aggregate(strategy,cutgeo)

T = Float64
order = 1
reffe = LagrangianRefFE(T,QUAD,order)
V0 = TestFESpace(reffe=reffe, model=model,dof_space=:physical)
V = AgFEMSpace(V0,aggregates)
U = TrialFESpace(V)
```
Here, all cells are well-posed (no cut cells).
![image](https://user-images.githubusercontent.com/53897026/100200780-356c1980-2f53-11eb-89f0-640cfae1edf8.png)

We can see an error **BoundsError: attempt to access 81-element Array{Int32,1} at index [0]** in the line
`V = AgFEMSpace(V0,aggregates)`.

Contributor guide

No contributing guide indexed for this repository

Assessment

This issue has not been assessed yet.

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.