| /petsc/src/mat/impls/blockmat/seq/ |
| H A D | blockmat.c | 416 Mat_BlockMat *bmat = (Mat_BlockMat *)mat->data; in MatDestroy_BlockMat() local
420 PetscCall(VecDestroy(&bmat->right)); in MatDestroy_BlockMat()
421 PetscCall(VecDestroy(&bmat->left)); in MatDestroy_BlockMat()
422 PetscCall(VecDestroy(&bmat->middle)); in MatDestroy_BlockMat()
423 PetscCall(VecDestroy(&bmat->workb)); in MatDestroy_BlockMat()
424 if (bmat->diags) { in MatDestroy_BlockMat()
425 for (i = 0; i < mat->rmap->n / mat->rmap->bs; i++) PetscCall(MatDestroy(&bmat->diags[i])); in MatDestroy_BlockMat()
427 if (bmat->a) { in MatDestroy_BlockMat()
428 for (i = 0; i < bmat->nz; i++) PetscCall(MatDestroy(&bmat->a[i])); in MatDestroy_BlockMat()
430 PetscCall(MatSeqXAIJFreeAIJ(mat, (PetscScalar **)&bmat->a, &bmat->j, &bmat->i)); in MatDestroy_BlockMat()
[all …]
|
| /petsc/src/ts/impls/implicit/glle/ |
| H A D | glle.c | 133 …PetscScalar *ImV, *H, *bmat, *workscalar, *c = scheme->c, *a = scheme->a, *b = scheme->b, *u = sch… in TSGLLESchemeCreate() local
137 …PetscCall(PetscMalloc7(PetscSqr(r), &ImV, 3 * s, &H, 3 * ss, &bmat, lwork, &workscalar, 5 * (3 + r… in TSGLLESchemeCreate()
225 bmat[0 + 0 * ss] = 1.; in TSGLLESchemeCreate()
226 bmat[0 + 1 * ss] = 0.; in TSGLLESchemeCreate()
227 bmat[0 + 2 * ss] = 0.; in TSGLLESchemeCreate()
228 bmat[1 + 0 * ss] = 1.; in TSGLLESchemeCreate()
229 bmat[1 + 1 * ss] = 1.; in TSGLLESchemeCreate()
230 bmat[1 + 2 * ss] = 0.; in TSGLLESchemeCreate()
231 bmat[2 + 0 * ss] = 0.; in TSGLLESchemeCreate()
232 bmat[2 + 1 * ss] = 0.; in TSGLLESchemeCreate()
[all …]
|
| /petsc/src/mat/impls/baij/mpi/ |
| H A D | mpibaij.c | 692 Mat_SeqBAIJ *amat = (Mat_SeqBAIJ *)baij->A->data, *bmat = (Mat_SeqBAIJ *)baij->B->data; in MatNorm_MPIBAIJ() local
708 v = bmat->a; in MatNorm_MPIBAIJ()
709 nz = bmat->nz * bs2; in MatNorm_MPIBAIJ()
738 v = bmat->a; in MatNorm_MPIBAIJ()
739 jj = bmat->j; in MatNorm_MPIBAIJ()
740 for (i = 0; i < bmat->nz; i++) { in MatNorm_MPIBAIJ()
748 PetscCall(VecSetValuesBlocked(col, bmat->nbs, garray, array, ADD_VALUES)); in MatNorm_MPIBAIJ()
771 v = bmat->a + bs2 * bmat->i[j]; in MatNorm_MPIBAIJ()
772 nz = bmat->i[j + 1] - bmat->i[j]; in MatNorm_MPIBAIJ()
|
| /petsc/src/mat/impls/aij/mpi/ |
| H A D | mpiov.c | 287 Mat amat, bmat; in MatIncreaseOverlap_MPIAIJ_Send_Scalable() local
295 PetscCall(MatMPIAIJGetSeqAIJ(mat, &amat, &bmat, &gcols)); in MatIncreaseOverlap_MPIAIJ_Send_Scalable()
299 PetscCall(MatGetRowIJ(bmat, 0, PETSC_FALSE, PETSC_FALSE, &bn, &bi, &bj, &done)); in MatIncreaseOverlap_MPIAIJ_Send_Scalable()
398 PetscCall(MatRestoreRowIJ(bmat, 0, PETSC_FALSE, PETSC_FALSE, &bn, &bi, &bj, &done)); in MatIncreaseOverlap_MPIAIJ_Send_Scalable()
408 Mat amat, bmat; in MatIncreaseOverlap_MPIAIJ_Local_Scalable() local
416 PetscCall(MatMPIAIJGetSeqAIJ(mat, &amat, &bmat, &gcols)); in MatIncreaseOverlap_MPIAIJ_Local_Scalable()
419 PetscCall(MatGetRowIJ(bmat, 0, PETSC_FALSE, PETSC_FALSE, &bn, &bi, &bj, &done)); in MatIncreaseOverlap_MPIAIJ_Local_Scalable()
465 PetscCall(MatRestoreRowIJ(bmat, 0, PETSC_FALSE, PETSC_FALSE, &bn, &bi, &bj, &done)); in MatIncreaseOverlap_MPIAIJ_Local_Scalable()
|
| H A D | mpiaij.c | 1809 Mat_SeqAIJ *amat = (Mat_SeqAIJ *)aij->A->data, *bmat = (Mat_SeqAIJ *)aij->B->data; in MatNorm_MPIAIJ() local
1827 for (i = 0; i < bmat->nz; i++) { in MatNorm_MPIAIJ()
1833 PetscCall(PetscLogFlops(2.0 * amat->nz + 2.0 * bmat->nz)); in MatNorm_MPIAIJ()
1850 jj = bmat->j; in MatNorm_MPIAIJ()
1851 for (j = 0; j < bmat->nz; j++) array[*jj++] += PetscAbsScalar(*v++); in MatNorm_MPIAIJ()
1868 v = PetscSafePointerPlusOffset(bmata, bmat->i[j]); in MatNorm_MPIAIJ()
1869 for (i = 0; i < bmat->i[j + 1] - bmat->i[j]; i++) { in MatNorm_MPIAIJ()
1876 PetscCall(PetscLogFlops(PetscMax(amat->nz + bmat->nz - 1, 0))); in MatNorm_MPIAIJ()
|
| /petsc/src/mat/impls/sbaij/mpi/ |
| H A D | mpisbaij.c | 709 Mat_SeqBAIJ *bmat = (Mat_SeqBAIJ *)baij->B->data; in MatNorm_MPISBAIJ() local
738 v = bmat->a; in MatNorm_MPISBAIJ()
739 jj = bmat->j; in MatNorm_MPISBAIJ()
742 nz = bmat->i[brow + 1] - bmat->i[brow]; in MatNorm_MPISBAIJ()
|