Dear John, The problem is solved. As you mentioned, the subdomain id was not well assigned. Technically, the problem was related to not typecasting PetscScalar to double in my code, not directly involving with libMesh.
Thanks! With best regards, Pilhwa _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users