Lines Matching refs:Brhs
719 PetscErrorCode TSGetRHSMats_Private(TS ts, Mat *Arhs, Mat *Brhs) in TSGetRHSMats_Private() argument
726 if (Brhs) *Brhs = NULL; in TSGetRHSMats_Private()
741 if (flg && !ijacobian && ts->Arhs == ts->Brhs) { in TSGetRHSMats_Private()
749 if (Brhs) { in TSGetRHSMats_Private()
750 if (!ts->Brhs) { in TSGetRHSMats_Private()
753 PetscCall(MatDuplicate(B, MAT_DO_NOT_COPY_VALUES, &ts->Brhs)); in TSGetRHSMats_Private()
755 ts->Brhs = B; in TSGetRHSMats_Private()
760 ts->Brhs = ts->Arhs; in TSGetRHSMats_Private()
763 *Brhs = ts->Brhs; in TSGetRHSMats_Private()
849 PetscCheck(B == ts->Brhs, PetscObjectComm((PetscObject)ts), PETSC_ERR_SUP, "Invalid Bmat"); in TSRecoverRHSJacobian()
853 if (B && B == ts->Brhs && A != B) { in TSRecoverRHSJacobian()
945 Mat Arhs = NULL, Brhs = NULL; in TSComputeIJacobian() local
948 if (rhsjacobian) PetscCall(TSGetRHSMats_Private(ts, &Arhs, &Brhs)); in TSComputeIJacobian()
993 PetscCall(TSComputeRHSJacobian(ts, t, U, Arhs, Brhs)); in TSComputeIJacobian()
995 if (A != B) PetscCall(MatAXPY(B, -1, Brhs, ts->axpy_pattern)); in TSComputeIJacobian()
1168 PetscCall(MatDestroy(&ts->Brhs)); in TSSetRHSJacobian()
1169 ts->Brhs = Pmat; in TSSetRHSJacobian()
2533 if (Pmat && Pmat == ts->Brhs) { in TSSetUp()
2534 PetscCall(MatDuplicate(ts->Brhs, MAT_COPY_VALUES, &Pmat)); in TSSetUp()
2599 PetscCall(MatDestroy(&ts->Brhs)); in TSReset()
4647 Mat Arhs, Brhs; in TSComputeRHSFunctionLinear() local
4650 PetscCall(TSGetRHSMats_Private(ts, &Arhs, &Brhs)); in TSComputeRHSFunctionLinear()
4652 PetscCall(TSRecoverRHSJacobian(ts, Arhs, Brhs)); in TSComputeRHSFunctionLinear()
4653 PetscCall(TSComputeRHSJacobian(ts, t, U, Arhs, Brhs)); in TSComputeRHSFunctionLinear()