Lines Matching refs:user
77 AppCtx user; /* user-defined application context */ in main() local
120 user.Nvglobal = 16; /* Global # of vertices */ in main()
121 user.Neglobal = 18; /* Global # of elements */ in main()
123 PetscCall(PetscOptionsGetInt(NULL, NULL, "-vert", &user.Nvglobal, NULL)); in main()
124 PetscCall(PetscOptionsGetInt(NULL, NULL, "-elem", &user.Neglobal, NULL)); in main()
126 user.non_lin_param = 0.06; in main()
128 PetscCall(PetscOptionsGetReal(NULL, NULL, "-nl_par", &user.non_lin_param, NULL)); in main()
130 user.lin_param = -1.0; in main()
132 PetscCall(PetscOptionsGetReal(NULL, NULL, "-lin_par", &user.lin_param, NULL)); in main()
134 user.Nvlocal = 0; in main()
135 user.Nelocal = 0; in main()
166 PetscCall(PetscMalloc1(user.Nvglobal, &user.gloInd)); in main()
168 for (inode = 0; inode < user.Nvglobal; inode++) { in main()
171 user.v2p[inode] = dtmp; in main()
172 if (user.v2p[inode] == rank) { in main()
173 …fptr1, "Node %" PetscInt_FMT " belongs to processor %" PetscInt_FMT "\n", inode, user.v2p[inode])); in main()
175 user.gloInd[user.Nvlocal] = inode; in main()
180 user.itot[user.Nvlocal] = nbrs; in main()
182 for (i = 0; i < user.itot[user.Nvlocal]; i++) { in main()
188 user.AdjM[user.Nvlocal][i] = dtmp; in main()
190 … PetscCall(PetscFPrintf(PETSC_COMM_SELF, fptr1, "%" PetscInt_FMT " ", user.AdjM[user.Nvlocal][i])); in main()
193 user.Nvlocal++; in main()
196 …Printf(PETSC_COMM_SELF, fptr1, "Total # of Local Vertices is %" PetscInt_FMT " \n", user.Nvlocal)); in main()
208 PetscCallMPI(MPI_Scan(&user.Nvlocal, &rstart, 1, MPIU_INT, MPI_SUM, PETSC_COMM_WORLD)); in main()
209 rstart -= user.Nvlocal; in main()
210 PetscCall(PetscMalloc1(user.Nvlocal, &pordering)); in main()
212 for (i = 0; i < user.Nvlocal; i++) pordering[i] = rstart + i; in main()
217 PetscCall(AOCreateBasic(MPI_COMM_WORLD, user.Nvlocal, user.gloInd, pordering, &ao)); in main()
223 PetscCall(PetscMalloc1(user.Nvlocal, &user.locInd)); in main()
231 for (i = 0; i < user.Nvlocal; i++) { in main()
232 PetscCall(PetscFPrintf(PETSC_COMM_SELF, fptr1, " %" PetscInt_FMT " ", user.gloInd[i])); in main()
234 user.locInd[i] = user.gloInd[i]; in main()
238 for (i = 0; i < user.Nvlocal; i++) { in main()
239 …intf(PETSC_COMM_SELF, fptr1, "Neghbors of local vertex %" PetscInt_FMT " are : ", user.gloInd[i])); in main()
240 for (j = 0; j < user.itot[i]; j++) { in main()
241 PetscCall(PetscFPrintf(PETSC_COMM_SELF, fptr1, "%" PetscInt_FMT " ", user.AdjM[i][j])); in main()
243 tmp[j + jstart] = user.AdjM[i][j]; in main()
245 jstart += user.itot[i]; in main()
252 PetscCall(AOApplicationToPetsc(ao, user.Nvlocal, user.locInd)); in main()
257 …for (i = 0; i < user.Nvlocal; i++) PetscCall(PetscFPrintf(PETSC_COMM_SELF, fptr1, " %" PetscInt_FM… in main()
261 for (i = 0; i < user.Nvlocal; i++) { in main()
262 …intf(PETSC_COMM_SELF, fptr1, "Neghbors of local vertex %" PetscInt_FMT " are : ", user.locInd[i])); in main()
263 for (j = 0; j < user.itot[i]; j++) { in main()
264 user.AdjM[i][j] = tmp[j + jstart]; in main()
266 PetscCall(PetscFPrintf(PETSC_COMM_SELF, fptr1, "%" PetscInt_FMT " ", user.AdjM[i][j])); in main()
268 jstart += user.itot[i]; in main()
291 PetscCall(PetscMalloc1(user.Nvglobal, &vertices)); in main()
292 PetscCall(PetscCalloc1(user.Nvglobal, &verticesmask)); in main()
298 for (i = 0; i < user.Nvlocal; i++) { in main()
299 vertices[nvertices++] = user.locInd[i]; in main()
300 verticesmask[user.locInd[i]] = nvertices; in main()
306 for (i = 0; i < user.Nvlocal; i++) { in main()
307 for (j = 0; j < user.itot[i]; j++) { in main()
308 nb = user.AdjM[i][j]; in main()
327 for (i = 0; i < user.Nvlocal; i++) { in main()
329 for (j = 0; j < user.itot[i]; j++) { in main()
330 nb = user.AdjM[i][j]; in main()
331 user.AdjM[i][j] = verticesmask[nb] - 1; in main()
333 PetscCall(PetscFPrintf(PETSC_COMM_SELF, fptr1, "%" PetscInt_FMT " ", user.AdjM[i][j])); in main()
338 N = user.Nvglobal; in main()
348 PetscCall(VecSetSizes(x, user.Nvlocal, N)); in main()
351 PetscCall(VecCreateSeq(MPI_COMM_SELF, bs * nvertices, &user.localX)); in main()
352 PetscCall(VecDuplicate(user.localX, &user.localF)); in main()
360 PetscCall(VecScatterCreate(x, isglobal, user.localX, islocal, &user.scatter)); in main()
390 PetscCall(SNESSetFunction(snes, r, FormFunction, (void *)&user)); in main()
394 PetscCall(SNESSetJacobian(snes, Jac, Jac, FormJacobian, (void *)&user)); in main()
400 PetscCall(FormJacobian(snes, x, Jac, Jac, &user)); in main()
408 PetscCall(MatFDColoringSetFunction(matfdcoloring, (MatFDColoringFn *)FormFunction, &user)); in main()
432 PetscCall(FormInitialGuess(&user, x)); in main()
438 …for (inode = 0; inode < user.Nvlocal; inode++) PetscCall(PetscFPrintf(PETSC_COMM_SELF, fptr1, "Ini… in main()
456 …for (inode = 0; inode < user.Nvlocal; inode++) PetscCall(PetscFPrintf(PETSC_COMM_SELF, fptr1, "Sol… in main()
466 PetscCall(PetscFree(user.gloInd)); in main()
467 PetscCall(PetscFree(user.locInd)); in main()
471 PetscCall(VecScatterDestroy(&user.scatter)); in main()
475 PetscCall(VecDestroy(&user.localX)); in main()
476 PetscCall(VecDestroy(&user.localF)); in main()
496 PetscErrorCode FormInitialGuess(AppCtx *user, Vec X) in FormInitialGuess() argument
506 Nvglobal = user->Nvglobal; in FormInitialGuess()
507 Neglobal = user->Neglobal; in FormInitialGuess()
508 lambda = user->non_lin_param; in FormInitialGuess()
509 alpha = user->lin_param; in FormInitialGuess()
513 Nvlocal = user->Nvlocal; in FormInitialGuess()
514 gloInd = user->gloInd; in FormInitialGuess()
550 AppCtx *user = (AppCtx *)ptr; in FormFunction() local
555 Vec localX = user->localX; in FormFunction()
563 Nvglobal = user->Nvglobal; in FormFunction()
564 Neglobal = user->Neglobal; in FormFunction()
565 gloInd = user->gloInd; in FormFunction()
569 Nvlocal = user->Nvlocal; in FormFunction()
570 lambda = user->non_lin_param; in FormFunction()
571 alpha = user->lin_param; in FormFunction()
572 scatter = user->scatter; in FormFunction()
599 f[i] = (user->itot[i] - alpha) * x[i] - lambda * x[i] * x[i]; in FormFunction()
600 for (PetscInt j = 0; j < user->itot[i]; j++) f[i] -= x[user->AdjM[i][j]]; in FormFunction()
628 AppCtx *user = (AppCtx *)ptr; in FormJacobian() local
631 Vec localX = user->localX; in FormJacobian()
639 Nvglobal = user->Nvglobal; in FormJacobian()
640 Neglobal = user->Neglobal; in FormJacobian()
641 gloInd = user->gloInd; in FormJacobian()
646 Nvlocal = user->Nvlocal; in FormJacobian()
647 lambda = user->non_lin_param; in FormJacobian()
648 alpha = user->lin_param; in FormJacobian()
649 scatter = user->scatter; in FormJacobian()
670 value[0] = user->itot[i] - 2.0 * lambda * x[i] - alpha; in FormJacobian()
671 for (PetscInt j = 0; j < user->itot[i]; j++) { in FormJacobian()
672 col[j + 1] = user->AdjM[i][j]; in FormJacobian()
681 PetscCall(MatSetValuesLocal(jac, 1, &i, 1 + user->itot[i], col, value, INSERT_VALUES)); in FormJacobian()