Index: /issm/trunk/src/c/objects/Bamg/AdjacentTriangle.h
===================================================================
--- /issm/trunk/src/c/objects/Bamg/AdjacentTriangle.h	(revision 9370)
+++ /issm/trunk/src/c/objects/Bamg/AdjacentTriangle.h	(revision 9371)
@@ -18,5 +18,5 @@
 
 			//Constructors
-			AdjacentTriangle(){};
+			AdjacentTriangle():a(0),t(NULL) {};
 			AdjacentTriangle(Triangle* tt,int  aa): t(tt),a(aa &3) {};
 
Index: /issm/trunk/src/c/objects/Bamg/EigenMetric.cpp
===================================================================
--- /issm/trunk/src/c/objects/Bamg/EigenMetric.cpp	(revision 9370)
+++ /issm/trunk/src/c/objects/Bamg/EigenMetric.cpp	(revision 9371)
@@ -10,5 +10,5 @@
 	/*Constructor*/
 	/*FUNCTION EigenMetric::EigenMetric(const Metric M){{{1*/
-	EigenMetric::EigenMetric(const Metric M){
+	EigenMetric::EigenMetric(const Metric& M){
 		/*From a metric (a11,a21,a22), get eigen values lambda1 and lambda2 and one eigen vector v*/
 
@@ -90,5 +90,5 @@
 	/*}}}1*/
 	/*FUNCTION EigenMetric::EigenMetric(double r1,double r2,const D2 vp1){{{1*/
-	EigenMetric::EigenMetric(double r1,double r2,const D2 vp1): lambda1(r1),lambda2(r2),v(vp1){
+	EigenMetric::EigenMetric(double r1,double r2,const D2& vp1): lambda1(r1),lambda2(r2),v(vp1){
 
 	}/*}}}*/
Index: /issm/trunk/src/c/objects/Bamg/Metric.cpp
===================================================================
--- /issm/trunk/src/c/objects/Bamg/Metric.cpp	(revision 9370)
+++ /issm/trunk/src/c/objects/Bamg/Metric.cpp	(revision 9371)
@@ -14,5 +14,5 @@
 
 	/*Constructor/Destructor*/
-	Metric::Metric(){ };
+	/*FUNCTION Metric::Metric(){{{1*/
 	/*FUNCTION Metric::Metric(double a){{{1*/
 	Metric::Metric(double a): a11(1/(a*a)),a21(0),a22(1/(a*a)){
@@ -24,5 +24,5 @@
 	}/*}}}*/
 	/*FUNCTION Metric::Metric(const double  a[3],const  Metric& m0, const  Metric& m1,const  Metric&  m2 ){{{1*/ 
-	Metric::Metric(const double  a[3],const  Metric m0, const  Metric m1,const  Metric m2 ){
+	Metric::Metric(const double  a[3],const  Metric& m0, const  Metric& m1,const  Metric& m2 ){
 		/*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Metric.cpp/Metric)*/
 
@@ -45,5 +45,5 @@
 	/*}}}1*/
 	/*FUNCTION Metric::Metric(double  a,const  Metric& ma, double  b,const  Metric& mb){{{1*/
-	Metric::Metric(double  a,const  Metric ma, double  b,const  Metric mb) { 
+	Metric::Metric(double  a,const  Metric& ma, double  b,const  Metric& mb) { 
 		/*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Metric.cpp/EigenMetric)*/
 
@@ -139,5 +139,5 @@
 	/*Intermediary*/
 	/*FUNCTION LengthInterpole{{{1*/
-	double LengthInterpole(const Metric Ma,const  Metric Mb, R2 AB) {
+	double LengthInterpole(const Metric& Ma,const  Metric& Mb, R2 AB) {
 		/*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Metric.cpp/LengthInterpole)*/
 
@@ -310,5 +310,5 @@
 	/*}}}1*/
 	/*FUNCTION abscisseInterpole{{{1*/
-	double abscisseInterpole(const Metric Ma,const  Metric Mb, R2 AB,double s,int optim) { 
+	double abscisseInterpole(const Metric& Ma,const  Metric& Mb, R2 AB,double s,int optim) { 
 		/*Original code from Frederic Hecht <hecht@ann.jussieu.fr> (BAMG v1.01, Metric.cpp/abscisseInterpole)*/
 
Index: /issm/trunk/src/c/objects/Bamg/Metric.h
===================================================================
--- /issm/trunk/src/c/objects/Bamg/Metric.h	(revision 9370)
+++ /issm/trunk/src/c/objects/Bamg/Metric.h	(revision 9371)
@@ -25,10 +25,10 @@
 
 			//functions
-			Metric();
-			Metric(const EigenMetric);
+			Metric():a11(0),a21(0),a22(0){};
+			Metric(const EigenMetric&);
 			Metric(double a);
 			Metric(double a,double b,double c);
-			Metric( double  a,const  Metric ma, double  b,const  Metric mb);
-			Metric(const double  a[3],const  Metric m0,const  Metric m1,const  Metric m2 );
+			Metric( double  a,const  Metric& ma, double  b,const  Metric& mb);
+			Metric(const double  a[3],const  Metric& m0,const  Metric& m1,const  Metric& m2 );
 			void        Echo();
 			R2          mul(const R2 x)const;
@@ -62,6 +62,6 @@
 
 			//functions
-			EigenMetric(const Metric );
-			EigenMetric(double r1,double r2,const D2 vp1);
+			EigenMetric(const Metric& );
+			EigenMetric(double r1,double r2,const D2& vp1);
 			void   Echo();
 			void   Abs();
@@ -97,6 +97,6 @@
 	//Functions
 	void  SimultaneousMatrixReduction( Metric M1,  Metric M2,D2xD2 &V);
-	double LengthInterpole(const Metric Ma,const  Metric Mb, R2 AB);
-	double abscisseInterpole(const Metric Ma,const  Metric Mb, R2 AB,double s,int optim=0);
+	double LengthInterpole(const Metric& Ma,const  Metric& Mb, R2 AB);
+	double abscisseInterpole(const Metric& Ma,const  Metric& Mb, R2 AB,double s,int optim=0);
 
 	//inlines
@@ -115,5 +115,5 @@
 		}
 	}
-	inline Metric::Metric(const EigenMetric M) {
+	inline Metric::Metric(const EigenMetric& M) {
 		double v00=M.v.x*M.v.x;
 		double v11=M.v.y*M.v.y;
Index: /issm/trunk/src/c/objects/Bamg/R2.h
===================================================================
--- /issm/trunk/src/c/objects/Bamg/R2.h	(revision 9370)
+++ /issm/trunk/src/c/objects/Bamg/R2.h	(revision 9371)
@@ -68,5 +68,5 @@
 			  P2<R,RR>    ty() {return P2<R,RR>(x.y,y.y);} 
 			  //Operators
-			  P2<R,RR>     operator*(const P2<R,RR> c) const {return P2<R,RR>(x.x*c.x + x.y*c.y, y.x*c.x + y.y*c.y);}
+			  P2<R,RR>     operator*(const P2<R,RR>& c) const {return P2<R,RR>(x.x*c.x + x.y*c.y, y.x*c.x + y.y*c.y);}
 			  P2xP2<R,RR>  operator*(P2xP2<R,RR> c) const{
 				  return  P2xP2<R,RR>(x.x*c.x.x + x.y*c.y.x,
Index: /issm/trunk/src/c/objects/Bamg/Triangle.cpp
===================================================================
--- /issm/trunk/src/c/objects/Bamg/Triangle.cpp	(revision 9370)
+++ /issm/trunk/src/c/objects/Bamg/Triangle.cpp	(revision 9371)
@@ -9,4 +9,9 @@
 
 	/*Constructors/Destructors*/
+	/*FUNCTION Triangle(){{{1*/
+	Triangle::Triangle(void){
+
+	}
+	/*}}}*/
 	/*FUNCTION Triangle(Mesh *Th,long i,long j,long k) {{{1*/
 	Triangle::Triangle(Mesh *Th,long i,long j,long k) {
Index: /issm/trunk/src/c/objects/Bamg/Triangle.h
===================================================================
--- /issm/trunk/src/c/objects/Bamg/Triangle.h	(revision 9370)
+++ /issm/trunk/src/c/objects/Bamg/Triangle.h	(revision 9371)
@@ -29,5 +29,5 @@
 
 			//Constructors/Destructors
-			Triangle(){}
+			Triangle();
 			Triangle(Mesh *Th,long i,long j,long k);
 			Triangle(BamgVertex *v0,BamgVertex *v1,BamgVertex *v2);
Index: /issm/trunk/src/c/objects/Bamg/VertexOnEdge.h
===================================================================
--- /issm/trunk/src/c/objects/Bamg/VertexOnEdge.h	(revision 9370)
+++ /issm/trunk/src/c/objects/Bamg/VertexOnEdge.h	(revision 9371)
@@ -19,6 +19,6 @@
 
 			//Constructors
-			VertexOnEdge(BamgVertex * w, Edge *bw,double s) :v(w),be(bw),abcisse(s) {}
-			VertexOnEdge(){}
+			VertexOnEdge(BamgVertex * w, Edge *bw,double s) :v(w),be(bw),abcisse(s) {};
+			VertexOnEdge(){};
 
 			//Operators
Index: /issm/trunk/src/c/objects/Bamg/VertexOnGeom.cpp
===================================================================
--- /issm/trunk/src/c/objects/Bamg/VertexOnGeom.cpp	(revision 9370)
+++ /issm/trunk/src/c/objects/Bamg/VertexOnGeom.cpp	(revision 9371)
@@ -16,5 +16,5 @@
 		curvilincoord=0;
 		gv=0;
-		//ge=0;
+		ge=0;
 	} 
 	/*}}}*/
@@ -24,5 +24,5 @@
 		curvilincoord=-1;
 		gv=&g;
-		//ge=0;
+		ge=0;
 	}
 	/*}}}*/
@@ -32,5 +32,5 @@
 		curvilincoord=s;
 		ge=&g;
-		//gv=0;
+		gv=0;
 	}
 	/*}}}*/
Index: /issm/trunk/src/c/objects/Bamg/VertexOnVertex.cpp
===================================================================
--- /issm/trunk/src/c/objects/Bamg/VertexOnVertex.cpp	(revision 9370)
+++ /issm/trunk/src/c/objects/Bamg/VertexOnVertex.cpp	(revision 9371)
@@ -12,6 +12,6 @@
 	/*FUNCTION VertexOnVertex::VertexOnVertex(){{{1*/
 	VertexOnVertex::VertexOnVertex() {
-		//v=NULL;
-		//bv=NULL;
+		v=NULL;
+		bv=NULL;
 	};/*}}}*/
 	/*FUNCTION VertexOnVertex::VertexOnVertex(BamgVertex * w,BamgVertex *bw){{{1*/