Home
last modified time | relevance | path

Searched refs:NOBSERVATIONS (Results 1 – 4 of 4) sorted by relevance

/petsc/src/binding/petsc4py/demo/legacy/taosolve/
H A Dchwirut.py16 NOBSERVATIONS = 100
20 x = np.random.rand(NOBSERVATIONS)
21 e = np.random.rand(NOBSERVATIONS)
25 self.NOBSERVATIONS = NOBSERVATIONS
34 F.setSizes(self.NOBSERVATIONS)
/petsc/src/tao/leastsquares/tests/
H A Dchwirut2.c27 #define NOBSERVATIONS 214 macro
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 */ in main()
58 PetscInt w_col[NOBSERVATIONS]; in main()
59 PetscReal w_vals[NOBSERVATIONS]; in main()
69 PetscCall(VecCreateSeq(MPI_COMM_SELF, NOBSERVATIONS, &f)); in main()
75 for (i = 0; i < NOBSERVATIONS; i++) { in main()
[all …]
/petsc/src/tao/leastsquares/tutorials/
H A Dchwirut1.c26 #define NOBSERVATIONS 214 macro
32 PetscReal t[NOBSERVATIONS]; /* array of independent variables of observation */
33 PetscReal y[NOBSERVATIONS]; /* array of dependent variables */
34 PetscReal j[NOBSERVATIONS][NPARAMETERS]; /* dense jacobian matrix array*/
35 PetscInt idm[NOBSERVATIONS]; /* Matrix indices for jacobian */
60 PetscCall(VecCreateSeq(MPI_COMM_SELF, NOBSERVATIONS, &f)); in main()
63 PetscCall(MatCreateSeqDense(MPI_COMM_SELF, NOBSERVATIONS, NPARAMETERS, NULL, &J)); in main()
65 for (i = 0; i < NOBSERVATIONS; i++) user.idm[i] = i; in main()
114 …for (i = 0; i < NOBSERVATIONS; i++) f[i] = y[i] - PetscExpScalar(-x[0] * t[i]) / (x[1] + x[2] * t[… in EvaluateFunction()
117 PetscCall(PetscLogFlops(6 * NOBSERVATIONS)); in EvaluateFunction()
[all …]
H A Dchwirut2.c26 #define NOBSERVATIONS 214 macro
35 PetscReal t[NOBSERVATIONS]; /* array of independent variables of observation */
36 PetscReal y[NOBSERVATIONS]; /* array of dependent variables */
66 PetscCall(VecCreateSeq(PETSC_COMM_SELF, NOBSERVATIONS, &f)); in main()
111 for (i = 0; i < NOBSERVATIONS; i++) PetscCall(RunSimulation(x, i, &f[i], user)); in EvaluateFunction()
123 while (finishedtasks < NOBSERVATIONS || checkedin < user->size - 1) { in EvaluateFunction()
133 if (next_task < NOBSERVATIONS) { in EvaluateFunction()
145 PetscCall(PetscLogFlops(6 * NOBSERVATIONS)); in EvaluateFunction()