def f_strain_spline3_abrupt_lmfit(alt, pars, t):
w_strain = pars[2:int(pars[0])+2]
N_abscisses = len(w_strain) - 3.
z = alt * N_abscisses / t
index = 0
strain = ones(len(z))
for i in z:
strain[index] = cubicSpline(i, w_strain) / 100.
index = index + 1
return strain
def f_DW_spline3_abrupt_lmfit(alt, pars, t):
w_DW = pars[-int(pars[1]):]
N_abscisses = len(w_DW) - 3.
z = alt * N_abscisses / t
index = 0
DW = ones(len(z))
for i in z:
DW[index] = cubicSpline(i, w_DW)
index = index + 1
return DW
def f_strain_spline3_abrupt(alt, sp, t):
w_strain = sp[:]
N_abscisses = len(w_strain) - 3.
z = alt * N_abscisses / t
index = 0
strain = ones(len(z))
for i in z:
strain[index] = cubicSpline(i, w_strain) / 100.
index = index + 1
return strain
def f_DW_spline3_abrupt_lmfit(alt, pars, t):
w_DW = pars[-int(pars[1]):]
N_abscisses = len(w_DW) - 3.
z = alt * N_abscisses / t
index = 0
DW = ones(len(z))
for i in z:
DW[index] = cubicSpline(i, w_DW)
index = index + 1
return DW
def f_DW_spline3_abrupt(alt, dwp, t):
w_DW = dwp[:]
N_abscisses = len(w_DW) - 3.
z = alt * N_abscisses / t
index = 0
DW = ones(len(z))
for i in z:
DW[index] = cubicSpline(i, w_DW)
index = index + 1
return DW
def f_strain_spline3_abrupt_lmfit(alt, pars, t):
w_strain = pars[2:int(pars[0]) + 2]
N_abscisses = len(w_strain) - 3.
z = alt * N_abscisses / t
index = 0
strain = ones(len(z))
for i in z:
strain[index] = cubicSpline(i, w_strain) / 100.
index = index + 1
return strain
def f_DW_spline3_abrupt(alt, dwp, t):
w_DW = dwp[:]
N_abscisses = len(w_DW) - 3.
z = alt * N_abscisses / t
index = 0
DW = ones(len(z))
for i in z:
DW[index] = cubicSpline(i, w_DW)
index = index + 1
return DW
def f_strain_spline3_abrupt(alt, sp, t):
w_strain = sp[:]
N_abscisses = len(w_strain) - 3.
z = alt * N_abscisses / t
index = 0
strain = ones(len(z))
for i in z:
strain[index] = cubicSpline(i, w_strain) / 100.
index = index + 1
return strain
def f_strain_spline3_smooth_lmfit(alt, pars, t):
w_strain_free = pars[2:int(pars[0])+2]
w_strain = array([0.0, 0.0, 0.0])
w_strain = append(w_strain, w_strain_free)
N_abscisses = len(w_strain) - 3.
z = alt * N_abscisses / t
index = 0
strain = ones(len(z))
for i in z:
strain[index] = cubicSpline(i, w_strain) / 100.
index = index + 1
return strain
def f_strain_spline3_smooth(alt, sp, t):
w_strain_free = sp[:]
w_strain = array([0.0, 0.0, 0.0])
w_strain = append(w_strain, w_strain_free)
N_abscisses = len(w_strain) - 3.
z = alt * N_abscisses / t
index = 0
strain = ones(len(z))
for i in z:
strain[index] = cubicSpline(i, w_strain) / 100.
index = index + 1
return strain
def f_DW_spline3_smooth_lmfit(alt, pars, t):
w_DW_free = pars[-int(pars[1]):]
w_DW = array([1.0, 1.0, 1.0])
w_DW = append(w_DW, w_DW_free)
N_abscisses = len(w_DW) - 3.
z = alt * N_abscisses / t
index = 0
DW = ones(len(z))
for i in z:
DW[index] = cubicSpline(i, w_DW)
index = index + 1
return DW
def f_DW_spline3_smooth(alt, dwp, t):
w_DW_free = dwp[:]
w_DW = array([1.0, 1.0, 1.0])
w_DW = append(w_DW, w_DW_free)
N_abscisses = len(w_DW) - 3.
z = alt * N_abscisses / t
index = 0
DW = ones(len(z))
for i in z:
DW[index] = cubicSpline(i, w_DW)
index = index + 1
return DW
def f_DW_spline3_smooth_lmfit(alt, pars, t):
w_DW_free = pars[-int(pars[1]):]
w_DW = array([1.0, 1.0, 1.0])
w_DW = append(w_DW, w_DW_free)
N_abscisses = len(w_DW) - 3.
z = alt * N_abscisses / t
index = 0
DW = ones(len(z))
for i in z:
DW[index] = cubicSpline(i, w_DW)
index = index + 1
return DW
def f_DW_spline3_smooth(alt, dwp, t):
w_DW_free = dwp[:]
w_DW = array([1.0, 1.0, 1.0])
w_DW = append(w_DW, w_DW_free)
N_abscisses = len(w_DW) - 3.
z = alt * N_abscisses / t
index = 0
DW = ones(len(z))
for i in z:
DW[index] = cubicSpline(i, w_DW)
index = index + 1
return DW
def f_strain_spline3_smooth_lmfit(alt, pars, t):
w_strain_free = pars[2:int(pars[0]) + 2]
w_strain = array([0.0, 0.0, 0.0])
w_strain = append(w_strain, w_strain_free)
N_abscisses = len(w_strain) - 3.
z = alt * N_abscisses / t
index = 0
strain = ones(len(z))
for i in z:
strain[index] = cubicSpline(i, w_strain) / 100.
index = index + 1
return strain
def f_strain_spline3_smooth(alt, sp, t):
w_strain_free = sp[:]
w_strain = array([0.0, 0.0, 0.0])
w_strain = append(w_strain, w_strain_free)
N_abscisses = len(w_strain) - 3.
z = alt * N_abscisses / t
index = 0
strain = ones(len(z))
for i in z:
strain[index] = cubicSpline(i, w_strain) / 100.
index = index + 1
return strain
