68 const Teuchos::SerialDenseMatrix<ordinal_type, value_type>& Input,
69 Teuchos::SerialDenseMatrix<ordinal_type, value_type>& Result,
70 ordinal_type m)
const {
71 ordinal_type n=Input.numRows();
72 Teuchos::SerialDenseMatrix<ordinal_type, value_type> G(Teuchos::Copy,
A);
73 Teuchos::SerialDenseMatrix<ordinal_type, value_type> z(n,1);
74 for (ordinal_type
j=0;
j<m;
j++){
76 for (ordinal_type i=0; i<n; i++)
77 z(i,0)=Input(i,0)/
A(i,i);
81 for (ordinal_type i=0; i<n; i++){
82 for (ordinal_type
j=0;
j<n;
j++){
86 G(i,
j)=-
A(i,
j)/
A(i,i);
92 Result.multiply(Teuchos::NO_TRANS,Teuchos::NO_TRANS,1.0, G, z, 1.0);