Index: /issm/trunk-jpl/src/c/classes/matrix/ElementMatrix.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/matrix/ElementMatrix.cpp	(revision 26916)
+++ /issm/trunk-jpl/src/c/classes/matrix/ElementMatrix.cpp	(revision 26917)
@@ -260,4 +260,22 @@
 	_printf_("   sglobaldoflist  (" << sglobaldoflist << "): ");
 	for(i=0;i<nrows;i++)_printf_(" " << sglobaldoflist[i]); _printf_(" \n");
+}
+/*}}}*/
+bool ElementMatrix::HasDof(int dof,int set){/*{{{*/
+
+	if(set==FsetEnum){
+		for(int i=0;i<this->nrows;i++) if(this->fglobaldoflist[i] == dof) return true;
+	}
+	else if(set==GsetEnum){
+		for(int i=0;i<this->nrows;i++) if(this->gglobaldoflist[i] == dof) return true;
+	}
+	else if(set==SsetEnum){
+		for(int i=0;i<this->nrows;i++) if(this->sglobaldoflist[i] == dof) return true;
+	}
+	else{
+		_error_("not supported yet");
+	}
+
+	return false;
 }
 /*}}}*/
Index: /issm/trunk-jpl/src/c/classes/matrix/ElementMatrix.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/matrix/ElementMatrix.h	(revision 26916)
+++ /issm/trunk-jpl/src/c/classes/matrix/ElementMatrix.h	(revision 26917)
@@ -46,4 +46,5 @@
 		void StaticCondensation(int numindices,int* indices);
 		void Transpose(void);
+		bool HasDof(int dof,int set);
 };
 #endif //#ifndef _ELEMENT_MATRIX_H_