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