Hello,
I'm using the python interface of GetFEM. My goal is to compute the laplace eigenvalues on a region of a spherical mesh, as follows (suppose the mesh, integration method IM and meshfem V is already declared):
# Assemble the stiffness and mass matrices
md = gf.Model("real")
md.add_fem_variable("u", V)
Kg = gf.asm_generic(IM, 2, "Grad_Test2_u.Grad_Test_u", region, md)
Mg = gf.asm_generic(IM, 2, "Test2_u*Test_u", region, md)
# Convert to scipy sparse
K = getfem_to_scipy_sparse(Kg)
M = getfem_to_scipy_sparse(Mg)
# Compute the eigenvalues
eVal, eVec = scipy.linalg.sparse.eigsh(K, k+max_multiplicity+1, M, sigma=0, which='LM', v0=v0)
however, the size of the matrices K and M are the same as V.nbdof() whereas it should have less since the region isn't the whole mesh. Moreover, when I try to use the function "basic_dof_on_region", I get the following error :
NameError: name 'basic_dof_on_region' is not defined
I saw the things about the extension matrices but this is not clear to me how to construct it from the label of the region. So my question is what is the best way to compute eigenvalues on a region ?
Thank you for your help and have a good day.
Eloi MARTINET.