[N,psi_0]=gf_spmat_get(C,'dirichlet_nullspace',R);
psi_0 = psi_0'; %psi_0 generated by the gf_spmat command is a row. Convert to column. HH0 = N' * H0 * N; SS = N' * S * N;
I find that : 1. HH0 has no imag entries! Hence, psi has no imag entries. The eigvalues are wrong. 2. N is of size nbdof X (nbdof-2) though I have imposed only one condition! If I had imposed the condition \psi(0) = (some real) * \psi(1), then N is of size nbdof X nbdof -1.
Is it possible to impose Bloch periodic conditions in getfem?