def get_constant(lin_op):
 constant = mul(lin_op, {})
 const_size = (constant.shape[0] * constant.shape[1])
 return np.reshape(constant, const_size, 'F')

def get_constant ( lin_op ) : constant = mul ( lin_op , { } ) const_size = ( constant . shape [ 0 ] * constant . shape [ 1 ] ) return np . reshape ( constant , const_size , 'F' )
returns the constant term in the expression .