Now, I am wondering how algebraic method constrains hanging nodes. Whether does it set the values at hanging nodes to zero? I output the matrix assembled in libmesh using PetSC function, I find lots of zero values in matrix. Because the matrix in Petsc is stored using sparse compressed strage format, zero values should not appear. Could you give me some hints?

thanks a lot.

Regards,

Yujie