Index: /issm/trunk-jpl/src/c/classes/matrix/Vector.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/matrix/Vector.h	(revision 14432)
+++ /issm/trunk-jpl/src/c/classes/matrix/Vector.h	(revision 14433)
@@ -158,5 +158,5 @@
 		/*Vector specific routines*/
 		/*FUNCTION Echo{{{*/
-		void Echo(void){
+		void Echo(void){_assert_(this);
 
 			if(type==PetscVecType){
@@ -175,5 +175,5 @@
 		/*}}}*/
 		/*FUNCTION Assemble{{{*/
-		void Assemble(void){
+		void Assemble(void){_assert_(this);
 
 			if(type==PetscVecType){
@@ -192,6 +192,5 @@
 		/*}}}*/
 		/*FUNCTION SetValues{{{*/
-		void SetValues(int ssize, int* list, doubletype* values, InsMode mode){
-
+		void SetValues(int ssize, int* list, doubletype* values, InsMode mode){ _assert_(this);
 			if(type==PetscVecType){
 				#ifdef _HAVE_PETSC_
@@ -209,5 +208,5 @@
 		/*}}}*/
 		/*FUNCTION SetValue{{{*/
-		void SetValue(int dof, doubletype value, InsMode mode){
+		void SetValue(int dof, doubletype value, InsMode mode){_assert_(this);
 
 			if(type==PetscVecType){
@@ -226,5 +225,5 @@
 		/*}}}*/
 		/*FUNCTION GetValue{{{*/
-		void GetValue(doubletype* pvalue,int dof){
+		void GetValue(doubletype* pvalue,int dof){_assert_(this);
 
 			if(type==PetscVecType){
@@ -243,5 +242,5 @@
 		/*}}}*/
 		/*FUNCTION GetSize{{{*/
-		void GetSize(int* pM){
+		void GetSize(int* pM){_assert_(this);
 
 			if(type==PetscVecType){
@@ -260,5 +259,5 @@
 		/*}}}*/
 		/*FUNCTION IsEmpty{{{*/
-		bool IsEmpty(void){
+		bool IsEmpty(void){_assert_(this);
 
 			int M;
@@ -273,5 +272,5 @@
 		/*}}}*/
 		/*FUNCTION GetLocalSize{{{*/
-		void GetLocalSize(int* pM){
+		void GetLocalSize(int* pM){_assert_(this);
 
 			if(type==PetscVecType){
@@ -290,5 +289,5 @@
 		/*}}}*/
 		/*FUNCTION Duplicate{{{*/
-		Vector* Duplicate(void){
+		Vector* Duplicate(void){_assert_(this);
 
 			Vector* output=NULL;
@@ -313,5 +312,5 @@
 		/*}}}*/
 		/*FUNCTION Set{{{*/
-		void Set(doubletype value){
+		void Set(doubletype value){_assert_(this);
 
 			if(type==PetscVecType){
@@ -330,5 +329,5 @@
 		/*}}}*/
 		/*FUNCTION AXPY{{{*/
-		void AXPY(Vector* X, doubletype a){
+		void AXPY(Vector* X, doubletype a){_assert_(this);
 
 			if(type==PetscVecType){
@@ -347,5 +346,5 @@
 		/*}}}*/
 		/*FUNCTION AYPX{{{*/
-		void AYPX(Vector* X, doubletype a){
+		void AYPX(Vector* X, doubletype a){_assert_(this);
 
 			if(type==PetscVecType){
@@ -364,5 +363,5 @@
 		/*}}}*/
 		/*FUNCTION ToMPISerial{{{*/
-		doubletype* ToMPISerial(void){
+		doubletype* ToMPISerial(void){_assert_(this);
 
 			doubletype* vec_serial=NULL;
@@ -385,5 +384,5 @@
 		/*}}}*/
 		/*FUNCTION Copy{{{*/
-		void Copy(Vector* to){
+		void Copy(Vector* to){_assert_(this);
 
 			if(type==PetscVecType){
@@ -402,5 +401,5 @@
 		/*}}}*/
 		/*FUNCTION Norm{{{*/
-		doubletype Norm(NormMode norm_type){
+		doubletype Norm(NormMode norm_type){_assert_(this);
 
 			doubletype norm=0;
@@ -422,5 +421,5 @@
 		/*}}}*/
 		/*FUNCTION Scale{{{*/
-		void Scale(doubletype scale_factor){
+		void Scale(doubletype scale_factor){_assert_(this);
 
 			if(type==PetscVecType){
@@ -439,5 +438,5 @@
 		/*}}}*/
 		/*FUNCTION Dot{{{*/
-		doubletype Dot(Vector* vector){
+		doubletype Dot(Vector* vector){_assert_(this);
 
 			doubletype dot;
@@ -459,5 +458,5 @@
 		/*}}}*/
 		/*FUNCTION PointwiseDivide{{{*/
-		void PointwiseDivide(Vector* x,Vector* y){
+		void PointwiseDivide(Vector* x,Vector* y){_assert_(this);
 
 			if(type==PetscVecType){