double LatticeTransportMaterial :: giveCharacteristicValue(MatResponseMode mode, GaussPoint *gp, TimeStep *tStep) { LatticeTransportMaterialStatus *status = static_cast< LatticeTransportMaterialStatus * >( this->giveStatus(gp) ); double suction = status->giveTempField().at(1); if ( mode == Capacity ) { return computeCapacity(suction, gp); } else if ( mode == Conductivity ) { return computeConductivity(suction, gp, tStep); } else { OOFEM_ERROR("unknown mode"); } return 0; // to make compiler happy }
Real DiffTensorKernel::computeQpJacobian() { RealTensorValue k = computeConductivity(_t, _qp); return k * _grad_test[_i][_qp] * _grad_phi[_j][_qp]; }
Real DiffTensorKernel::computeQpResidual() { RealTensorValue k = computeConductivity(_t, _qp); return k * _grad_test[_i][_qp] * _grad_u[_qp] - 4.0; }