leastsquares/tests/chwirut2.c

27 #define NOBSERVATIONS 214
33   PetscReal t[NOBSERVATIONS];              /* array of independent variables of observation */
34   PetscReal y[NOBSERVATIONS];              /* array of dependent variables */
35   PetscReal j[NOBSERVATIONS][NPARAMETERS]; /* dense jacobian matrix array*/
36   PetscInt  idm[NOBSERVATIONS];            /* Matrix indices for jacobian */
57   PetscInt  w_row[NOBSERVATIONS]; /* explicit weights */
58   PetscInt  w_col[NOBSERVATIONS];
59   PetscReal w_vals[NOBSERVATIONS];
69   PetscCall(VecCreateSeq(MPI_COMM_SELF, NOBSERVATIONS, &f));
75   for (i = 0; i < NOBSERVATIONS; i++) {
82   PetscCall(MatCreateSeqDense(MPI_COMM_SELF, NOBSERVATIONS, NPARAMETERS, NULL, &J));
84   for (i = 0; i < NOBSERVATIONS; i++) user.idm[i] = i;
100     PetscCall(TaoSetResidualWeights(tao, NULL, NOBSERVATIONS, w_row, w_col, w_vals));
137   for (i = 0; i < NOBSERVATIONS; i++) f[i] = y[i] - PetscExpScalar(-x[0] * t[i]) / (x[1] + x[2] * t[i]);
140   PetscCall(PetscLogFlops(6 * NOBSERVATIONS));
156   for (i = 0; i < NOBSERVATIONS; i++) {
165   PetscCall(MatSetValues(J, NOBSERVATIONS, user->idm, NPARAMETERS, user->idn, (PetscReal *)user->j, INSERT_VALUES));
170   PetscCall(PetscLogFlops(NOBSERVATIONS * 13));