added functions to find average segment length, tortuoisity in a network

pranathivemuri
1 parent b3a38641
Showing 1 changed file with 67 additions and 0 deletions Show diff stats
stim/biomodels/network.h
@@ -90,6 +90,7 @@ class network{
 				return ss.str();
 			}
  
+
 	};
  
 protected:
@@ -109,6 +110,72 @@ public:
 		return V.size();
 	}
  
+	std::vector<vertex> operator*(T s){
+		for (unsigned i=0; i< vertices; i ++ ){
+			V[i] = V[i] * s;
+		}
+		return V;
+	}
+
+	std::vector<vertex> operator*(vec<T> s){
+		for (unsigned i=0; i< vertices; i ++ ){
+			for (unsigned dim = 0 ; dim< 3; dim ++){
+				V[i][dim] = V[i][dim] * s[dim];
+			}
+		}
+		return V;
+	}
+
+	// Returns an average of fiber/edge lengths in the network
+	double Lengths(){
+		stim::vec<T> L;double sumLength = 0;
+		for(unsigned e = 0; e < E.size(); e++){				//for each edge in the network
+			L.push_back(E[e].length());						//append the edge length
+			sumLength = sumLength + E[e].length();
+		}
+		double avg = sumLength / E.size();
+		return avg;
+	}
+
+
+	// Returns an average of tortuosities in the network
+	double Tortuosities(){
+		stim::vec<T> t;
+		stim::vec<T> id1, id2;                        // starting and ending vertices of the edge
+		//std::vector< stim::vec<T> > vert;	             // list of visited vertices
+		//std::vector<stim::vec<T>>::iterator it;			//create an iterator for searching the id2vert array
+		//std::vector<unsigned> daughters;               // to count number of daughter vessels at branch points
+		//stim::vec<T> visitedVert;
+		//unsigned visited;
+		double distance;double tortuosity;double sumTortuosity = 0;
+		for(unsigned e = 0; e < E.size(); e++){				//for each edge in the network
+			id1 = E[e][0];									//get the edge starting point
+			id2 = E[e][E[e].size() - 1];					//get the edge ending point
+			distance = (id1 - id2).len();                   //displacement between the starting and ending points
+			//it = find(vert.begin(), vert.end(), id1);	    //look for the first node has a daughter vessel already
+			//visitedVert = std::distance(vert.begin(), it);    //index where it already exists
+			//if(visitedVert == id1){
+			//	visited += 1;
+			//	daughters[e] = visited;
+			//}
+			//else{
+			//	vert.push_back(id1);
+			//	daughters[e] = 1;
+			//}
+			if(distance > 0){
+				tortuosity = E[e].length()/	distance	;		// tortuoisty = edge length / edge displacement
+			}
+			else{
+				tortuosity = 0;}
+			t.push_back(tortuosity);
+			sumTortuosity += tortuosity;
+		}
+		double avg = sumTortuosity / E.size();
+		//std::cout<<"number of unique branch points"<<daughters.size()<<std::endl;
+		return avg;
+	}
+
+
 	stim::cylinder<T> get_cylinder(unsigned f){
 		return E[f];					//return the specified edge (casting it to a fiber)
 	}