Hong <[email protected]> writes: > Jed : > >> >> Is it better this way or as a fallback when !A->ops->rart? MatPtAP >> handles other combinations like MAIJ. >> > > Do you mean > if ( !A->ops->rart) { > Mat Rt; > ierr = MatTranspose(R,MAT_INITIAL_MATRIX,&Rt);CHKERRQ(ierr); > ierr = MatMatMatMult(R,A,Rt,scall,fill,C);CHKERRQ(ierr); > ierr = MatDestroy(&Rt);CHKERRQ(ierr); > } > This does NOT work for most matrix formats because we do not have fallbacks > for MatTranspose() and MatMatMult().
That's fine; they'll trigger an error and we'll be able to see from the stack that it can be made to work by either implementing the appropriate MatRARt or MatTranspose and MatMatMatMult.
signature.asc
Description: PGP signature
