def updatePolarity(self):
    '\n        Fills in the polarity attribute of the interpolation, determining whether\n        the "plus" (True) or "minus" (False) solution of the system of equations\n        should be used for each sector.  Needs to be called in __init__.\n\n        Parameters\n        ----------\n        none\n\n        Returns\n        -------\n        none\n        '
    x_temp = (0.5 * (self.x_values[0:(self.x_n - 1), 0:(self.y_n - 1)] + self.x_values[1:self.x_n, 1:self.y_n]))
    y_temp = (0.5 * (self.y_values[0:(self.x_n - 1), 0:(self.y_n - 1)] + self.y_values[1:self.x_n, 1:self.y_n]))
    size = ((self.x_n - 1) * (self.y_n - 1))
    x_temp = np.reshape(x_temp, size)
    y_temp = np.reshape(y_temp, size)
    y_pos = np.tile(np.arange(0, (self.y_n - 1)), (self.x_n - 1))
    x_pos = np.reshape(np.tile(np.arange(0, (self.x_n - 1)), ((self.y_n - 1), 1)).transpose(), size)
    self.polarity = np.ones(((self.x_n - 1), (self.y_n - 1)), dtype=bool)
    (alpha, beta) = self.findCoords(x_temp, y_temp, x_pos, y_pos)
    polarity = np.logical_and(np.logical_and((alpha > 0), (alpha < 1)), np.logical_and((beta > 0), (beta < 1)))
    self.polarity = np.reshape(polarity, ((self.x_n - 1), (self.y_n - 1)))
Fills in the polarity attribute of the interpolation, determining whether
the "plus" (True) or "minus" (False) solution of the system of equations
should be used for each sector.  Needs to be called in __init__.

Parameters
----------
none

Returns
-------
none
def updatePolarity(self):
    '\n        Fills in the polarity attribute of the interpolation, determining whether\n        the "plus" (True) or "minus" (False) solution of the system of equations\n        should be used for each sector.  Needs to be called in __init__.\n\n        Parameters\n        ----------\n        none\n\n        Returns\n        -------\n        none\n        '
    x_temp = (0.5 * (self.x_values[0:(self.x_n - 1), 0:(self.y_n - 1)] + self.x_values[1:self.x_n, 1:self.y_n]))
    y_temp = (0.5 * (self.y_values[0:(self.x_n - 1), 0:(self.y_n - 1)] + self.y_values[1:self.x_n, 1:self.y_n]))
    size = ((self.x_n - 1) * (self.y_n - 1))
    x_temp = np.reshape(x_temp, size)
    y_temp = np.reshape(y_temp, size)
    y_pos = np.tile(np.arange(0, (self.y_n - 1)), (self.x_n - 1))
    x_pos = np.reshape(np.tile(np.arange(0, (self.x_n - 1)), ((self.y_n - 1), 1)).transpose(), size)
    self.polarity = np.ones(((self.x_n - 1), (self.y_n - 1)), dtype=bool)
    (alpha, beta) = self.findCoords(x_temp, y_temp, x_pos, y_pos)
    polarity = np.logical_and(np.logical_and((alpha > 0), (alpha < 1)), np.logical_and((beta > 0), (beta < 1)))
    self.polarity = np.reshape(polarity, ((self.x_n - 1), (self.y_n - 1)))