// 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
}
// 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
}
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
}
