1 #ifndef __TAO_TRON_H 2 #define __TAO_TRON_H 3 4 #include <petsc/private/taoimpl.h> 5 #include <petscis.h> 6 7 typedef struct { 8 9 /* Parameters */ 10 PetscReal pg_ftol; 11 PetscReal actred; 12 PetscReal f_new; 13 14 PetscReal eta1,eta2,eta3,eta4; 15 PetscReal sigma1,sigma2,sigma3; 16 17 PetscInt maxgpits; 18 19 /* Problem variables, vectors and index sets */ 20 PetscReal stepsize; 21 PetscReal pgstepsize; 22 23 /* Problem statistics */ 24 25 PetscInt n; /* Dimension of the Problem */ 26 PetscReal delta; /* Trust region size */ 27 PetscReal gnorm; 28 PetscReal f; 29 30 PetscInt total_gp_its; 31 PetscInt gp_iterates; 32 33 34 Vec X_New; 35 Vec G_New; 36 Vec Work; 37 38 /* Subvectors and submatrices */ 39 Vec DXFree; 40 Vec R; 41 Vec rmask; 42 Vec diag; 43 Mat H_sub; 44 Mat Hpre_sub; 45 46 IS Free_Local; /* Indices of local variables equal to lower bound */ 47 VecScatter scatter; 48 49 PetscInt n_free; /* Number of free variables */ 50 PetscInt n_free_last; 51 52 } TAO_TRON; 53 54 #endif 55 56