Bloch boundary conditions

General Discussion
1

Is there a way to implement Bloch (i.e. “twisted” periodic) boundary conditions in FreeFEM?

2

yes ,

on square

mesh Th=square(10,10);
fespace Vh(Th,P1,periodic=[[4,y],[2,1-y]]);

…