ralf/Calcium_Model
1
0
Fork 0
Code Issues 2 Pull Requests Actions Packages Projects Releases Wiki Activity

Compare commits

base: ralf:a2bb97fe22c9266e519b057da81ee785a5050a21
Branches Tags
ralf:main
...
compare: ralf:06a2b217d49a194c52d9f801da61783644c19f82
Branches Tags
ralf:main

2 Commits
a2bb97fe22 ... 06a2b217d4

Author SHA1 Message Date
ralf
06a2b217d4 formatting 2024-08-13 11:28:05 +03:00
ralf
2a687571c6 adjusted tolerances, with little to no effect on results 2024-08-13 11:15:44 +03:00
1 changed files with 10 additions and 1 deletions
Show Stats Expand all files Collapse all files

11
fitter.py
View File

@ -100,6 +100,7 @@ class Fitter:
def optimize(self, init_parameters=None): def optimize(self, init_parameters=None):
t0 = time.time() t0 = time.time()
self.iteration = 0 self.iteration = 0
if init_parameters is None: if init_parameters is None:
m = self.model() m = self.model()
K_pc_half = m.K_pc_half K_pc_half = m.K_pc_half
@ -118,7 +119,15 @@ class Fitter:
(10, 100, 100, 100, 10, 10), (10, 100, 100, 100, 10, 10),
) )
res = least_squares(self.cost_func, init_parameters, bounds=bounds, xtol=1e-10) res = least_squares(
self.cost_func,
init_parameters,
bounds=bounds,
xtol=1e-10,
ftol=1e-12,
verbose=2,
)
print() print()
print(" Parameters: [gGaL, ECal, K_pc_half, tau_xfer, tau_RC, offset]") print(" Parameters: [gGaL, ECal, K_pc_half, tau_xfer, tau_RC, offset]")
print(" Initial:", init_parameters.tolist()) print(" Initial:", init_parameters.tolist())