Index: /issm/trunk/src/c/modules/Reduceloadx/Reduceloadx.cpp
===================================================================
--- /issm/trunk/src/c/modules/Reduceloadx/Reduceloadx.cpp	(revision 7190)
+++ /issm/trunk/src/c/modules/Reduceloadx/Reduceloadx.cpp	(revision 7191)
@@ -17,4 +17,5 @@
 	Vec         Kfsy_s = NULL;
 	int         Kfsm,Kfsn;
+	int         global_m,global_n;
 	PetscScalar a;
 	bool        fromlocalsize = true;
@@ -22,6 +23,7 @@
 
 	_printf_(VerboseModule(),"   Dirichlet lifting applied to load vector\n");
-	
-	if(pf && Kfs){
+
+	MatGetSize(Kfs,&global_m,&global_n);
+	if(pf && global_m*global_n){
 
 		/*pf = pf - Kfs * y_s;*/
Index: /issm/trunk/src/c/toolkits/petsc/patches/MatMultPatch.cpp
===================================================================
--- /issm/trunk/src/c/toolkits/petsc/patches/MatMultPatch.cpp	(revision 7190)
+++ /issm/trunk/src/c/toolkits/petsc/patches/MatMultPatch.cpp	(revision 7191)
@@ -24,4 +24,6 @@
 	int m,n;
 	Vec X_rel=NULL;
+
+	_assert_(A); _assert_(X);
 
 	if (MatMultCompatible(A,X)){