def __init__(self, parameterizer):
"""no more constraints than parameters"""
LogUniND.__init__(self,
vinfs=parameterizer.MINF,
vsups=parameterizer.MSUP,
k=self.k,
nanbehavior=self.nanbehavior)
def __init__(self, parameterizer, dvsmin, dvsmax):
assert dvsmin < dvsmax
self.p = parameterizer
NLAYER = self.p.NLAYER
MINF = np.concatenate((self.p.MINF, np.ones(NLAYER - 1) * dvsmin))
MSUP = np.concatenate((self.p.MSUP, np.ones(NLAYER - 1) * dvsmax))
LogUniND.__init__(self,
vinfs=MINF,
vsups=MSUP,
k=self.k,
nanbehavior=self.nanbehavior)
def __call__(self, m):
_, VP, VS, RH = self.p.inv(m)
DVP = VS[1:] - VS[:-1]
DVS = VP[1:] - VP[:-1]
DRH = RH[1:] - RH[:-1]
extended_m = np.concatenate((m, DVP, DVS, DRH))
return LogUniND.__call__(self, extended_m)
def __call__(self, m):
_, _, VS, _ = self.p.inv(m)
DVS = VS[1:] - VS[:-1]
extended_m = np.concatenate((m, DVS))
return LogUniND.__call__(self, extended_m)