1*a4963045SJacob Faibussowitsch #pragma once 2af0996ceSBarry Smith #include <petsc/private/taoimpl.h> 3a7e14dcfSSatish Balay 4a7e14dcfSSatish Balay typedef struct { 5a7e14dcfSSatish Balay PetscReal mu_ic; 6a7e14dcfSSatish Balay PetscReal mu_oc; 7a7e14dcfSSatish Balay PetscReal mu_r; 8a7e14dcfSSatish Balay PetscReal mu_e; 9a7e14dcfSSatish Balay 10a8d3b578SPierre Jolivet PetscReal lambda; /* starting point delta for finding starting simplex */ 11a7e14dcfSSatish Balay 12a7e14dcfSSatish Balay PetscInt N; 13a7e14dcfSSatish Balay PetscReal oneOverN; 14a7e14dcfSSatish Balay Vec Xbar, Xmuc, Xmur, Xmue; 15a7e14dcfSSatish Balay Vec G; 16a7e14dcfSSatish Balay Vec *simplex; 17a7e14dcfSSatish Balay 18a7e14dcfSSatish Balay PetscReal *f_values; 19a7e14dcfSSatish Balay PetscInt *indices; 20a7e14dcfSSatish Balay 21a7e14dcfSSatish Balay PetscInt nshrink; 22a7e14dcfSSatish Balay PetscInt nexpand; 23a7e14dcfSSatish Balay PetscInt nreflect; 24a7e14dcfSSatish Balay PetscInt nincontract; 25a7e14dcfSSatish Balay PetscInt noutcontract; 26a7e14dcfSSatish Balay 27a7e14dcfSSatish Balay } TAO_NelderMead; 28