Taking the Neumann boundary conditions into account

 

Computation of the right-hand side :

The vector which must be assembled is defined by :

with :

- i are shape functions defined on the reference interval
- is the number of nodes per edge
- n1 and n2 are the start and end nodes of the edge
- g2 is the value of the boundary condition on the edge

 


 

Assembling process :

For each border domain of type 2

For each edge of the domain

Computation of f(e2)

For each node of the edge, which has i as global number and k as local number