original_docstring: 
          This method plots brachial, radial and ulnar WSS (POISEUILLE) signal.
          
docstring: This method plots brachial, radial and ulnar WSS (POISEUILLE) signal.
docstring_tokens: This method plots brachial radial and ulnar WSS ( POISEUILLE ) signal .
code: def PlotPWSSComparative(self):
          '''
          This method plots brachial, radial and ulnar WSS (POISEUILLE) signal.
          '''
          try:
              from matplotlib.pyplot import plot, xlabel, ylabel, title, legend, savefig, close
          except ImportError:
              MatPlotLibError()
          FlowBrachial = 0
          FlowRadial = 0
          FlowUlnar = 0
          WSSBrachial = 0
          WSSRadial = 0
          WSSUlnar = 0
          i = 0
          j = 0
          k = 0
          for element in self.NetworkMesh.Elements:
              if element.Type is not 'Windkessel' and element.Name.find('rachial') != -1:
                  dofs = element.GetPoiseuilleDofs()
                  FlowBrachial = (self.Solutions[(self.DofMap.DofMap[element.Id, dofs[0]]),self.CardiacFreq*(self.Cycles-1):] - self.Solutions[(self.DofMap.DofMap[element.Id, dofs[1]]),self.CardiacFreq*(self.Cycles-1):])/element.R
                  WSSBrachial += ((4.0*element.eta)/pi) * (FlowBrachial/(mean(element.Radius))**3)
                  i += 1               
          
              if element.Type is not 'Windkessel' and element.Name.find('adial') != -1:
                  dofs = element.GetPoiseuilleDofs()
                  FlowRadial = (self.Solutions[(self.DofMap.DofMap[element.Id, dofs[0]]),self.CardiacFreq*(self.Cycles-1):] - self.Solutions[(self.DofMap.DofMap[element.Id, dofs[1]]),self.CardiacFreq*(self.Cycles-1):])/element.R
                  WSSRadial += ((4.0*element.eta)/pi) * (FlowRadial/(mean(element.Radius))**3)
                  j += 1
                  
              if element.Type is not 'Windkessel' and element.Name.find('lnar') != -1:
                  dofs = element.GetPoiseuilleDofs()
                  FlowUlnar = (self.Solutions[(self.DofMap.DofMap[element.Id, dofs[0]]),self.CardiacFreq*(self.Cycles-1):] - self.Solutions[(self.DofMap.DofMap[element.Id, dofs[1]]),self.CardiacFreq*(self.Cycles-1):])/element.R
                  WSSUlnar += ((4.0*element.eta)/pi) * (FlowUlnar/(mean(element.Radius))**3)
                  k += 1
          
          if i != 0:
              WSSBrachial = WSSBrachial/i
              print "Brachial"
              print "Wss(mean) = ", mean(WSSBrachial), " Pa", " Wss(peak) = ", max(WSSBrachial), " Pa"
              plot(self.t, WSSBrachial, 'r-', linewidth = 3, label = 'brachial')   #red line
          if i == 0:
              WSSBrachial = 0
              print "There isn't any element named brachial, check your network xml file"  
          if j != 0:
              WSSRadial = WSSRadial/j
              print "Radial"
              print "Wss(mean) = ", mean(WSSRadial), " Pa", " Wss(peak) = ", max(WSSRadial), " Pa"
              plot(self.t, WSSRadial, 'b-', linewidth = 3, label = 'radial')   #blue line
          if j == 0:
              WSSRadial = 0
              print "There isn't any element named radial, check your network xml file"
          if k != 0:
              WSSUlnar = WSSUlnar/k
              print "Ulnar"
              print "Wss(mean) = ", mean(WSSUlnar), " Pa", " Wss(peak) = ", max(WSSUlnar), " Pa"
              plot(self.t, WSSUlnar, 'g-', linewidth = 3, label = 'ulnar')   #green line
          if k == 0:
              WSSUlnar = 0
              print "There isn't any element named ulnar, check your network xml file"
          
          xlabel('Time ($s$)')
          ylabel('Wall Shear Stress ($Pa$)')
          title ('Brachial, Radial and Ulnar Wall Shear Stress')    
          legend()
          savefig(self.images+'brach_rad_uln_wss.png')    
          close()
