Пример #1
0
 def getInverse(self, s):
     """
     returns the value of the inverse of the mapping for s
     """
     # self.logger.info("m0:"+str(log(s[0]**2+s[1]**2)/2))
     return (log(s[0]**2 + s[1]**2) / 2 -
             self.Mr) * [1., 0] + (atan2(s[1], s[0]) - self.Mi) * [0, 1.]
Пример #2
0
    def __init__(self, V_prior, Q_prior):
        """
        initializes the mapping

        :param V_prior: a-priori p-wave velocity
        :param Q_prior: a-priori Q-index (must be positive)
        """
        over2Q = 1.0 / (2 * Q_prior)
        # sigma_prior=1/(V_prior*(1-I*over2Q))**2 = 1/( V_prior * (1+over2Q**2)) **2 * ( (1-over2Q**2) + I * 2* over2Q )
        self.Mr = log(sqrt((1 - over2Q ** 2) ** 2 + (2 * over2Q) ** 2) / (V_prior * (1 + over2Q ** 2)) ** 2)
        self.Mi = atan2(2 * over2Q, 1 - over2Q ** 2)
Пример #3
0
    def __init__(self, V_prior, Q_prior):
        """
        initializes the mapping

        :param V_prior: a-priori p-wave velocity
        :param Q_prior: a-priori Q-index (must be positive)
        """
        over2Q = 1. / (2 * Q_prior)
        # sigma_prior=1/(V_prior*(1-I*over2Q))**2 = 1/( V_prior * (1+over2Q**2)) **2 * ( (1-over2Q**2) + I * 2* over2Q )
        self.Mr = log(
            sqrt((1 - over2Q**2)**2 + (2 * over2Q)**2) / (V_prior *
                                                          (1 + over2Q**2))**2)
        self.Mi = atan2(2 * over2Q, 1 - over2Q**2)
Пример #4
0
 def getInverse(self, s):
     """
     returns the value of the inverse of the mapping for s
     """
     # self.logger.info("m0:"+str(log(s[0]**2+s[1]**2)/2))
     return (log(s[0] ** 2 + s[1] ** 2) / 2 - self.Mr) * [1.0, 0] + (atan2(s[1], s[0]) - self.Mi) * [0, 1.0]