Lines Matching refs:stencil
795 MatStencil stencil[6], rowstencil; in IJacobianByHand() local
827 stencil[0].k = 0; in IJacobianByHand()
828 stencil[1].k = 0; in IJacobianByHand()
829 stencil[2].k = 0; in IJacobianByHand()
830 stencil[3].k = 0; in IJacobianByHand()
831 stencil[4].k = 0; in IJacobianByHand()
832 stencil[5].k = 0; in IJacobianByHand()
838 stencil[0].j = j - 1; in IJacobianByHand()
839 stencil[1].j = j + 1; in IJacobianByHand()
840 stencil[2].j = j; in IJacobianByHand()
841 stencil[3].j = j; in IJacobianByHand()
842 stencil[4].j = j; in IJacobianByHand()
843 stencil[5].j = j; in IJacobianByHand()
857 stencil[0].i = i; in IJacobianByHand()
858 stencil[0].c = 0; in IJacobianByHand()
860 stencil[1].i = i; in IJacobianByHand()
861 stencil[1].c = 0; in IJacobianByHand()
863 stencil[2].i = i - 1; in IJacobianByHand()
864 stencil[2].c = 0; in IJacobianByHand()
866 stencil[3].i = i + 1; in IJacobianByHand()
867 stencil[3].c = 0; in IJacobianByHand()
869 stencil[4].i = i; in IJacobianByHand()
870 stencil[4].c = 0; in IJacobianByHand()
872 stencil[5].i = i; in IJacobianByHand()
873 stencil[5].c = 1; in IJacobianByHand()
878 PetscCall(MatSetValuesStencil(A, 1, &rowstencil, 6, stencil, entries, INSERT_VALUES)); in IJacobianByHand()
879 …if (appctx->aijpc) PetscCall(MatSetValuesStencil(B, 1, &rowstencil, 6, stencil, entries, INSERT_VA… in IJacobianByHand()
880 stencil[0].c = 1; in IJacobianByHand()
882 stencil[1].c = 1; in IJacobianByHand()
884 stencil[2].c = 1; in IJacobianByHand()
886 stencil[3].c = 1; in IJacobianByHand()
888 stencil[4].c = 1; in IJacobianByHand()
890 stencil[5].c = 0; in IJacobianByHand()
893 PetscCall(MatSetValuesStencil(A, 1, &rowstencil, 6, stencil, entries, INSERT_VALUES)); in IJacobianByHand()
894 …if (appctx->aijpc) PetscCall(MatSetValuesStencil(B, 1, &rowstencil, 6, stencil, entries, INSERT_VA… in IJacobianByHand()
925 MatStencil stencil[6], rowstencil; in RHSJacobianByHand() local
957 stencil[0].k = 0; in RHSJacobianByHand()
958 stencil[1].k = 0; in RHSJacobianByHand()
959 stencil[2].k = 0; in RHSJacobianByHand()
960 stencil[3].k = 0; in RHSJacobianByHand()
961 stencil[4].k = 0; in RHSJacobianByHand()
962 stencil[5].k = 0; in RHSJacobianByHand()
969 stencil[0].j = j - 1; in RHSJacobianByHand()
970 stencil[1].j = j + 1; in RHSJacobianByHand()
971 stencil[2].j = j; in RHSJacobianByHand()
972 stencil[3].j = j; in RHSJacobianByHand()
973 stencil[4].j = j; in RHSJacobianByHand()
974 stencil[5].j = j; in RHSJacobianByHand()
988 stencil[0].i = i; in RHSJacobianByHand()
989 stencil[0].c = 0; in RHSJacobianByHand()
991 stencil[1].i = i; in RHSJacobianByHand()
992 stencil[1].c = 0; in RHSJacobianByHand()
994 stencil[2].i = i - 1; in RHSJacobianByHand()
995 stencil[2].c = 0; in RHSJacobianByHand()
997 stencil[3].i = i + 1; in RHSJacobianByHand()
998 stencil[3].c = 0; in RHSJacobianByHand()
1000 stencil[4].i = i; in RHSJacobianByHand()
1001 stencil[4].c = 0; in RHSJacobianByHand()
1003 stencil[5].i = i; in RHSJacobianByHand()
1004 stencil[5].c = 1; in RHSJacobianByHand()
1009 PetscCall(MatSetValuesStencil(A, 1, &rowstencil, 6, stencil, entries, INSERT_VALUES)); in RHSJacobianByHand()
1011 stencil[0].c = 1; in RHSJacobianByHand()
1013 stencil[1].c = 1; in RHSJacobianByHand()
1015 stencil[2].c = 1; in RHSJacobianByHand()
1017 stencil[3].c = 1; in RHSJacobianByHand()
1019 stencil[4].c = 1; in RHSJacobianByHand()
1021 stencil[5].c = 0; in RHSJacobianByHand()
1025 PetscCall(MatSetValuesStencil(A, 1, &rowstencil, 6, stencil, entries, INSERT_VALUES)); in RHSJacobianByHand()