code_tokens: def PlotPWSSComparative ( self ) : try : from matplotlib . pyplot import plot , xlabel , ylabel , title , legend , savefig , close except ImportError : MatPlotLibError ( ) FlowBrachial = 0 FlowRadial = 0 FlowUlnar = 0 WSSBrachial = 0 WSSRadial = 0 WSSUlnar = 0 i = 0 j = 0 k = 0 for element in self . NetworkMesh . Elements : if element . Type is not 'Windkessel' and element . Name . find ( 'rachial' ) != - 1 : dofs = element . GetPoiseuilleDofs ( ) FlowBrachial = ( self . Solutions [ ( self . DofMap . DofMap [ element . Id , dofs [ 0 ] ] ) , self . CardiacFreq * ( self . Cycles - 1 ) : ] - self . Solutions [ ( self . DofMap . DofMap [ element . Id , dofs [ 1 ] ] ) , self . CardiacFreq * ( self . Cycles - 1 ) : ] ) / element . R WSSBrachial += ( ( 4.0 * element . eta ) / pi ) * ( FlowBrachial / ( mean ( element . Radius ) ) ** 3 ) i += 1 if element . Type is not 'Windkessel' and element . Name . find ( 'adial' ) != - 1 : dofs = element . GetPoiseuilleDofs ( ) FlowRadial = ( self . Solutions [ ( self . DofMap . DofMap [ element . Id , dofs [ 0 ] ] ) , self . CardiacFreq * ( self . Cycles - 1 ) : ] - self . Solutions [ ( self . DofMap . DofMap [ element . Id , dofs [ 1 ] ] ) , self . CardiacFreq * ( self . Cycles - 1 ) : ] ) / element . R WSSRadial += ( ( 4.0 * element . eta ) / pi ) * ( FlowRadial / ( mean ( element . Radius ) ) ** 3 ) j += 1 if element . Type is not 'Windkessel' and element . Name . find ( 'lnar' ) != - 1 : dofs = element . GetPoiseuilleDofs ( ) FlowUlnar = ( self . Solutions [ ( self . DofMap . DofMap [ element . Id , dofs [ 0 ] ] ) , self . CardiacFreq * ( self . Cycles - 1 ) : ] - self . Solutions [ ( self . DofMap . DofMap [ element . Id , dofs [ 1 ] ] ) , self . CardiacFreq * ( self . Cycles - 1 ) : ] ) / element . R WSSUlnar += ( ( 4.0 * element . eta ) / pi ) * ( FlowUlnar / ( mean ( element . Radius ) ) ** 3 ) k += 1 if i != 0 : WSSBrachial = WSSBrachial / i print "Brachial" print "Wss(mean) = " , mean ( WSSBrachial ) , " Pa" , " Wss(peak) = " , max ( WSSBrachial ) , " Pa" plot ( self . t , WSSBrachial , 'r-' , linewidth = 3 , label = 'brachial' ) if i == 0 : WSSBrachial = 0 print "There isn't any element named brachial, check your network xml file" if j != 0 : WSSRadial = WSSRadial / j print "Radial" print "Wss(mean) = " , mean ( WSSRadial ) , " Pa" , " Wss(peak) = " , max ( WSSRadial ) , " Pa" plot ( self . t , WSSRadial , 'b-' , linewidth = 3 , label = 'radial' ) if j == 0 : WSSRadial = 0 print "There isn't any element named radial, check your network xml file" if k != 0 : WSSUlnar = WSSUlnar / k print "Ulnar" print "Wss(mean) = " , mean ( WSSUlnar ) , " Pa" , " Wss(peak) = " , max ( WSSUlnar ) , " Pa" plot ( self . t , WSSUlnar , 'g-' , linewidth = 3 , label = 'ulnar' ) if k == 0 : WSSUlnar = 0 print "There isn't any element named ulnar, check your network xml file" xlabel ( 'Time ($s$)' ) ylabel ( 'Wall Shear Stress ($Pa$)' ) title ( 'Brachial, Radial and Ulnar Wall Shear Stress' ) legend ( ) savefig ( self . images + 'brach_rad_uln_wss.png' ) close ( )
short_docstring: This method plots brachial, radial and ulnar WSS (POISEUILLE) signal.
short_docstring_tokens: This method plots brachial radial and ulnar WSS ( POISEUILLE ) signal .
comment: '''
          This method plots brachial, radial and ulnar WSS (POISEUILLE) signal.
          '''$SEP$#red line$SEP$#blue line$SEP$#green line