| /petsc/src/ts/tutorials/ |
| H A D | ex20opt_ic.c | 209 const PetscScalar *x_ptr; in FormFunctionGradient() local 221 PetscCall(VecGetArrayRead(user_ptr->U, &x_ptr)); in FormFunctionGradient() 223 …*f = (x_ptr[0] - user_ptr->ob[0]) * (x_ptr[0] - user_ptr->ob[0]) + (x_ptr[1] - user_ptr->ob[… in FormFunctionGradient() 224 y_ptr[0] = 2. * (x_ptr[0] - user_ptr->ob[0]); in FormFunctionGradient() 225 y_ptr[1] = 2. * (x_ptr[1] - user_ptr->ob[1]); in FormFunctionGradient() 227 PetscCall(VecRestoreArrayRead(user_ptr->U, &x_ptr)); in FormFunctionGradient() 239 PetscScalar *x_ptr; in FormHessian() local 245 PetscCall(VecGetArray(user_ptr->Dir, &x_ptr)); in FormHessian() 246 x_ptr[0] = 1.; in FormHessian() 247 x_ptr[1] = 0.; in FormHessian() [all …]
|
| H A D | ex23fwdadj.c | 93 PetscScalar *x_ptr; in main() local 127 PetscCall(VecGetArrayWrite(user.x, &x_ptr)); in main() 128 x_ptr[0] = user.a; in main() 129 PetscCall(VecRestoreArrayWrite(user.x, &x_ptr)); in main() 148 PetscCall(VecGetArray(user.x, &x_ptr)); in main() 149 PetscCall(PetscPrintf(PETSC_COMM_WORLD, "\n ode solution %g\n", (double)PetscRealPart(x_ptr[0]))); in main() 150 PetscCall(VecRestoreArray(user.x, &x_ptr)); in main() 160 PetscCall(VecGetArrayWrite(user.lambda[0], &x_ptr)); in main() 161 x_ptr[0] = 1.0; in main() 162 PetscCall(VecRestoreArrayWrite(user.lambda[0], &x_ptr)); in main() [all …]
|
| H A D | ex20opt_p.c | 367 PetscScalar *x_ptr; in main() local 439 PetscCall(VecGetArray(user.U, &x_ptr)); in main() 440 x_ptr[0] = 2.0; in main() 441 x_ptr[1] = -2.0 / 3.0 + 10.0 / (81.0 * user.mu) - 292.0 / (2187.0 * user.mu * user.mu); in main() 442 PetscCall(VecRestoreArray(user.U, &x_ptr)); in main() 466 PetscCall(VecGetArray(P, &x_ptr)); in main() 467 x_ptr[0] = PetscRealConstant(1.2); in main() 468 PetscCall(VecRestoreArray(P, &x_ptr)); in main() 529 PetscScalar *x_ptr, *g; in FormFunctionGradient() local 541 PetscCall(VecGetArray(user_ptr->U, &x_ptr)); in FormFunctionGradient() [all …]
|
| H A D | ex20adj.c | 241 PetscScalar *x_ptr, *y_ptr, derp; in main() local 316 PetscCall(VecGetArray(user.U, &x_ptr)); in main() 317 x_ptr[0] = 2.0; in main() 318 x_ptr[1] = -2.0 / 3.0 + 10.0 / (81.0 * user.mu) - 292.0 / (2187.0 * user.mu * user.mu); in main() 319 PetscCall(VecRestoreArray(user.U, &x_ptr)); in main() 352 PetscCall(VecGetArray(user.mup[0], &x_ptr)); in main() 353 x_ptr[0] = 0.0; in main() 354 PetscCall(VecRestoreArray(user.mup[0], &x_ptr)); in main() 356 PetscCall(VecGetArray(user.mup[1], &x_ptr)); in main() 357 x_ptr[0] = 0.0; in main() [all …]
|
| H A D | ex49.c | 78 PetscScalar *x_ptr; in main() local 123 PetscCall(VecGetArray(x, &x_ptr)); in main() 124 x_ptr[0] = 2.0; in main() 125 x_ptr[1] = -6.6e-01; in main() 126 PetscCall(VecRestoreArray(x, &x_ptr)); in main()
|
| H A D | ex19.c | 160 PetscScalar *x_ptr; in main() local 204 PetscCall(VecGetArray(x, &x_ptr)); in main() 205 x_ptr[0] = -2; in main() 206 x_ptr[1] = -2.355301397608119909925287735864250951918; in main() 207 PetscCall(VecRestoreArray(x, &x_ptr)); in main()
|
| H A D | ex16fwd.c | 144 PetscScalar *x_ptr; in main() local 201 PetscCall(VecGetArray(x, &x_ptr)); in main() 203 x_ptr[0] = 2; in main() 204 x_ptr[1] = 0.66666654321; in main() 205 PetscCall(VecRestoreArray(x, &x_ptr)); in main()
|
| H A D | ex20.c | 122 PetscScalar *x_ptr; in main() local 175 PetscCall(VecGetArray(x, &x_ptr)); in main() 176 x_ptr[0] = 2.0; in main() 177 x_ptr[1] = -2.0 / 3.0 + 10.0 / (81.0 * user.mu) - 292.0 / (2187.0 * user.mu * user.mu); in main() 178 PetscCall(VecRestoreArray(x, &x_ptr)); in main()
|
| H A D | ex16.c | 188 PetscScalar *x_ptr; in main() local 237 PetscCall(VecGetArray(x, &x_ptr)); in main() 238 x_ptr[0] = 2.0; in main() 239 x_ptr[1] = -2.0 / 3.0 + 10.0 / (81.0 * user.mu) - 292.0 / (2187.0 * user.mu * user.mu); in main() 240 PetscCall(VecRestoreArray(x, &x_ptr)); in main()
|
| H A D | ex20fwd.c | 126 PetscScalar *x_ptr; in main() local 177 PetscCall(VecGetArray(user.x, &x_ptr)); in main() 178 x_ptr[0] = 2.0; in main() 179 x_ptr[1] = -0.66666654321; in main() 180 PetscCall(VecRestoreArray(user.x, &x_ptr)); in main()
|
| H A D | ex20td.c | 244 PetscScalar *x_ptr, *y_ptr, *u_ptr; in main() local 327 PetscCall(VecGetArray(user.U, &x_ptr)); in main() 328 x_ptr[0] = 2.0; in main() 329 x_ptr[1] = -2.0 / 3.0; in main() 330 PetscCall(VecRestoreArray(user.U, &x_ptr)); in main()
|
| /petsc/src/ts/tutorials/autodiff/ |
| H A D | ex16opt_ic.cxx | 128 PetscScalar *x_ptr; in main() local 172 PetscCall(VecGetArray(user.x, &x_ptr)); in main() 173 x_ptr[0] = 2.0; in main() 174 x_ptr[1] = 0.66666654321; in main() 175 PetscCall(VecRestoreArray(user.x, &x_ptr)); in main() 218 PetscCall(VecGetArray(user.x, &x_ptr)); in main() 219 user.x_ob[0] = x_ptr[0]; in main() 220 user.x_ob[1] = x_ptr[1]; in main() 221 PetscCall(VecRestoreArray(user.x, &x_ptr)); in main() 231 PetscCall(VecGetArray(ic, &x_ptr)); in main() [all …]
|
| H A D | ex16adj.cxx | 175 PetscScalar *x_ptr; in main() local 228 PetscCall(VecGetArray(x, &x_ptr)); in main() 229 x_ptr[0] = 2; in main() 230 x_ptr[1] = 0.66666654321; in main() 231 PetscCall(VecRestoreArray(x, &x_ptr)); in main() 277 PetscCall(VecGetArray(lambda[0], &x_ptr)); in main() 278 x_ptr[0] = 1.0; in main() 279 x_ptr[1] = 0.0; in main() 280 PetscCall(VecRestoreArray(lambda[0], &x_ptr)); in main() 281 PetscCall(VecGetArray(lambda[1], &x_ptr)); in main() [all …]
|
| H A D | ex16adj_tl.cxx | 195 PetscScalar *x_ptr; in main() local 248 PetscCall(VecGetArray(x, &x_ptr)); in main() 249 x_ptr[0] = 2; in main() 250 x_ptr[1] = 0.66666654321; in main() 251 PetscCall(VecRestoreArray(x, &x_ptr)); in main() 287 PetscCall(VecGetArray(lambda[0], &x_ptr)); in main() 288 x_ptr[0] = 1.0; in main() 289 x_ptr[1] = 0.0; in main() 290 PetscCall(VecRestoreArray(lambda[0], &x_ptr)); in main() 291 PetscCall(VecGetArray(lambda[1], &x_ptr)); in main() [all …]
|
| /petsc/src/ts/tutorials/power_grid/ |
| H A D | ex3opt.c | 43 PetscScalar *x_ptr; in main() local 167 PetscCall(VecGetArray(p, &x_ptr)); in main() 168 x_ptr[0] = ctx.Pm; in main() 169 PetscCall(VecRestoreArray(p, &x_ptr)); in main() 178 PetscCall(VecGetArray(lowerb, &x_ptr)); in main() 179 x_ptr[0] = 0.; in main() 180 PetscCall(VecRestoreArray(lowerb, &x_ptr)); in main() 181 PetscCall(VecGetArray(upperb, &x_ptr)); in main() 182 x_ptr[0] = 1.1; in main() 183 PetscCall(VecRestoreArray(upperb, &x_ptr)); in main() [all …]
|
| H A D | ex3opt_fd.c | 48 PetscScalar *x_ptr; in main() local 105 PetscCall(VecGetArray(p, &x_ptr)); in main() 106 x_ptr[0] = ctx.Pm; in main() 107 PetscCall(VecRestoreArray(p, &x_ptr)); in main() 117 PetscCall(VecGetArray(lowerb, &x_ptr)); in main() 118 x_ptr[0] = 0.; in main() 119 PetscCall(VecRestoreArray(lowerb, &x_ptr)); in main() 120 PetscCall(VecGetArray(upperb, &x_ptr)); in main() 121 x_ptr[0] = 1.1; in main() 122 PetscCall(VecRestoreArray(upperb, &x_ptr)); in main() [all …]
|
| H A D | ex9opt.c | 170 PetscScalar *x_ptr; in main() local 280 PetscCall(VecGetArray(p, &x_ptr)); in main() 281 x_ptr[0] = ctx.Pm; in main() 282 PetscCall(VecRestoreArray(p, &x_ptr)); in main() 292 PetscCall(VecGetArray(lowerb, &x_ptr)); in main() 293 x_ptr[0] = 0.; in main() 294 PetscCall(VecRestoreArray(lowerb, &x_ptr)); in main() 295 PetscCall(VecGetArray(upperb, &x_ptr)); in main() 296 x_ptr[0] = PetscRealConstant(1.1); in main() 297 PetscCall(VecRestoreArray(upperb, &x_ptr)); in main() [all …]
|
| H A D | ex3sa.c | 36 PetscScalar *x_ptr, *y_ptr, *s_ptr; in main() local 209 PetscCall(VecGetArray(mu[0], &x_ptr)); in main() 210 x_ptr[0] = 0.0; in main() 211 PetscCall(VecRestoreArray(mu[0], &x_ptr)); in main() 220 PetscCall(VecGetArray(q, &x_ptr)); in main() 221 PetscCall(PetscPrintf(PETSC_COMM_WORLD, "\n cost function=%g\n", (double)(x_ptr[0] - ctx.Pm))); in main() 222 PetscCall(VecRestoreArray(q, &x_ptr)); in main() 224 PetscCall(VecGetArray(mu[0], &x_ptr)); in main() 225 PetscCall(PetscPrintf(PETSC_COMM_WORLD, "\n gradient=%g\n", (double)x_ptr[0])); in main() 226 PetscCall(VecRestoreArray(mu[0], &x_ptr)); in main()
|
| H A D | ex9adj.c | 192 PetscScalar *x_ptr, *y_ptr; in main() local 305 PetscCall(VecGetArray(mu[0], &x_ptr)); in main() 306 x_ptr[0] = -1.0; in main() 307 PetscCall(VecRestoreArray(mu[0], &x_ptr)); in main() 347 PetscCall(VecGetArray(mu[0], &x_ptr)); in main() 348 x_ptr[0] = -1.0; in main() 349 PetscCall(VecRestoreArray(mu[0], &x_ptr)); in main() 358 PetscCall(VecGetArray(q, &x_ptr)); in main() 359 PetscCall(PetscPrintf(PETSC_COMM_WORLD, "\n cost function=%g\n", (double)(x_ptr[0] - ctx.Pm))); in main() 360 PetscCall(VecRestoreArray(q, &x_ptr)); in main()
|
| /petsc/src/ts/tutorials/eimex/ |
| H A D | ct_vdp_imex.c | 47 PetscScalar *x_ptr; in main() local 78 PetscCall(VecGetArray(ref, &x_ptr)); in main() 96 x_ptr[0] = 1.619084329683235; in main() 97 x_ptr[1] = -0.803530465176385; in main() 116 PetscCall(VecGetArray(x, &x_ptr)); in main() 117 x_ptr[0] = 2.; in main() 118 …x_ptr[1] = -2. / 3. + 10. / 81. * (user.mu) - 292. / 2187. * (user.mu) * (user.mu) - 1814. / 19683… in main() 143 PetscCall(VecGetArray(x, &x_ptr)); in main() 144 …wcol[1], (double)dt, (double)norm, (double)PetscRealPart(x_ptr[0]), (double)PetscRealPart(x_ptr[1]… in main() 145 PetscCall(VecRestoreArray(x, &x_ptr)); in main()
|
| /petsc/src/ksp/pc/impls/tfs/ |
| H A D | xxt.c | 174 PetscScalar *x_ptr; in xxt_generate() local 287 x_ptr = x; in xxt_generate() 293 PetscCallBLAS("BLASdot", uu[k] = BLASdot_(&dlen, u + off, &i1, x_ptr, &i1)); in xxt_generate() 294 x_ptr += len; in xxt_generate() 302 x_ptr = x; in xxt_generate() 308 PetscCallBLAS("BLASaxpy", BLASaxpy_(&dlen, &uu[k], x_ptr, &i1, z + off, &i1)); in xxt_generate() 309 x_ptr += len; in xxt_generate() 356 x_ptr = (PetscScalar *)malloc(xxt_max_nnz * sizeof(PetscScalar)); in xxt_generate() 357 PetscCall(PCTFS_rvec_copy(x_ptr, x, xxt_nnz)); in xxt_generate() 359 x = x_ptr; in xxt_generate() [all …]
|
| H A D | xyt.c | 227 PetscScalar *x_ptr, *y_ptr; in xyt_generate() local 371 x_ptr = x; in xyt_generate() 377 PetscCallBLAS("BLASaxpy", BLASaxpy_(&dlen, &uu[k], x_ptr, &i1, z + off, &i1)); in xyt_generate() 378 x_ptr += len; in xyt_generate() 426 x_ptr = (PetscScalar *)malloc(xt_max_nnz * sizeof(PetscScalar)); in xyt_generate() 427 PetscCall(PCTFS_rvec_copy(x_ptr, x, xt_nnz)); in xyt_generate() 429 x = x_ptr; in xyt_generate() 430 x_ptr += xt_nnz; in xyt_generate() 433 PetscCall(PCTFS_rvec_copy(x_ptr, v + off, len)); in xyt_generate() 437 xcol_vals[i] = x_ptr; in xyt_generate() [all …]
|
| /petsc/src/ts/tutorials/power_grid/stability_9bus/ |
| H A D | ex9busoptfd.c | 825 PetscScalar *x_ptr; in main() local 907 PetscCall(VecGetArray(p, &x_ptr)); in main() 908 x_ptr[0] = PG[0]; in main() 909 x_ptr[1] = PG[1]; in main() 910 x_ptr[2] = PG[2]; in main() 911 PetscCall(VecRestoreArray(p, &x_ptr)); in main() 921 PetscCall(VecGetArray(lowerb, &x_ptr)); in main() 922 x_ptr[0] = 0.5; in main() 923 x_ptr[1] = 0.5; in main() 924 x_ptr[2] = 0.5; in main() [all …]
|
| H A D | ex9busopt.c | 1032 PetscScalar *x_ptr; in main() local 1147 PetscCall(VecGetArray(p, &x_ptr)); in main() 1148 x_ptr[0] = PG[0]; in main() 1149 x_ptr[1] = PG[1]; in main() 1150 x_ptr[2] = PG[2]; in main() 1151 PetscCall(VecRestoreArray(p, &x_ptr)); in main() 1160 PetscCall(VecGetArray(lowerb, &x_ptr)); in main() 1161 x_ptr[0] = 0.5; in main() 1162 x_ptr[1] = 0.5; in main() 1163 x_ptr[2] = 0.5; in main() [all …]
|
| /petsc/src/mat/impls/sbaij/seq/ |
| H A D | sbaij2.c | 677 const PetscScalar *x, *x_ptr, *xb; in MatMult_SeqSBAIJ_N() local 689 x_ptr = x; in MatMult_SeqSBAIJ_N() 708 xb = x_ptr + bs * (*idx++); in MatMult_SeqSBAIJ_N() 1192 const PetscScalar *x, *x_ptr, *xb; in MatMultAdd_SeqSBAIJ_N() local 1202 x_ptr = x; in MatMultAdd_SeqSBAIJ_N() 1222 xb = x_ptr + bs * (*idx++); in MatMultAdd_SeqSBAIJ_N()
|