PreviousSolutionFunction<Scalar>::PreviousSolutionFunction(TSolutionPtr<Scalar> soln, LinearTermPtr solnExpression, bool multiplyFluxesByCellParity) : TFunction<Scalar>(solnExpression->rank()) { _soln = soln; _solnExpression = solnExpression; _overrideMeshCheck = false; if ((solnExpression->termType() == FLUX) && multiplyFluxesByCellParity) { TFunctionPtr<double> parity = TFunction<double>::sideParity(); _solnExpression = parity * solnExpression; } }
PreviousSolutionFunction::PreviousSolutionFunction(SolutionPtr soln, LinearTermPtr solnExpression, bool multiplyFluxesByCellParity) : Function(solnExpression->rank()) { _soln = soln; _solnExpression = solnExpression; _overrideMeshCheck = false; if ((solnExpression->termType() == FLUX) && multiplyFluxesByCellParity) { FunctionPtr parity = Teuchos::rcp( new SideParityFunction ); _solnExpression = parity * solnExpression; } }