HajimeKawahara

HajimeKawahara

HajimeKawahara

HajimeKawahara

dat1 = pd.read_csv("/home/kawahara/gascell/data/SLD300K2.SPE", skiprows=8)
wav1 = dat1["WaveLength(nm) "]
val1 = dat1["Level(mW)     "]

HajimeKawahara

HajimeKawahara

def diad_merge_grids(grid1, grid2):
    """merge two different grids into one grid using diad
    
    Args:
        grid1: grid 1
        grid2: grid 2
        
    Returns:
        merged grid (len(grid1)*len(grid2),2)
        

HajimeKawahara

"""
a+b+c=3
"""

def invf(x,y):
    if x is not None and y is not None:
        return 3 - x - y
    else:
        return None
    

HajimeKawahara

import GPy
import numpy as np
import matplotlib.pyplot as plt
np.random.seed(seed=1)


kernel = GPy.kern.Matern32(2, ARD=True)
N = 100
X = np.random.uniform(-3.,3.,(N, 2))
Z = np.sin(X[:,0:1]) * np.sin(X[:,1:2]) + np.random.randn(N,1)*0.05

HajimeKawahara

import pandas as pd
import numpy as np
import matplotlib.pyplot as plt
import tqdm
from matplotlib.ticker import MultipleLocator, FormatStrFormatter

#filename_0="CTLv8_tyotto.csv"
#filename_1="CTLv8_tyotto.csv"

filename_0="CTLv8_0.csv"

HajimeKawahara

import pandas as pd
import numpy as np
from exojax.dynamics.rvfunc import rvf
import jax.numpy as jnp
from jax import random
from jax import vmap, jit
import matplotlib.pyplot as plt
to1=2459370 #fit offset
to2=2450000 #table Toffset
Toff=9370