def __init__(self):
self.model = vmodel.model1d()
self.indata = data.data1d()
self.eigkR = eigenkernel.eigkernel()
self.eigkL = eigenkernel.eigkernel()
self.hArr = np.array([], dtype=np.float32)
self.disprefR = False
self.disprefL = False
return
def __init__(self):
self.model = vmodel.model1d()
self.data = data.data1d()
self.fs = 40.
self.slowness = 0.06
self.gausswidth = 2.5
self.amplevel = 0.005
self.t0 = 0.
return
def __init__(self,
factor=1.,
thresh=0.5,
waterdepth=-1.,
vpwater=1.5,
stdfactor=2.):
self.data = data.data1d()
self.factor = factor
self.thresh = thresh
self.stdfactor = stdfactor
# models
self.avg_model = vmodel.model1d()
self.min_model = vmodel.model1d()
self.init_model = vmodel.model1d()
self.real_model = vmodel.model1d()
self.temp_model = vmodel.model1d()
#
self.vprfwrd = vprofile.vprofile1d()
self.waterdepth = waterdepth
self.vpwater = vpwater
#
self.avg_misfit = 0.
self.code = ''
return
import data
import numpy as np
import multiprocessing
from functools import partial
disp = data.disp()
disp.readdisptxt('../disp_-112.0_36.0_ray.txt')
disp.pvelp = np.float32(np.random.randn(20)) + disp.pvelo
disp.readaziphitxt('../aziphi_-112.0_36.0.txt')
disp.pphip = np.float32(np.random.randn(20)) + disp.pphio
disp.readaziamptxt('../aziamp_-112.0_36.0.txt')
disp.pampp = np.float32(np.random.randn(20)) + disp.pampo
rf = data.rf()
rf.readrftxt('../old_code/TEST/in.rf')
data1d = data.data1d()
data1d.rfr.readrftxt('../old_code/TEST/in.rf')
data1d.dispR.readdisptxt('../disp_-112.0_36.0_ray.txt')
rf.tp = rf.to[:]
rf.rfp = rf.rfo[:]
#
# import time
# def testmp(i, rf):
# print i
# rf.get_misfit_incompatible()
#
# ilst = (np.arange(10000000)).tolist()
# TESTmp = partial(testmp, rf=rf)
def __init__(self):
self.model = vmodel.model1d()
self.indata = data.data1d()
self.eigkR = eigenkernel.eigkernel()
self.eigkL = eigenkernel.eigkernel()
return
import data indata = data.data1d() data.readdisp(infname='old_code/TEST/Q22A.com.txt', indisp=indata.dispR)