original_docstring: 
          This method plots brachial, radial and ulnar mean flow.
          If cycle is not specified, default cycle is the last one
          
docstring: This method plots brachial, radial and ulnar mean flow.
  If cycle is not specified, default cycle is the last one
docstring_tokens: This method plots brachial radial and ulnar mean flow . If cycle is not specified default cycle is the last one
code: def PlotFlowComparative(self):
          '''
          This method plots brachial, radial and ulnar mean flow.
          If cycle is not specified, default cycle is the last one
          '''
          try:
              from matplotlib.pyplot import plot, xlabel, ylabel, title, legend, savefig, close
          except ImportError:
              MatPlotLibError()
          FlowBrachial = 0
          FlowRadial = 0
          FlowUlnar = 0
          i = 0
          j = 0
          k = 0
                      
          for element in self.NetworkMesh.Elements:
              if element.Type is not 'Windkessel' and element.Name.find('brachial') != -1:
                  dofs = element.GetPoiseuilleDofs()
                  p1 = self.Solutions[(self.DofMap.DofMap[element.Id, dofs[0]]),self.CardiacFreq*(self.Cycles-1):]
                  p2 = self.Solutions[(self.DofMap.DofMap[element.Id, dofs[1]]),self.CardiacFreq*(self.Cycles-1):]
                  if mean(p1)>=mean(p2):
                      FlowBrachial += (p1-p2)/element.R
                  else:
                      FlowBrachial += (p2-p1)/element.R
                  i += 1               
          
              if element.Type is not 'Windkessel' and element.Name.find('radial') != -1:   
                  dofs = element.GetPoiseuilleDofs()
                  p1 = self.Solutions[(self.DofMap.DofMap[element.Id, dofs[0]]),self.CardiacFreq*(self.Cycles-1):]
                  p2 = self.Solutions[(self.DofMap.DofMap[element.Id, dofs[1]]),self.CardiacFreq*(self.Cycles-1):]
                  if mean(p1)>=mean(p2):
                      FlowRadial += (p1-p2)/element.R
                  else:
                      FlowRadial += (p2-p1)/element.R
                  j += 1
                  
              if element.Type is not 'Windkessel' and element.Name.find('ulnar') != -1: 
                  dofs = element.GetPoiseuilleDofs()
                  p1 = self.Solutions[(self.DofMap.DofMap[element.Id, dofs[0]]),self.CardiacFreq*(self.Cycles-1):]
                  p2 = self.Solutions[(self.DofMap.DofMap[element.Id, dofs[1]]),self.CardiacFreq*(self.Cycles-1):]
                  if mean(p1)>=mean(p2):
                      FlowUlnar += (p1-p2)/element.R
                  else:
                      FlowUlnar += (p2-p1)/element.R
                  k += 1
  
          if i != 0:
              FlowBrachial = (FlowBrachial*6e7)/i
              print "Brachial Flow = ", mean(FlowBrachial), "mL/min"
              plot(self.t, FlowBrachial, 'r-', linewidth = 3, label = 'brachial')   #red line
          if i == 0:
              FlowBrachial = 0
              print "There isn't any element named brachial, check your network xml file"  
          if j != 0:
              FlowRadial = (FlowRadial*6e7)/j 
              print "Radial Flow = ", mean(FlowRadial), "mL/min"
              plot(self.t, FlowRadial, 'b-', linewidth = 3, label = 'radial')   #blue line
          if j == 0:
              FlowRadial = 0
              print "There isn't any element named brachial, check your network xml file"  
          if k != 0:
              FlowUlnar = (FlowUlnar*6e7)/k
              print "Ulnar Flow = " , mean(FlowUlnar), "mL/min"
              plot(self.t, FlowUlnar, 'g-', linewidth = 3, label = 'ulnar')   #green line
          if k == 0:
              FlowUlnar = 0
              print "There isn't any element named brachial, check your network xml file"       
  
          xlabel('Time ($s$)')
          ylabel('Flow ($mL/min$)')
          title ('Brachial, Radial and Ulnar Flow Output')    
          legend()
          savefig(self.images+'brach_rad_uln_flow.png')
          close()
code_tokens: def PlotFlowComparative ( self ) : try : from matplotlib . pyplot import plot , xlabel , ylabel , title , legend , savefig , close except ImportError : MatPlotLibError ( ) FlowBrachial = 0 FlowRadial = 0 FlowUlnar = 0 i = 0 j = 0 k = 0 for element in self . NetworkMesh . Elements : if element . Type is not 'Windkessel' and element . Name . find ( 'brachial' ) != - 1 : dofs = element . GetPoiseuilleDofs ( ) p1 = self . Solutions [ ( self . DofMap . DofMap [ element . Id , dofs [ 0 ] ] ) , self . CardiacFreq * ( self . Cycles - 1 ) : ] p2 = self . Solutions [ ( self . DofMap . DofMap [ element . Id , dofs [ 1 ] ] ) , self . CardiacFreq * ( self . Cycles - 1 ) : ] if mean ( p1 ) >= mean ( p2 ) : FlowBrachial += ( p1 - p2 ) / element . R else : FlowBrachial += ( p2 - p1 ) / element . R i += 1 if element . Type is not 'Windkessel' and element . Name . find ( 'radial' ) != - 1 : dofs = element . GetPoiseuilleDofs ( ) p1 = self . Solutions [ ( self . DofMap . DofMap [ element . Id , dofs [ 0 ] ] ) , self . CardiacFreq * ( self . Cycles - 1 ) : ] p2 = self . Solutions [ ( self . DofMap . DofMap [ element . Id , dofs [ 1 ] ] ) , self . CardiacFreq * ( self . Cycles - 1 ) : ] if mean ( p1 ) >= mean ( p2 ) : FlowRadial += ( p1 - p2 ) / element . R else : FlowRadial += ( p2 - p1 ) / element . R j += 1 if element . Type is not 'Windkessel' and element . Name . find ( 'ulnar' ) != - 1 : dofs = element . GetPoiseuilleDofs ( ) p1 = self . Solutions [ ( self . DofMap . DofMap [ element . Id , dofs [ 0 ] ] ) , self . CardiacFreq * ( self . Cycles - 1 ) : ] p2 = self . Solutions [ ( self . DofMap . DofMap [ element . Id , dofs [ 1 ] ] ) , self . CardiacFreq * ( self . Cycles - 1 ) : ] if mean ( p1 ) >= mean ( p2 ) : FlowUlnar += ( p1 - p2 ) / element . R else : FlowUlnar += ( p2 - p1 ) / element . R k += 1 if i != 0 : FlowBrachial = ( FlowBrachial * 6e7 ) / i print "Brachial Flow = " , mean ( FlowBrachial ) , "mL/min" plot ( self . t , FlowBrachial , 'r-' , linewidth = 3 , label = 'brachial' ) if i == 0 : FlowBrachial = 0 print "There isn't any element named brachial, check your network xml file" if j != 0 : FlowRadial = ( FlowRadial * 6e7 ) / j print "Radial Flow = " , mean ( FlowRadial ) , "mL/min" plot ( self . t , FlowRadial , 'b-' , linewidth = 3 , label = 'radial' ) if j == 0 : FlowRadial = 0 print "There isn't any element named brachial, check your network xml file" if k != 0 : FlowUlnar = ( FlowUlnar * 6e7 ) / k print "Ulnar Flow = " , mean ( FlowUlnar ) , "mL/min" plot ( self . t , FlowUlnar , 'g-' , linewidth = 3 , label = 'ulnar' ) if k == 0 : FlowUlnar = 0 print "There isn't any element named brachial, check your network xml file" xlabel ( 'Time ($s$)' ) ylabel ( 'Flow ($mL/min$)' ) title ( 'Brachial, Radial and Ulnar Flow Output' ) legend ( ) savefig ( self . images + 'brach_rad_uln_flow.png' ) close ( )
short_docstring: This method plots brachial, radial and ulnar mean flow.
short_docstring_tokens: This method plots brachial radial and ulnar mean flow .
comment: '''
          This method plots brachial, radial and ulnar mean flow.
          If cycle is not specified, default cycle is the last one
          '''$SEP$#red line$SEP$#blue line$SEP$#green line