Home
last modified time | relevance | path

Searched refs:acols (Results 1 – 6 of 6) sorted by relevance

/petsc/src/mat/impls/aij/seq/crl/
H A Dcrl.c19 if (aijcrl) PetscCall(PetscFree2(aijcrl->acols, aijcrl->icols)); in MatDestroy_SeqAIJCRL()
39 PetscScalar *acols; in MatSeqAIJCRL_create_aijcrl() local
46 PetscCall(PetscFree2(aijcrl->acols, aijcrl->icols)); in MatSeqAIJCRL_create_aijcrl()
47 PetscCall(PetscMalloc2(rmax * m, &aijcrl->acols, rmax * m, &aijcrl->icols)); in MatSeqAIJCRL_create_aijcrl()
48 acols = aijcrl->acols; in MatSeqAIJCRL_create_aijcrl()
52 acols[j * m + i] = *aa++; in MatSeqAIJCRL_create_aijcrl()
56 acols[j * m + i] = 0.0; in MatSeqAIJCRL_create_aijcrl()
91 PetscScalar *acols = aijcrl->acols; in MatMult_AIJCRL() local
115 fortranmultcrl_(&m, &rmax, x, y, icols, acols); in MatMult_AIJCRL()
119 for (j = 0; j < m; j++) y[j] = acols[j] * x[icols[j]]; in MatMult_AIJCRL()
[all …]
H A Dcrl.h11 PetscScalar *acols; /* values of nonzeros, stored as icols */ member
/petsc/src/mat/impls/aij/mpi/crl/
H A Dmcrl.c23 PetscCall(PetscFree2(aijcrl->acols, aijcrl->icols)); in MatDestroy_MPIAIJCRL()
44 PetscScalar *aa = Aij->a, *ba = Bij->a, *acols, *array; in MatMPIAIJCRL_create_aijcrl() local
53 PetscCall(PetscFree2(aijcrl->acols, aijcrl->icols)); in MatMPIAIJCRL_create_aijcrl()
54 PetscCall(PetscMalloc2(rmax * m, &aijcrl->acols, rmax * m, &aijcrl->icols)); in MatMPIAIJCRL_create_aijcrl()
55 acols = aijcrl->acols; in MatMPIAIJCRL_create_aijcrl()
59 acols[j * m + i] = *aa++; in MatMPIAIJCRL_create_aijcrl()
63 acols[j * m + i] = *ba++; in MatMPIAIJCRL_create_aijcrl()
67 acols[j * m + i] = 0.0; in MatMPIAIJCRL_create_aijcrl()
/petsc/src/mat/impls/aij/seq/crl/ftn-kernels/
H A Dfmultcrl.F9011 PetscScalar, intent(in) :: x(0:m - 1), acols(m, rmax)
16 y(1:m) = acols(1:m, 1)*x(icols(1:m, 1))
18 y(1:m) = y(1:m) + acols(1:m, i)*x(icols(1:m, i))
/petsc/src/ksp/ksp/tutorials/
H A Dex67.c197 MatStencil *acols; in FormMatrix() local
203 PetscCall(PetscMalloc1(M, &acols)); in FormMatrix()
206 acols[i].i = i; in FormMatrix()
209 PetscCall(MatSetValuesStencil(jac, 1, &row, M, acols, avals, ADD_VALUES)); in FormMatrix()
210 PetscCall(PetscFree(acols)); in FormMatrix()
/petsc/src/mat/impls/aij/mpi/
H A Dmpiaij.c1575 PetscInt *acols = dnnz, *bcols = onnz; /* Repurpose now-unneeded arrays */ in MatPermute_MPIAIJ() local
1579 for (; j < PetscMin(rowlen, j0 + m); j++) acols[j - j0] = cdest[aj[ai[i] + j]]; in MatPermute_MPIAIJ()
1580 PetscCall(MatSetValues(Aperm, 1, &rdest[i], j - j0, acols, aa + ai[i] + j0, INSERT_VALUES)); in MatPermute_MPIAIJ()