コード例 #1
0
ファイル: omega.go プロジェクト: tflovorn/scExplorer
// Calculate omega_+(k) by finding zeros of 1 - lambda_+
func OmegaPlus(env *tempAll.Environment, k vec.Vector) (float64, error) {
	lp := lambdaPlusFn(env, k)
	initOmega, epsAbs, epsRel := 0.01, 1e-9, 1e-9
	root, err := solve.OneDimDiffRoot(lp, initOmega, epsAbs, epsRel)
	return root, err
}
コード例 #2
0
ファイル: omega.go プロジェクト: tflovorn/scExplorer
// Calculate omega_-(k) by finding zeros of 1 - lambda_-
func OmegaMinus(env *tempAll.Environment, k vec.Vector) (float64, error) {
	lm := lambdaMinusFn(env, k)
	root, err := solve.OneDimDiffRoot(lm, 0.01, 1e-9, 1e-9)
	return root, err
}
コード例 #3
0
ファイル: omegaPair.go プロジェクト: tflovorn/scExplorer
func OmegaPair(env *tempAll.Environment, k vec.Vector, r, s int) (float64, error) {
	lfn := LambdaFn(env, k, r, s)
	initOmega, epsAbs, epsRel := 0.01, 1e-9, 1e-9
	root, err := solve.OneDimDiffRoot(lfn, initOmega, epsAbs, epsRel)
	return root, err
}