| /petsc/src/mat/utils/ |
| H A D | freespace.c | 139 PetscInt row, nnz, *uj, *array, total; in PetscFreeSpaceContiguous_Cholesky() local 152 uj = space + ui[row]; in PetscFreeSpaceContiguous_Cholesky() 153 PetscCall(PetscArraycpy(uj, array + 1, nnz)); in PetscFreeSpaceContiguous_Cholesky() 154 uj[nnz] = array[0]; /* diagonal */ in PetscFreeSpaceContiguous_Cholesky()
|
| /petsc/src/mat/impls/aij/seq/bas/ |
| H A D | basfactor.c | 13 PetscInt ncols, *cols, *uj; in MatICCFactorSymbolic_SeqAIJ_Bas() local 31 PetscCall(PetscMalloc1(ui[am] + 1, &uj)); in MatICCFactorSymbolic_SeqAIJ_Bas() 32 cols = uj; in MatICCFactorSymbolic_SeqAIJ_Bas() 56 PetscCall(spbas_matrix_to_crs(Pattern_P, NULL, &ui, &uj)); in MatICCFactorSymbolic_SeqAIJ_Bas() 65 b->j = uj; in MatICCFactorSymbolic_SeqAIJ_Bas()
|
| /petsc/src/ts/tutorials/multirate/ |
| H A D | ex7.c | 549 PetscScalar *u, *uj, xj, xi; in FVSample() local 558 PetscCall(PetscMalloc1(dof, &uj)); in FVSample() 570 …physics.sample)(ctx->physics.user, ctx->initial, ctx->bctype, ctx->xmin, ctx->xmax, time, xj, uj)); in FVSample() 571 for (k = 0; k < dof; k++) u[i * dof + k] += ((j == 0 || j == N) ? 0.5 : 1.0) * uj[k] / N; in FVSample() 579 …physics.sample)(ctx->physics.user, ctx->initial, ctx->bctype, ctx->xmin, ctx->xmax, time, xj, uj)); in FVSample() 580 for (k = 0; k < dof; k++) u[i * dof + k] += ((j == 0 || j == N) ? 0.5 : 1.0) * uj[k] / N; in FVSample() 588 …physics.sample)(ctx->physics.user, ctx->initial, ctx->bctype, ctx->xmin, ctx->xmax, time, xj, uj)); in FVSample() 589 for (k = 0; k < dof; k++) u[i * dof + k] += ((j == 0 || j == N) ? 0.5 : 1.0) * uj[k] / N; in FVSample() 594 PetscCall(PetscFree(uj)); in FVSample()
|
| H A D | ex8.c | 125 PetscScalar *u, *uj, xj, xi; in FVSample_3WaySplit() local 135 PetscCall(PetscMalloc1(dof, &uj)); in FVSample_3WaySplit() 147 …sics2.sample2)(ctx->physics2.user, ctx->initial, ctx->bctype, ctx->xmin, ctx->xmax, time, xj, uj)); in FVSample_3WaySplit() 148 for (k = 0; k < dof; k++) u[i * dof + k] += ((j == 0 || j == N) ? 0.5 : 1.0) * uj[k] / N; in FVSample_3WaySplit() 156 …sics2.sample2)(ctx->physics2.user, ctx->initial, ctx->bctype, ctx->xmin, ctx->xmax, time, xj, uj)); in FVSample_3WaySplit() 157 for (k = 0; k < dof; k++) u[i * dof + k] += ((j == 0 || j == N) ? 0.5 : 1.0) * uj[k] / N; in FVSample_3WaySplit() 165 …sics2.sample2)(ctx->physics2.user, ctx->initial, ctx->bctype, ctx->xmin, ctx->xmax, time, xj, uj)); in FVSample_3WaySplit() 166 for (k = 0; k < dof; k++) u[i * dof + k] += ((j == 0 || j == N) ? 0.5 : 1.0) * uj[k] / N; in FVSample_3WaySplit() 174 …sics2.sample2)(ctx->physics2.user, ctx->initial, ctx->bctype, ctx->xmin, ctx->xmax, time, xj, uj)); in FVSample_3WaySplit() 175 for (k = 0; k < dof; k++) u[i * dof + k] += ((j == 0 || j == N) ? 0.5 : 1.0) * uj[k] / N; in FVSample_3WaySplit() [all …]
|
| H A D | ex6.c | 132 PetscScalar *u, *uj, xj, xi; in FVSample_2WaySplit() local 142 PetscCall(PetscMalloc1(dof, &uj)); in FVSample_2WaySplit() 152 …sics2.sample2)(ctx->physics2.user, ctx->initial, ctx->bctype, ctx->xmin, ctx->xmax, time, xj, uj)); in FVSample_2WaySplit() 153 for (k = 0; k < dof; k++) u[i * dof + k] += ((j == 0 || j == N) ? 0.5 : 1.0) * uj[k] / N; in FVSample_2WaySplit() 161 …sics2.sample2)(ctx->physics2.user, ctx->initial, ctx->bctype, ctx->xmin, ctx->xmax, time, xj, uj)); in FVSample_2WaySplit() 162 for (k = 0; k < dof; k++) u[i * dof + k] += ((j == 0 || j == N) ? 0.5 : 1.0) * uj[k] / N; in FVSample_2WaySplit() 170 …sics2.sample2)(ctx->physics2.user, ctx->initial, ctx->bctype, ctx->xmin, ctx->xmax, time, xj, uj)); in FVSample_2WaySplit() 171 for (k = 0; k < dof; k++) u[i * dof + k] += ((j == 0 || j == N) ? 0.5 : 1.0) * uj[k] / N; in FVSample_2WaySplit() 176 PetscCall(PetscFree(uj)); in FVSample_2WaySplit()
|
| H A D | ex4.c | 476 PetscScalar *u, *uj, xj, xi; in FVSample_2WaySplit() local 486 PetscCall(PetscMalloc1(dof, &uj)); in FVSample_2WaySplit() 496 …sics2.sample2)(ctx->physics2.user, ctx->initial, ctx->bctype, ctx->xmin, ctx->xmax, time, xj, uj)); in FVSample_2WaySplit() 497 for (k = 0; k < dof; k++) u[i * dof + k] += ((j == 0 || j == N) ? 0.5 : 1.0) * uj[k] / N; in FVSample_2WaySplit() 505 …sics2.sample2)(ctx->physics2.user, ctx->initial, ctx->bctype, ctx->xmin, ctx->xmax, time, xj, uj)); in FVSample_2WaySplit() 506 for (k = 0; k < dof; k++) u[i * dof + k] += ((j == 0 || j == N) ? 0.5 : 1.0) * uj[k] / N; in FVSample_2WaySplit() 514 …sics2.sample2)(ctx->physics2.user, ctx->initial, ctx->bctype, ctx->xmin, ctx->xmax, time, xj, uj)); in FVSample_2WaySplit() 515 for (k = 0; k < dof; k++) u[i * dof + k] += ((j == 0 || j == N) ? 0.5 : 1.0) * uj[k] / N; in FVSample_2WaySplit() 520 PetscCall(PetscFree(uj)); in FVSample_2WaySplit()
|
| H A D | finitevolume1d.c | 622 PetscScalar *u, *uj; in FVSample() local 630 PetscCall(PetscMalloc1(dof, &uj)); in FVSample() 638 …physics.sample)(ctx->physics.user, ctx->initial, ctx->bctype, ctx->xmin, ctx->xmax, time, xj, uj)); in FVSample() 639 for (k = 0; k < dof; k++) u[i * dof + k] += ((j == 0 || j == N) ? 0.5 : 1.0) * uj[k] / N; in FVSample() 643 PetscCall(PetscFree(uj)); in FVSample()
|
| /petsc/src/mat/impls/baij/seq/ |
| H A D | baijfact.c | 1025 …PetscInt nlnk, *lnk, *lnk_lvl = NULL, ncols, ncols_upper, *cols, *cols_lvl, *uj, **uj_pt… in MatICCFactorSymbolic_SeqBAIJ() local 1060 uj = b->j; in MatICCFactorSymbolic_SeqBAIJ() 1063 for (j = 0; j < ui[i]; j++) *uj++ = *aj++; in MatICCFactorSymbolic_SeqBAIJ() 1185 PetscCall(PetscMalloc1(ui[am] + 1, &uj)); in MatICCFactorSymbolic_SeqBAIJ() 1186 PetscCall(PetscFreeSpaceContiguous(&free_space, uj)); in MatICCFactorSymbolic_SeqBAIJ() 1200 b->j = uj; in MatICCFactorSymbolic_SeqBAIJ() 1255 PetscInt nlnk, *lnk, ncols, ncols_upper, *cols, *uj, **ui_ptr, *uj_ptr; in MatCholeskyFactorSymbolic_SeqBAIJ() local 1369 PetscCall(PetscMalloc1(ui[mbs] + 1, &uj)); in MatCholeskyFactorSymbolic_SeqBAIJ() 1370 PetscCall(PetscFreeSpaceContiguous(&free_space, uj)); in MatCholeskyFactorSymbolic_SeqBAIJ() 1383 b->j = uj; in MatCholeskyFactorSymbolic_SeqBAIJ()
|
| /petsc/src/mat/impls/sbaij/seq/ |
| H A D | sbaijfact.c | 229 PetscInt nlnk, *lnk, ncols, *cols, *uj, **ui_ptr, *uj_ptr, *udiag; in MatCholeskyFactorSymbolic_SeqSBAIJ() local 336 PetscCall(PetscShmgetAllocateArray(ui[mbs], sizeof(PetscInt), (void **)&uj)); in MatCholeskyFactorSymbolic_SeqSBAIJ() 337 …PetscCall(PetscFreeSpaceContiguous_Cholesky(&free_space, uj, mbs, ui, udiag)); /* store matrix fac… in MatCholeskyFactorSymbolic_SeqSBAIJ() 347 b->j = uj; in MatCholeskyFactorSymbolic_SeqSBAIJ() 394 PetscInt nlnk, *lnk, ncols, *cols, *uj, **ui_ptr, *uj_ptr; in MatCholeskyFactorSymbolic_SeqSBAIJ_inplace() local 511 PetscCall(PetscShmgetAllocateArray(ui[mbs] + 1, sizeof(PetscInt), (void **)&uj)); in MatCholeskyFactorSymbolic_SeqSBAIJ_inplace() 512 PetscCall(PetscFreeSpaceContiguous(&free_space, uj)); in MatCholeskyFactorSymbolic_SeqSBAIJ_inplace() 522 b->j = uj; in MatCholeskyFactorSymbolic_SeqSBAIJ_inplace()
|
| H A D | sbaijfact2.c | 2478 PetscInt nlnk, *lnk, *lnk_lvl = NULL, ncols, *uj, **uj_ptr, **uj_lvl_ptr; in MatICCFactorSymbolic_SeqSBAIJ() local 2510 PetscCall(PetscMalloc1(ui[am], &uj)); in MatICCFactorSymbolic_SeqSBAIJ() 2511 cols = uj; in MatICCFactorSymbolic_SeqSBAIJ() 2565 uj = uj_lvl_ptr[prow] + i; /* levels of cols */ in MatICCFactorSymbolic_SeqSBAIJ() 2566 j = *(uj - 1); in MatICCFactorSymbolic_SeqSBAIJ() 2567 PetscCall(PetscICCLLAddSorted(ncols, cols, levels, uj, am, &nlnk, lnk, lnk_lvl, lnkbt, j)); in MatICCFactorSymbolic_SeqSBAIJ() 2618 PetscCall(PetscMalloc1(ui[am], &uj)); in MatICCFactorSymbolic_SeqSBAIJ() 2619 …PetscCall(PetscFreeSpaceContiguous_Cholesky(&free_space, uj, am, ui, udiag)); /* store matrix fact… in MatICCFactorSymbolic_SeqSBAIJ() 2633 b->j = uj; in MatICCFactorSymbolic_SeqSBAIJ() 2680 PetscInt nlnk, *lnk, *lnk_lvl = NULL, ncols, *cols_lvl, *uj, **uj_ptr, **uj_lvl_ptr; in MatICCFactorSymbolic_SeqSBAIJ_inplace() local [all …]
|
| /petsc/src/mat/impls/aij/seq/ |
| H A D | aijfact.c | 2060 PetscInt ncols, ncols_upper, *cols, *ajtmp, *uj, **uj_ptr, **uj_lvl_ptr; in MatICCFactorSymbolic_SeqAIJ() local 2086 PetscCall(PetscMalloc1(ui[am] + 1, &uj)); in MatICCFactorSymbolic_SeqAIJ() 2087 cols = uj; in MatICCFactorSymbolic_SeqAIJ() 2147 uj = uj_lvl_ptr[prow] + i; /* levels of cols */ in MatICCFactorSymbolic_SeqAIJ() 2148 j = *(uj - 1); in MatICCFactorSymbolic_SeqAIJ() 2149 PetscCall(PetscICCLLAddSorted(ncols, cols, levels, uj, am, &nlnk, lnk, lnk_lvl, lnkbt, j)); in MatICCFactorSymbolic_SeqAIJ() 2202 PetscCall(PetscShmgetAllocateArray(ui[am] + 1, sizeof(PetscInt), (void **)&uj)); in MatICCFactorSymbolic_SeqAIJ() 2203 …PetscCall(PetscFreeSpaceContiguous_Cholesky(&free_space, uj, am, ui, udiag)); /* store matrix fact… in MatICCFactorSymbolic_SeqAIJ() 2214 b->j = uj; in MatICCFactorSymbolic_SeqAIJ() 2261 PetscInt nlnk, *lnk, ncols, ncols_upper, *cols, *uj, **ui_ptr, *uj_ptr, *udiag; in MatCholeskyFactorSymbolic_SeqAIJ() local [all …]
|
| /petsc/share/petsc/datafiles/meshes/ |
| H A D | square_periodic_bin.msh | 14 …D�Q�;8@R�� r@E:6���v@k��,I@F�r�Zd�"@�uj����?G���…
|
| /petsc/src/ts/tutorials/ |
| H A D | ex9.c | 1311 PetscScalar *u, *uj; in FVSample() local 1319 PetscCall(PetscMalloc1(dof, &uj)); in FVSample() 1327 …physics.sample)(ctx->physics.user, ctx->initial, ctx->bctype, ctx->xmin, ctx->xmax, time, xj, uj)); in FVSample() 1328 for (k = 0; k < dof; k++) u[i * dof + k] += ((j == 0 || j == N) ? 0.5 : 1.0) * uj[k] / N; in FVSample() 1332 PetscCall(PetscFree(uj)); in FVSample()
|