Lines Matching refs:uik
566 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_N() local
581 PetscCall(PetscMalloc3(bs2, &dk, bs2, &uik, bs, &work)); in MatCholeskyFactorNumeric_SeqSBAIJ_N()
652 PetscCall(PetscArrayzero(uik, bs2)); in MatCholeskyFactorNumeric_SeqSBAIJ_N()
653 PetscKernel_A_gets_A_minus_B_times_C(bs, uik, diag, u); in MatCholeskyFactorNumeric_SeqSBAIJ_N()
656 PetscKernel_A_gets_A_plus_Btranspose_times_C(bs, dk, uik, u); in MatCholeskyFactorNumeric_SeqSBAIJ_N()
660 PetscCall(PetscArraycpy(ba + ili * bs2, uik, bs2)); in MatCholeskyFactorNumeric_SeqSBAIJ_N()
670 PetscKernel_A_gets_A_plus_Btranspose_times_C(bs, rtmp_ptr, uik, u); in MatCholeskyFactorNumeric_SeqSBAIJ_N()
715 PetscCall(PetscFree3(dk, uik, work)); in MatCholeskyFactorNumeric_SeqSBAIJ_N()
739 MatScalar *ba = b->a, *aa, *ap, *dk, *uik; in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering() local
751 PetscCall(PetscMalloc3(bs2, &dk, bs2, &uik, bs, &work)); in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering()
784 PetscCall(PetscArrayzero(uik, bs2)); in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering()
785 PetscKernel_A_gets_A_minus_B_times_C(bs, uik, diag, u); in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering()
788 PetscKernel_A_gets_A_plus_Btranspose_times_C(bs, dk, uik, u); in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering()
792 PetscCall(PetscArraycpy(ba + ili * bs2, uik, bs2)); in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering()
802 PetscKernel_A_gets_A_plus_Btranspose_times_C(bs, rtmp_ptr, uik, u); in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering()
847 PetscCall(PetscFree3(dk, uik, work)); in MatCholeskyFactorNumeric_SeqSBAIJ_N_NaturalOrdering()
873 MatScalar *u, *diag, *rtmp, *rtmp_ptr, dk[4], uik[4]; in MatCholeskyFactorNumeric_SeqSBAIJ_2() local
960 uik[0] = -(diag[0] * u[0] + diag[2] * u[1]); in MatCholeskyFactorNumeric_SeqSBAIJ_2()
961 uik[1] = -(diag[1] * u[0] + diag[3] * u[1]); in MatCholeskyFactorNumeric_SeqSBAIJ_2()
962 uik[2] = -(diag[0] * u[2] + diag[2] * u[3]); in MatCholeskyFactorNumeric_SeqSBAIJ_2()
963 uik[3] = -(diag[1] * u[2] + diag[3] * u[3]); in MatCholeskyFactorNumeric_SeqSBAIJ_2()
966 dk[0] += uik[0] * u[0] + uik[1] * u[1]; in MatCholeskyFactorNumeric_SeqSBAIJ_2()
967 dk[1] += uik[2] * u[0] + uik[3] * u[1]; in MatCholeskyFactorNumeric_SeqSBAIJ_2()
968 dk[2] += uik[0] * u[2] + uik[1] * u[3]; in MatCholeskyFactorNumeric_SeqSBAIJ_2()
969 dk[3] += uik[2] * u[2] + uik[3] * u[3]; in MatCholeskyFactorNumeric_SeqSBAIJ_2()
974 PetscCall(PetscArraycpy(ba + ili * 4, uik, 4)); in MatCholeskyFactorNumeric_SeqSBAIJ_2()
984 rtmp_ptr[0] += uik[0] * u[0] + uik[1] * u[1]; in MatCholeskyFactorNumeric_SeqSBAIJ_2()
985 rtmp_ptr[1] += uik[2] * u[0] + uik[3] * u[1]; in MatCholeskyFactorNumeric_SeqSBAIJ_2()
986 rtmp_ptr[2] += uik[0] * u[2] + uik[1] * u[3]; in MatCholeskyFactorNumeric_SeqSBAIJ_2()
987 rtmp_ptr[3] += uik[2] * u[2] + uik[3] * u[3]; in MatCholeskyFactorNumeric_SeqSBAIJ_2()
1051 MatScalar *ba = b->a, *aa, *ap, dk[8], uik[8]; in MatCholeskyFactorNumeric_SeqSBAIJ_2_NaturalOrdering() local
1102 uik[0] = -(diag[0] * u[0] + diag[2] * u[1]); in MatCholeskyFactorNumeric_SeqSBAIJ_2_NaturalOrdering()
1103 uik[1] = -(diag[1] * u[0] + diag[3] * u[1]); in MatCholeskyFactorNumeric_SeqSBAIJ_2_NaturalOrdering()
1104 uik[2] = -(diag[0] * u[2] + diag[2] * u[3]); in MatCholeskyFactorNumeric_SeqSBAIJ_2_NaturalOrdering()
1105 uik[3] = -(diag[1] * u[2] + diag[3] * u[3]); in MatCholeskyFactorNumeric_SeqSBAIJ_2_NaturalOrdering()
1108 dk[0] += uik[0] * u[0] + uik[1] * u[1]; in MatCholeskyFactorNumeric_SeqSBAIJ_2_NaturalOrdering()
1109 dk[1] += uik[2] * u[0] + uik[3] * u[1]; in MatCholeskyFactorNumeric_SeqSBAIJ_2_NaturalOrdering()
1110 dk[2] += uik[0] * u[2] + uik[1] * u[3]; in MatCholeskyFactorNumeric_SeqSBAIJ_2_NaturalOrdering()
1111 dk[3] += uik[2] * u[2] + uik[3] * u[3]; in MatCholeskyFactorNumeric_SeqSBAIJ_2_NaturalOrdering()
1116 PetscCall(PetscArraycpy(ba + ili * 4, uik, 4)); in MatCholeskyFactorNumeric_SeqSBAIJ_2_NaturalOrdering()
1126 rtmp_ptr[0] += uik[0] * u[0] + uik[1] * u[1]; in MatCholeskyFactorNumeric_SeqSBAIJ_2_NaturalOrdering()
1127 rtmp_ptr[1] += uik[2] * u[0] + uik[3] * u[1]; in MatCholeskyFactorNumeric_SeqSBAIJ_2_NaturalOrdering()
1128 rtmp_ptr[2] += uik[0] * u[2] + uik[1] * u[3]; in MatCholeskyFactorNumeric_SeqSBAIJ_2_NaturalOrdering()
1129 rtmp_ptr[3] += uik[2] * u[2] + uik[3] * u[3]; in MatCholeskyFactorNumeric_SeqSBAIJ_2_NaturalOrdering()