I have made a function that calculates the misfit(squared difference) between an observed signal at a target station and reconstructed signal from 10 arbitrary stations chosen from a dataset. There are a total of 105000 stations in the dataset. See the code below for clarity.
ds1 = ASDFDataSet("output1_custom/receivers.h5") #dataset from earthquake simulation Target = 500 #enter target station Total_trace = [Target,np.arrange(10)] #combination of target and reconstruction stations (needed as input for loaddata function) misfit=loaddata(Target,Total_trace,ds1) #loaddata function calculates the misfit
I was wondering if it is possible to find the combination of 10 stations (out of the 105000) that minimizes this misfit, with MCMC. I am brand new to the pymc3 package, and none of the tutorials appear to be relevant to what I am attempting.
Any comments/criticism is welcome!