from cmf import fit_retention_curve raises the follwoing error: A

Import issue and runtime error with fit_retention_curve.py about cmf HOT 1 CLOSED

philippkraft commented on June 7, 2024

Import issue and runtime error with fit_retention_curve.py

from cmf.

Comments (1)

philippkraft commented on June 7, 2024

What version are you using? I fixed that problem in Feb 2013 and reviewed it in Oct 2014.

Actual code reads as:

cmf/cmf/fit_retention_curve.py

Lines 62 to 83 in 1fcd70f

 def make_vgm(params): 

 try: 

 vgm = cmf.VanGenuchtenMualem(1.0) 

 vgm.Phi = params[0] 

 vgm.alpha = params[1] 

 vgm.n = params[2] 

 if len(params) > 3: 

 vgm.theta_r = params[3] 

 return vgm 

 except RuntimeError: 

 return None 

 def get_error_vgm(params, pF, theta): 

 vgm = make_vgm(params) 

 if vgm: 

 model_theta = vgm.theta(vgm.Wetness_pF(pF)) 

 err = np.sum((theta - model_theta) ** 2) 

 var = np.sum((theta - np.mean(theta)) ** 2) 

 return err / var 

 else: 

 return 100.0 # return a really bad error for illegal cases

from cmf.

Related Issues (20)

Remove inline
Remove cmfmemory.h
Describe how to build a custom cmf
Update Docs with new Solver names HOT 1
Run CVodeKLU with CUDA
Setup.py build libraries on clean command
Freundlich Adsorption not working
Can't create a SoluteWaterIntegrator HOT 1
new flux: Schmitt-trigger HOT 4
make potential of Dirichlet BC a time series HOT 1
Undefined behaviour when c'tor of flux_connection throws HOT 2
Remove integratables HOT 1
CVodeKLU is not suitable in a SoluteWaterIntegrator HOT 1
Add additional ET method with parameterized aerodynamic resistance HOT 5
Diffusive singularity protection is inconsistent HOT 1
cmf.timeseries should get a from_pandas factory function HOT 2
Modernize the doxygen documentation HOT 1
Allow for Python flux connections, even if they are slow HOT 1
Simple infiltration ignores waterbalance connections regardless of their direction
Try out scikit-build with ninja for CMF 2.0

Import issue and runtime error with fit_retention_curve.py about cmf HOT 1 CLOSED

Comments (1)

Related Issues (20)

Recommend Projects

React

Vue.js

Typescript

TensorFlow

Django

Laravel

D3

Recommend Topics

javascript

web

server

Machine learning

Visualization

Game

Recommend Org

Facebook

Microsoft

Google

Alibaba

D3

Tencent

	def make_vgm(params):
	try:
	vgm = cmf.VanGenuchtenMualem(1.0)
	vgm.Phi = params[0]
	vgm.alpha = params[1]
	vgm.n = params[2]
	if len(params) > 3:
	vgm.theta_r = params[3]
	return vgm
	except RuntimeError:
	return None


	def get_error_vgm(params, pF, theta):
	vgm = make_vgm(params)
	if vgm:
	model_theta = vgm.theta(vgm.Wetness_pF(pF))
	err = np.sum((theta - model_theta) ** 2)
	var = np.sum((theta - np.mean(theta)) ** 2)
	return err / var
	else:
	return 100.0 # return a really bad error for illegal cases