Hi guys, is it possible to build a edge mesh on FreeFem++?

For example, given a triangle, I wish to partition each edge in 2 and build a FEM space on each piece.

Hi guys, is it possible to build a edge mesh on FreeFem++?

For example, given a triangle, I wish to partition each edge in 2 and build a FEM space on each piece.

buildmeshL and then use periodic boundary conditions ?

I’ve tried like this:

```
load "msh3"
mesh Th = square(3,3);
meshL ThEdge = extract(Th);
```

but I can’t quite figure it out how to split ThEdge.

The extract function doesn’t extract the internal edges, is that right?

This came up with 12 dofs although I’m not sure if there is just someway to

save a meshL. Curious though if you wanted an actual closed rung.

That is, what would eigenvalues look like for some problem? Derivatives

zero at ends or merely solution same at both ends?

lari.edp (215 Bytes)

```
-- FreeFem++ v4.12 (Sat 20 May 2023 07:29:03 PM EDT - git no git)
file : lari.edp
Load: lg_fem lg_mesh lg_mesh3 eigenvalue
1 : load "msh3"
2 : mesh Th = square(3,3);
3 : meshL ThEdge = extract(Th);
4 : fespace Vh(ThEdge,P1);
5 : Vh vx=x;
6 : Vh vy=y;
7 : Vh vz=z;
8 : for(int i=0; i<vx[].n; ++i)
9 : {
10 : cout<<i<<" "<<vx[][i]<<" "<<vy[][i]<<" "<<vz[][i]<<endl;
11 :
12 : }
13 : cout.flush;
14 : sizestack + 1024 =1816 ( 792 )
-- Square mesh : nb vertices =16 , nb triangles = 18 , nb boundary edges 12
empty list label, extract all boundaries
-- FESpace: Nb of Nodes 12 Nb of DoF 12
0 0 0 0
1 0.333333 0 0
2 0.666667 0 0
3 1 0 0
4 1 0.333333 0
5 1 0.666667 0
6 1 1 0
7 0 1 0
8 0.333333 1 0
9 0.666667 1 0
10 0 0.333333 0
11 0 0.666667 0
times: compile 0.003532s, execution 0.000165s, mpirank:0
CodeAlloc : nb ptr 3839, size :513424 mpirank: 0
Ok: Normal End
```

But when it creates the mshL, it just extracts the square and not the inner boundary.

Is there some way to get all of then?

Probably mesh Th=square(1,1) would work.

I’ve tried with square(1,1) but no luck. I can only get external boundaries

```
load "msh3"
mesh Th = square(1,1);
meshL ThEdge = extract(Th);
ThEdge = trunc(ThEdge, 1, split=1);
plot(ThEdge, wait=1);
fespace Eh(ThEdge, P0);
Eh xx = x, yy=y;
int nDoFl0 = Eh.ndof;
int nDoFl0K = Eh.ndofK;
cout << nDoFl0 << "\n";
```

What problem are you trying to solve? I guess you could start with something

close and modify the matrix by hand.

I’m trying to implement a method to solve the Darcy equation

I’m gonna solve the Darcy equation on each element of the mesh and then assemble to the global mesh, but I need to build the location Matrix

For that, I need the global DoF numbers, so I can find the right places to assemble the values

The example you posted earlier gets the dof coords and for boundaries

in the 2d case you just used int1d over the labelled boundary. There

are a lot of examples for that.

Yes, but in thar example I can only get 4 boundaries, the ones that form the square

I need all five of them, the 4 boundaries from the square and the middle one

I guess you could add something like this to the mesh

border ba(t=0, 1){x=t; y=t; label=1;}

and then every element edge is a boundary.

or there is the “intAllEdge” if that helps.