#NEWLINE# def S_e( self , problem ) : #NEWLINE# #INDENT# '\n        Get the electrical field source\n        ' #NEWLINE# e_p = self . ePrimary ( problem ) #NEWLINE# Map_sigma_p = Maps . Vertical1DMap ( problem . mesh ) #NEWLINE# sigma_p = Map_sigma_p . _transform ( self . sigma1d ) #NEWLINE# if ( problem . mesh . dim == 1 ) : #NEWLINE# #INDENT# Mesigma = problem . mesh . getFaceInnerProduct ( problem . curModel . sigma ) #NEWLINE# Mesigma_p = problem . mesh . getFaceInnerProduct ( sigma_p ) #NEWLINE# #UNINDENT# if ( Mesigma . mesh . dim == 2 ) : #NEWLINE# #INDENT# pass #NEWLINE# #UNINDENT# if ( problem . mesh . dim == 3 ) : #NEWLINE# #INDENT# Mesigma = problem . MeSigma #NEWLINE# Mesigma_p = problem . mesh . getEdgeInnerProduct ( sigma_p ) #NEWLINE# #UNINDENT# return ( ( Mesigma - Mesigma_p ) * e_p )