Lines Matching refs:tl
30 TAO_NTL *tl = (TAO_NTL *)tao->data; in TaoSolve_NTL() local
66 tao->trust = PetscMax(tao->trust, tl->min_radius); in TaoSolve_NTL()
67 tao->trust = PetscMin(tao->trust, tl->max_radius); in TaoSolve_NTL()
74 tl->bfgs_pre = pc; in TaoSolve_NTL()
75 PetscCall(PCLMVMGetMatLMVM(tl->bfgs_pre, &tl->M)); in TaoSolve_NTL()
78 PetscCall(MatSetSizes(tl->M, n, n, N, N)); in TaoSolve_NTL()
79 PetscCall(MatLMVMAllocate(tl->M, tao->solution, tao->gradient)); in TaoSolve_NTL()
80 PetscCall(MatIsSymmetricKnown(tl->M, &sym_set, &is_symmetric)); in TaoSolve_NTL()
97 switch (tl->init_type) { in TaoSolve_NTL()
116 PetscCall(VecCopy(tao->solution, tl->W)); in TaoSolve_NTL()
117 PetscCall(VecAXPY(tl->W, -tao->trust / gnorm, tao->gradient)); in TaoSolve_NTL()
119 PetscCall(TaoComputeObjective(tao, tl->W, &ftrial)); in TaoSolve_NTL()
121 tau = tl->gamma1_i; in TaoSolve_NTL()
133 if ((PetscAbsScalar(actred) <= tl->epsilon) && (PetscAbsScalar(prered) <= tl->epsilon)) { in TaoSolve_NTL()
139 …tau_1 = tl->theta_i * gnorm * tao->trust / (tl->theta_i * gnorm * tao->trust + (1.0 - tl->theta_… in TaoSolve_NTL()
140 …tau_2 = tl->theta_i * gnorm * tao->trust / (tl->theta_i * gnorm * tao->trust - (1.0 + tl->theta_… in TaoSolve_NTL()
144 if (PetscAbsScalar(kappa - (PetscReal)1.0) <= tl->mu1_i) { in TaoSolve_NTL()
149 tau = tl->gamma3_i; in TaoSolve_NTL()
150 } else if (tau_max > tl->gamma4_i) { in TaoSolve_NTL()
151 tau = tl->gamma4_i; in TaoSolve_NTL()
152 } else if (tau_1 >= 1.0 && tau_1 <= tl->gamma4_i && tau_2 < 1.0) { in TaoSolve_NTL()
154 } else if (tau_2 >= 1.0 && tau_2 <= tl->gamma4_i && tau_1 < 1.0) { in TaoSolve_NTL()
159 } else if (PetscAbsScalar(kappa - (PetscReal)1.0) <= tl->mu2_i) { in TaoSolve_NTL()
163 if (tau_max < tl->gamma2_i) { in TaoSolve_NTL()
164 tau = tl->gamma2_i; in TaoSolve_NTL()
165 } else if (tau_max > tl->gamma3_i) { in TaoSolve_NTL()
166 tau = tl->gamma3_i; in TaoSolve_NTL()
173 tau = tl->gamma2_i; in TaoSolve_NTL()
174 } else if (tau_max < tl->gamma1_i) { in TaoSolve_NTL()
175 tau = tl->gamma1_i; in TaoSolve_NTL()
176 } else if ((tau_min < tl->gamma1_i) && (tau_max >= 1.0)) { in TaoSolve_NTL()
177 tau = tl->gamma1_i; in TaoSolve_NTL()
178 …} else if ((tau_1 >= tl->gamma1_i) && (tau_1 < 1.0) && ((tau_2 < tl->gamma1_i) || (tau_2 >= 1.0)))… in TaoSolve_NTL()
180 …} else if ((tau_2 >= tl->gamma1_i) && (tau_2 < 1.0) && ((tau_1 < tl->gamma1_i) || (tau_2 >= 1.0)))… in TaoSolve_NTL()
208 tao->trust = PetscMax(tao->trust, tl->min_radius); in TaoSolve_NTL()
209 tao->trust = PetscMin(tao->trust, tl->max_radius); in TaoSolve_NTL()
219 tl->ntrust = 0; in TaoSolve_NTL()
220 tl->newt = 0; in TaoSolve_NTL()
221 tl->bfgs = 0; in TaoSolve_NTL()
222 tl->grad = 0; in TaoSolve_NTL()
236 if (tl->bfgs_pre) { in TaoSolve_NTL()
238 PetscCall(MatLMVMUpdate(tl->M, tao->solution, tao->gradient)); in TaoSolve_NTL()
243 PetscCall(KSPCGSetRadius(tao->ksp, tl->max_radius)); in TaoSolve_NTL()
256 tao->trust = PetscMax(tao->trust, tl->min_radius); in TaoSolve_NTL()
257 tao->trust = PetscMin(tao->trust, tl->max_radius); in TaoSolve_NTL()
264 tao->trust = PetscMax(tao->trust, tl->min_radius); in TaoSolve_NTL()
265 tao->trust = PetscMin(tao->trust, tl->max_radius); in TaoSolve_NTL()
267 PetscCall(KSPCGSetRadius(tao->ksp, tl->max_radius)); in TaoSolve_NTL()
280 if ((KSP_DIVERGED_INDEFINITE_PC == ksp_reason) && (tl->bfgs_pre)) { in TaoSolve_NTL()
283 PetscCall(MatLMVMReset(tl->M, PETSC_FALSE)); in TaoSolve_NTL()
284 PetscCall(MatLMVMUpdate(tl->M, tao->solution, tao->gradient)); in TaoSolve_NTL()
290 if (NTL_UPDATE_REDUCTION == tl->update_type) { in TaoSolve_NTL()
297 tao->trust = tl->alpha1 * PetscMin(tao->trust, norm_d); in TaoSolve_NTL()
301 PetscCall(VecCopy(tao->solution, tl->W)); in TaoSolve_NTL()
302 PetscCall(VecAXPY(tl->W, 1.0, tao->stepdirection)); in TaoSolve_NTL()
303 PetscCall(TaoComputeObjective(tao, tl->W, &ftrial)); in TaoSolve_NTL()
306 tao->trust = tl->alpha1 * PetscMin(tao->trust, norm_d); in TaoSolve_NTL()
312 if ((PetscAbsScalar(actred) <= tl->epsilon) && (PetscAbsScalar(prered) <= tl->epsilon)) { in TaoSolve_NTL()
319 if (kappa < tl->eta1) { in TaoSolve_NTL()
321 tao->trust = tl->alpha1 * PetscMin(tao->trust, norm_d); in TaoSolve_NTL()
325 if (kappa < tl->eta2) { in TaoSolve_NTL()
327 tao->trust = tl->alpha2 * PetscMin(tao->trust, norm_d); in TaoSolve_NTL()
328 } else if (kappa < tl->eta3) { in TaoSolve_NTL()
330 tao->trust = tl->alpha3 * tao->trust; in TaoSolve_NTL()
331 } else if (kappa < tl->eta4) { in TaoSolve_NTL()
333 tao->trust = PetscMax(tl->alpha4 * norm_d, tao->trust); in TaoSolve_NTL()
336 tao->trust = PetscMax(tl->alpha5 * norm_d, tao->trust); in TaoSolve_NTL()
348 tao->trust = tl->gamma1 * PetscMin(tao->trust, norm_d); in TaoSolve_NTL()
351 PetscCall(VecCopy(tao->solution, tl->W)); in TaoSolve_NTL()
352 PetscCall(VecAXPY(tl->W, 1.0, tao->stepdirection)); in TaoSolve_NTL()
353 PetscCall(TaoComputeObjective(tao, tl->W, &ftrial)); in TaoSolve_NTL()
355 tao->trust = tl->gamma1 * PetscMin(tao->trust, norm_d); in TaoSolve_NTL()
362 if ((PetscAbsScalar(actred) <= tl->epsilon) && (PetscAbsScalar(prered) <= tl->epsilon)) { in TaoSolve_NTL()
368 tau_1 = tl->theta * gdx / (tl->theta * gdx - (1.0 - tl->theta) * prered + actred); in TaoSolve_NTL()
369 tau_2 = tl->theta * gdx / (tl->theta * gdx + (1.0 + tl->theta) * prered - actred); in TaoSolve_NTL()
373 if (kappa >= 1.0 - tl->mu1) { in TaoSolve_NTL()
376 tao->trust = PetscMax(tao->trust, tl->gamma3 * norm_d); in TaoSolve_NTL()
377 } else if (tau_max > tl->gamma4) { in TaoSolve_NTL()
378 tao->trust = PetscMax(tao->trust, tl->gamma4 * norm_d); in TaoSolve_NTL()
382 } else if (kappa >= 1.0 - tl->mu2) { in TaoSolve_NTL()
385 if (tau_max < tl->gamma2) { in TaoSolve_NTL()
386 tao->trust = tl->gamma2 * PetscMin(tao->trust, norm_d); in TaoSolve_NTL()
387 } else if (tau_max > tl->gamma3) { in TaoSolve_NTL()
388 tao->trust = PetscMax(tao->trust, tl->gamma3 * norm_d); in TaoSolve_NTL()
397 tao->trust = tl->gamma2 * PetscMin(tao->trust, norm_d); in TaoSolve_NTL()
398 } else if (tau_max < tl->gamma1) { in TaoSolve_NTL()
399 tao->trust = tl->gamma1 * PetscMin(tao->trust, norm_d); in TaoSolve_NTL()
400 } else if ((tau_min < tl->gamma1) && (tau_max >= 1.0)) { in TaoSolve_NTL()
401 tao->trust = tl->gamma1 * PetscMin(tao->trust, norm_d); in TaoSolve_NTL()
402 … } else if ((tau_1 >= tl->gamma1) && (tau_1 < 1.0) && ((tau_2 < tl->gamma1) || (tau_2 >= 1.0))) { in TaoSolve_NTL()
404 … } else if ((tau_2 >= tl->gamma1) && (tau_2 < 1.0) && ((tau_1 < tl->gamma1) || (tau_2 >= 1.0))) { in TaoSolve_NTL()
422 if (!tl->bfgs_pre) { in TaoSolve_NTL()
427 ++tl->grad; in TaoSolve_NTL()
431 PetscCall(MatSolve(tl->M, tao->gradient, tao->stepdirection)); in TaoSolve_NTL()
443 PetscCall(MatLMVMReset(tl->M, PETSC_FALSE)); in TaoSolve_NTL()
444 PetscCall(MatLMVMUpdate(tl->M, tao->solution, tao->gradient)); in TaoSolve_NTL()
445 PetscCall(MatSolve(tl->M, tao->gradient, tao->stepdirection)); in TaoSolve_NTL()
449 ++tl->grad; in TaoSolve_NTL()
452 PetscCall(MatLMVMGetUpdateCount(tl->M, &bfgsUpdates)); in TaoSolve_NTL()
455 ++tl->grad; in TaoSolve_NTL()
458 ++tl->bfgs; in TaoSolve_NTL()
465 ++tl->newt; in TaoSolve_NTL()
471 PetscCall(VecCopy(tao->solution, tl->Xold)); in TaoSolve_NTL()
472 PetscCall(VecCopy(tao->gradient, tl->Gold)); in TaoSolve_NTL()
481 PetscCall(VecCopy(tl->Xold, tao->solution)); in TaoSolve_NTL()
482 PetscCall(VecCopy(tl->Gold, tao->gradient)); in TaoSolve_NTL()
488 if (tl->bfgs_pre) { in TaoSolve_NTL()
492 ++tl->grad; in TaoSolve_NTL()
496 PetscCall(MatSolve(tl->M, tao->gradient, tao->stepdirection)); in TaoSolve_NTL()
504 PetscCall(MatLMVMReset(tl->M, PETSC_FALSE)); in TaoSolve_NTL()
505 PetscCall(MatLMVMUpdate(tl->M, tao->solution, tao->gradient)); in TaoSolve_NTL()
506 PetscCall(MatSolve(tl->M, tao->gradient, tao->stepdirection)); in TaoSolve_NTL()
509 ++tl->grad; in TaoSolve_NTL()
512 PetscCall(MatLMVMGetUpdateCount(tl->M, &bfgsUpdates)); in TaoSolve_NTL()
515 ++tl->grad; in TaoSolve_NTL()
518 ++tl->bfgs; in TaoSolve_NTL()
529 PetscCall(MatLMVMReset(tl->M, PETSC_FALSE)); in TaoSolve_NTL()
530 PetscCall(MatLMVMUpdate(tl->M, tao->solution, tao->gradient)); in TaoSolve_NTL()
531 PetscCall(MatSolve(tl->M, tao->gradient, tao->stepdirection)); in TaoSolve_NTL()
534 ++tl->grad; in TaoSolve_NTL()
550 PetscCall(VecCopy(tl->Xold, tao->solution)); in TaoSolve_NTL()
551 PetscCall(VecCopy(tl->Gold, tao->gradient)); in TaoSolve_NTL()
557 if (step < tl->nu1) { in TaoSolve_NTL()
559 tao->trust = tl->omega1 * PetscMin(norm_d, tao->trust); in TaoSolve_NTL()
560 } else if (step < tl->nu2) { in TaoSolve_NTL()
562 tao->trust = tl->omega2 * PetscMin(norm_d, tao->trust); in TaoSolve_NTL()
563 } else if (step < tl->nu3) { in TaoSolve_NTL()
565 if (tl->omega3 < 1.0) { in TaoSolve_NTL()
566 tao->trust = tl->omega3 * PetscMin(norm_d, tao->trust); in TaoSolve_NTL()
567 } else if (tl->omega3 > 1.0) { in TaoSolve_NTL()
568 tao->trust = PetscMax(tl->omega3 * norm_d, tao->trust); in TaoSolve_NTL()
570 } else if (step < tl->nu4) { in TaoSolve_NTL()
572 tao->trust = PetscMax(tl->omega4 * norm_d, tao->trust); in TaoSolve_NTL()
575 tao->trust = PetscMax(tl->omega5 * norm_d, tao->trust); in TaoSolve_NTL()
579 tao->trust = tl->omega1 * PetscMin(norm_d, tao->trust); in TaoSolve_NTL()
583 PetscCall(VecCopy(tl->W, tao->solution)); in TaoSolve_NTL()
586 ++tl->ntrust; in TaoSolve_NTL()
590 tao->trust = PetscMin(tao->trust, tl->max_radius); in TaoSolve_NTL()
606 TAO_NTL *tl = (TAO_NTL *)tao->data; in TaoSetUp_NTL() local
611 if (!tl->W) PetscCall(VecDuplicate(tao->solution, &tl->W)); in TaoSetUp_NTL()
612 if (!tl->Xold) PetscCall(VecDuplicate(tao->solution, &tl->Xold)); in TaoSetUp_NTL()
613 if (!tl->Gold) PetscCall(VecDuplicate(tao->solution, &tl->Gold)); in TaoSetUp_NTL()
614 tl->bfgs_pre = NULL; in TaoSetUp_NTL()
615 tl->M = NULL; in TaoSetUp_NTL()
621 TAO_NTL *tl = (TAO_NTL *)tao->data; in TaoDestroy_NTL() local
625 PetscCall(VecDestroy(&tl->W)); in TaoDestroy_NTL()
626 PetscCall(VecDestroy(&tl->Xold)); in TaoDestroy_NTL()
627 PetscCall(VecDestroy(&tl->Gold)); in TaoDestroy_NTL()
636 TAO_NTL *tl = (TAO_NTL *)tao->data; in TaoSetFromOptions_NTL() local
640 …adius initialization type", "", NTL_INIT, NTL_INIT_TYPES, NTL_INIT[tl->init_type], &tl->init_type,… in TaoSetFromOptions_NTL()
641 …radius update type", "", NTL_UPDATE, NTL_UPDATE_TYPES, NTL_UPDATE[tl->update_type], &tl->update_ty… in TaoSetFromOptions_NTL()
642 …PetscCall(PetscOptionsReal("-tao_ntl_eta1", "poor steplength; reduce radius", "", tl->eta1, &tl->e… in TaoSetFromOptions_NTL()
643 …Real("-tao_ntl_eta2", "reasonable steplength; leave radius alone", "", tl->eta2, &tl->eta2, NULL)); in TaoSetFromOptions_NTL()
644 …PetscCall(PetscOptionsReal("-tao_ntl_eta3", "good steplength; increase radius", "", tl->eta3, &tl-… in TaoSetFromOptions_NTL()
645 …("-tao_ntl_eta4", "excellent steplength; greatly increase radius", "", tl->eta4, &tl->eta4, NULL)); in TaoSetFromOptions_NTL()
646 PetscCall(PetscOptionsReal("-tao_ntl_alpha1", "", "", tl->alpha1, &tl->alpha1, NULL)); in TaoSetFromOptions_NTL()
647 PetscCall(PetscOptionsReal("-tao_ntl_alpha2", "", "", tl->alpha2, &tl->alpha2, NULL)); in TaoSetFromOptions_NTL()
648 PetscCall(PetscOptionsReal("-tao_ntl_alpha3", "", "", tl->alpha3, &tl->alpha3, NULL)); in TaoSetFromOptions_NTL()
649 PetscCall(PetscOptionsReal("-tao_ntl_alpha4", "", "", tl->alpha4, &tl->alpha4, NULL)); in TaoSetFromOptions_NTL()
650 PetscCall(PetscOptionsReal("-tao_ntl_alpha5", "", "", tl->alpha5, &tl->alpha5, NULL)); in TaoSetFromOptions_NTL()
651 …PetscCall(PetscOptionsReal("-tao_ntl_nu1", "poor steplength; reduce radius", "", tl->nu1, &tl->nu1… in TaoSetFromOptions_NTL()
652 …onsReal("-tao_ntl_nu2", "reasonable steplength; leave radius alone", "", tl->nu2, &tl->nu2, NULL)); in TaoSetFromOptions_NTL()
653 …PetscCall(PetscOptionsReal("-tao_ntl_nu3", "good steplength; increase radius", "", tl->nu3, &tl->n… in TaoSetFromOptions_NTL()
654 …eal("-tao_ntl_nu4", "excellent steplength; greatly increase radius", "", tl->nu4, &tl->nu4, NULL)); in TaoSetFromOptions_NTL()
655 PetscCall(PetscOptionsReal("-tao_ntl_omega1", "", "", tl->omega1, &tl->omega1, NULL)); in TaoSetFromOptions_NTL()
656 PetscCall(PetscOptionsReal("-tao_ntl_omega2", "", "", tl->omega2, &tl->omega2, NULL)); in TaoSetFromOptions_NTL()
657 PetscCall(PetscOptionsReal("-tao_ntl_omega3", "", "", tl->omega3, &tl->omega3, NULL)); in TaoSetFromOptions_NTL()
658 PetscCall(PetscOptionsReal("-tao_ntl_omega4", "", "", tl->omega4, &tl->omega4, NULL)); in TaoSetFromOptions_NTL()
659 PetscCall(PetscOptionsReal("-tao_ntl_omega5", "", "", tl->omega5, &tl->omega5, NULL)); in TaoSetFromOptions_NTL()
660 PetscCall(PetscOptionsReal("-tao_ntl_mu1_i", "", "", tl->mu1_i, &tl->mu1_i, NULL)); in TaoSetFromOptions_NTL()
661 PetscCall(PetscOptionsReal("-tao_ntl_mu2_i", "", "", tl->mu2_i, &tl->mu2_i, NULL)); in TaoSetFromOptions_NTL()
662 PetscCall(PetscOptionsReal("-tao_ntl_gamma1_i", "", "", tl->gamma1_i, &tl->gamma1_i, NULL)); in TaoSetFromOptions_NTL()
663 PetscCall(PetscOptionsReal("-tao_ntl_gamma2_i", "", "", tl->gamma2_i, &tl->gamma2_i, NULL)); in TaoSetFromOptions_NTL()
664 PetscCall(PetscOptionsReal("-tao_ntl_gamma3_i", "", "", tl->gamma3_i, &tl->gamma3_i, NULL)); in TaoSetFromOptions_NTL()
665 PetscCall(PetscOptionsReal("-tao_ntl_gamma4_i", "", "", tl->gamma4_i, &tl->gamma4_i, NULL)); in TaoSetFromOptions_NTL()
666 PetscCall(PetscOptionsReal("-tao_ntl_theta_i", "", "", tl->theta_i, &tl->theta_i, NULL)); in TaoSetFromOptions_NTL()
667 PetscCall(PetscOptionsReal("-tao_ntl_mu1", "", "", tl->mu1, &tl->mu1, NULL)); in TaoSetFromOptions_NTL()
668 PetscCall(PetscOptionsReal("-tao_ntl_mu2", "", "", tl->mu2, &tl->mu2, NULL)); in TaoSetFromOptions_NTL()
669 PetscCall(PetscOptionsReal("-tao_ntl_gamma1", "", "", tl->gamma1, &tl->gamma1, NULL)); in TaoSetFromOptions_NTL()
670 PetscCall(PetscOptionsReal("-tao_ntl_gamma2", "", "", tl->gamma2, &tl->gamma2, NULL)); in TaoSetFromOptions_NTL()
671 PetscCall(PetscOptionsReal("-tao_ntl_gamma3", "", "", tl->gamma3, &tl->gamma3, NULL)); in TaoSetFromOptions_NTL()
672 PetscCall(PetscOptionsReal("-tao_ntl_gamma4", "", "", tl->gamma4, &tl->gamma4, NULL)); in TaoSetFromOptions_NTL()
673 PetscCall(PetscOptionsReal("-tao_ntl_theta", "", "", tl->theta, &tl->theta, NULL)); in TaoSetFromOptions_NTL()
674 …nsReal("-tao_ntl_min_radius", "lower bound on initial radius", "", tl->min_radius, &tl->min_radius… in TaoSetFromOptions_NTL()
675 …tscOptionsReal("-tao_ntl_max_radius", "upper bound on radius", "", tl->max_radius, &tl->max_radius… in TaoSetFromOptions_NTL()
676 …"tolerance used when computing actual and predicted reduction", "", tl->epsilon, &tl->epsilon, NUL… in TaoSetFromOptions_NTL()
685 TAO_NTL *tl = (TAO_NTL *)tao->data; in TaoView_NTL() local
692 … PetscCall(PetscViewerASCIIPrintf(viewer, "Trust-region steps: %" PetscInt_FMT "\n", tl->ntrust)); in TaoView_NTL()
693 PetscCall(PetscViewerASCIIPrintf(viewer, "Newton search steps: %" PetscInt_FMT "\n", tl->newt)); in TaoView_NTL()
694 PetscCall(PetscViewerASCIIPrintf(viewer, "BFGS search steps: %" PetscInt_FMT "\n", tl->bfgs)); in TaoView_NTL()
695 … PetscCall(PetscViewerASCIIPrintf(viewer, "Gradient search steps: %" PetscInt_FMT "\n", tl->grad)); in TaoView_NTL()
742 TAO_NTL *tl; in TaoCreate_NTL() local
746 PetscCall(PetscNew(&tl)); in TaoCreate_NTL()
758 tao->data = (void *)tl; in TaoCreate_NTL()
761 tl->nu1 = 0.25; in TaoCreate_NTL()
762 tl->nu2 = 0.50; in TaoCreate_NTL()
763 tl->nu3 = 1.00; in TaoCreate_NTL()
764 tl->nu4 = 1.25; in TaoCreate_NTL()
766 tl->omega1 = 0.25; in TaoCreate_NTL()
767 tl->omega2 = 0.50; in TaoCreate_NTL()
768 tl->omega3 = 1.00; in TaoCreate_NTL()
769 tl->omega4 = 2.00; in TaoCreate_NTL()
770 tl->omega5 = 4.00; in TaoCreate_NTL()
773 tl->eta1 = 1.0e-4; in TaoCreate_NTL()
774 tl->eta2 = 0.25; in TaoCreate_NTL()
775 tl->eta3 = 0.50; in TaoCreate_NTL()
776 tl->eta4 = 0.90; in TaoCreate_NTL()
778 tl->alpha1 = 0.25; in TaoCreate_NTL()
779 tl->alpha2 = 0.50; in TaoCreate_NTL()
780 tl->alpha3 = 1.00; in TaoCreate_NTL()
781 tl->alpha4 = 2.00; in TaoCreate_NTL()
782 tl->alpha5 = 4.00; in TaoCreate_NTL()
785 tl->mu1 = 0.10; in TaoCreate_NTL()
786 tl->mu2 = 0.50; in TaoCreate_NTL()
788 tl->gamma1 = 0.25; in TaoCreate_NTL()
789 tl->gamma2 = 0.50; in TaoCreate_NTL()
790 tl->gamma3 = 2.00; in TaoCreate_NTL()
791 tl->gamma4 = 4.00; in TaoCreate_NTL()
793 tl->theta = 0.05; in TaoCreate_NTL()
796 tl->mu1_i = 0.35; in TaoCreate_NTL()
797 tl->mu2_i = 0.50; in TaoCreate_NTL()
799 tl->gamma1_i = 0.0625; in TaoCreate_NTL()
800 tl->gamma2_i = 0.5; in TaoCreate_NTL()
801 tl->gamma3_i = 2.0; in TaoCreate_NTL()
802 tl->gamma4_i = 5.0; in TaoCreate_NTL()
804 tl->theta_i = 0.25; in TaoCreate_NTL()
807 tl->min_radius = 1.0e-10; in TaoCreate_NTL()
808 tl->max_radius = 1.0e10; in TaoCreate_NTL()
809 tl->epsilon = 1.0e-6; in TaoCreate_NTL()
811 tl->init_type = NTL_INIT_INTERPOLATION; in TaoCreate_NTL()
812 tl->update_type = NTL_UPDATE_REDUCTION; in TaoCreate_NTL()