210   PetscInt         m, n;  in EvaluateFunction()  local
219   for (m = 0; m < M; m++) {  in EvaluateFunction()
220     f[m] = -b[m];  in EvaluateFunction()
221     for (n = 0; n < N; n++) f[m] += user->A[m][n] * x[n];  in EvaluateFunction()
234   PetscInt         m, n;  in EvaluateJacobian()  local
241   for (m = 0; m < M; ++m) {  in EvaluateJacobian()
242     for (n = 0; n < N; ++n) user->J[m][n] = user->A[m][n];  in EvaluateJacobian()
279   PetscInt   m, n, k;     /* loop index for M,N,K dimension. */  in InitializeUserData()  local
283   m      = 0;  in InitializeUserData()
284   b[m++] = 0.28;  in InitializeUserData()
285   b[m++] = 0.55;  in InitializeUserData()
286   b[m++] = 0.96;  in InitializeUserData()
293   m               = 0;  in InitializeUserData()
295   user->A[m][n++] = 0.81;  in InitializeUserData()
296   user->A[m][n++] = 0.91;  in InitializeUserData()
297   user->A[m][n++] = 0.28;  in InitializeUserData()
298   user->A[m][n++] = 0.96;  in InitializeUserData()
299   user->A[m][n++] = 0.96;  in InitializeUserData()
300   ++m;  in InitializeUserData()
302   user->A[m][n++] = 0.91;  in InitializeUserData()
303   user->A[m][n++] = 0.63;  in InitializeUserData()
304   user->A[m][n++] = 0.55;  in InitializeUserData()
305   user->A[m][n++] = 0.16;  in InitializeUserData()
306   user->A[m][n++] = 0.49;  in InitializeUserData()
307   ++m;  in InitializeUserData()
309   user->A[m][n++] = 0.13;  in InitializeUserData()
310   user->A[m][n++] = 0.10;  in InitializeUserData()
311   user->A[m][n++] = 0.96;  in InitializeUserData()
312   user->A[m][n++] = 0.97;  in InitializeUserData()
313   user->A[m][n++] = 0.80;  in InitializeUserData()