tmp<scalarField> faMatrix<scalar>::residual() const { scalarField boundaryDiag(psi_.size(), 0.0); addBoundaryDiag(boundaryDiag, 0); // Make a copy of interfaces: no longer a reference // HJ, 20/Nov/2007 lduInterfaceFieldPtrsList interfaces = psi_.boundaryField().interfaces(); tmp<scalarField> tres ( lduMatrix::residual ( psi_.internalField(), source_ - boundaryDiag*psi_.internalField(), boundaryCoeffs_, interfaces, 0 ) ); addBoundarySource(tres()); return tres; }
Foam::tmp<Foam::scalarField> Foam::fvMatrix<Foam::scalar>::residual() const { scalarField boundaryDiag(psi_.size(), 0.0); addBoundaryDiag(boundaryDiag, 0); tmp<scalarField> tres ( lduMatrix::residual ( psi_.primitiveField(), source_ - boundaryDiag*psi_.primitiveField(), boundaryCoeffs_, psi_.boundaryField().scalarInterfaces(), 0 ) ); addBoundarySource(tres.ref()); return tres; }