Commit edcfec7b authored by Philipp Frank's avatar Philipp Frank

option remove asperity

parent 002494dd
Pipeline #77544 passed with stages
in 13 minutes and 20 seconds
......@@ -77,19 +77,20 @@ class SimpleCorrelatedField(Operator):
sig_flex = vflex @ expander @ flex
xi = ducktape(dom, None, prefix + 'spectrum')
if asperity is None:
raise NotImplementedError
asp = LognormalTransform(*asperity, prefix + 'asperity', 0)
vasp = np.zeros(dom.shape)
vasp[0] = 1
vasp = makeOp(makeField(dom, vasp))
sig_asp = vasp @ expander @ asp
shift = np.ones(dom.shape)
shift[0] = _log_vol(pspace)**2/12.
shift = makeField(dom, shift)
asp = (Adder(shift) @ sig_asp).ptw("sqrt")
a = a + _SlopeRemover(pspace, 0) @ twolog @ (xi*sig_flex*asp)
if asperity is None:
asp = makeOp(shift.ptw("sqrt"))
a = a + _SlopeRemover(pspace, 0) @ twolog @ asp @ (xi*sig_flex)
else:
asp = LognormalTransform(*asperity, prefix + 'asperity', 0)
vasp = np.zeros(dom.shape)
vasp[0] = 1
vasp = makeOp(makeField(dom, vasp))
sig_asp = vasp @ expander @ asp
asp = (Adder(shift) @ sig_asp).ptw("sqrt")
a = a + _SlopeRemover(pspace, 0) @ twolog @ (xi*sig_flex*asp)
a = _Normalization(pspace, 0) @ a
maskzm = np.ones(pspace.shape)
maskzm[0] = 0
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment