def test_create_stratigraphy_PLI_pyamg(): data, bb = load_claudius() model = GeologicalModel(bb[0, :], bb[1, :]) model.set_model_data(data) s0 = model.create_and_add_foliation("s0", interpolatortype="PLI", nelements=1000, solver="pyamg", damp=True)
def test_remove_constraints_PLI(): data, bb = load_claudius() model = GeologicalModel(bb[0, :], bb[1, :]) model.set_model_data(data) s0 = model.create_and_add_foliation("s0", interpolatortype="FDI", nelements=1000, solver="cg", damp=False)
def test_create_stratigraphy_PLI_lu(): data, bb = load_claudius() model = GeologicalModel(bb[0, :], bb[1, :]) model.set_model_data(data) s0 = model.create_and_add_foliation('s0', interpolatortype='PLI', nelements=1000, solver='lu', damp=True)
def test_access_feature_model(): data, bb = load_claudius() model = GeologicalModel(bb[0, :], bb[1, :]) model.set_model_data(data) s0 = model.create_and_add_foliation("s0", interpolatortype="FDI", nelements=1000, solver="fake", damp=False) assert s0 == model["s0"]
from LoopStructural.visualisation import LavaVuModelViewer from LoopStructural.datasets import load_claudius #demo data from LoopStructural import log_to_file import pandas as pd import numpy as np ################################################################################################## # Specify a log file # ~~~~~~~~~~~~~~~~~~~~ log_to_file('logging_demo_log.log') ################################################################################################## # Create model # ~~~~~~~~~~~~~~~~~~~~ data, bb = load_claudius() model = GeologicalModel(bb[0,:],bb[1,:]) model.set_model_data(data) vals = [0,60,250,330,600] strat_column = {'strati':{}} for i in range(len(vals)-1): strat_column['strati']['unit_{}'.format(i)] = {'min':vals[i],'max':vals[i+1],'id':i} model.set_stratigraphic_column(strat_column) strati = model.create_and_add_foliation("strati", interpolatortype="FDI", # try changing this to 'PLI' nelements=1e4, # try changing between 1e3 and 5e4 buffer=0.3, solver='pyamg', damp=True )
Once you have the correct file for your system you will need to add an environment variable to windows :code:`SURFE` """ from LoopStructural import GeologicalModel from LoopStructural.visualisation import LavaVuModelViewer from LoopStructural.datasets import load_claudius #demo data import pandas as pd import glob import numpy as np import logging logging.getLogger().setLevel(logging.INFO) data, bb = load_claudius()#claudius.get_data() bb[1,0]+=200 bb[0,0]-=200 bb[1,1]+=200 bb[0,1]-=200 bb[1,2]+=200 bb[0,2]-=200 model = GeologicalModel(bb[0,:],bb[1,:]) data['random'] = np.random.random(data.shape[0]) model.set_model_data(data[data['random'] < 0.01])#[np.isnan(data['val'])]) strati = model.create_and_add_foliation("strati", interpolatortype="surfe", method='single_surface' ) print(strati.evaluate_value(model.regular_grid((10,10,10))))
def test_add_data(): data, bb = load_claudius() model = GeologicalModel(bb[0, :], bb[1, :]) model.set_model_data(data)
def test_create_model(): data, bb = load_claudius() model = GeologicalModel(bb[0, :], bb[1, :])