Hello, Pleas see attached following code, http://userweb.eng.gla.ac.uk/lukasz.kaczmarczyk/MoFem/html/_p_c_m_g_set_up_via_approx_orders_8cpp_source.html
I encounter following error, if restriction and interpolation operator is the same for multigrid solver. This happens when I use hierarchical approx. basis and my R matrix is has only values on diagonal. If I set PCMGSetRestriction it works great, but I believe that this is not optimal solution. Could you advise some solution for this? 176 // Petsc Development GIT revision: v3.5.3-1524-gee900cc GIT Date: 2015-01-31 17:44:15 -0600 178 // [0]PETSC ERROR: [0] MatPtAPSymbolic_MPIAIJ_MPIAIJ line 124 /opt/petsc/src/mat/impls/aij/mpi/mpiptap.c 179 // [0]PETSC ERROR: [0] MatPtAP_MPIAIJ_MPIAIJ line 80 /opt/petsc/src/mat/impls/aij/mpi/mpiptap.c 180 // [0]PETSC ERROR: [0] MatPtAP line 8458 /opt/petsc/src/mat/interface/matrix.c 181 // [0]PETSC ERROR: [0] PCSetUp_MG line 552 /opt/petsc/src/ksp/pc/impls/mg/mg.c 182 // [0]PETSC ERROR: [0] KSPSetUp line 220 /opt/petsc/src/ksp/ksp/interface/itfunc.c 183 184 // ==4284== Invalid read of size 8 185 // ==4284== at 0x5CAC873: MatPtAPSymbolic_MPIAIJ_MPIAIJ (mpiptap.c:154) 186 // ==4284== by 0x5CABA74: MatPtAP_MPIAIJ_MPIAIJ (mpiptap.c:83) 187 // ==4284== by 0x5D566D1: MatPtAP (matrix.c:8537) 188 // ==4284== by 0x61F27D6: PCSetUp_MG (mg.c:642) 189 // ==4284== by 0x612BD9C: PCSetUp (precon.c:909) 190 // ==4284== by 0x62C1A51: KSPSetUp (itfunc.c:306) 191 // ==4284== by 0xB98326: main (elasticity.cpp:403) 192 // ==4284== Address 0x208 is not stack'd, malloc'd or (recently) free'd Kind regards, Lukasz
