00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037 #ifndef multivariateScheme_H
00038 #define multivariateScheme_H
00039
00040 #include "multivariateSurfaceInterpolationScheme.H"
00041 #include "surfaceFields.H"
00042
00043
00044
00045 namespace Foam
00046 {
00047
00048
00049
00050
00051
00052 template<class Type, class Scheme>
00053 class multivariateScheme
00054 :
00055 public multivariateSurfaceInterpolationScheme<Type>,
00056 public Scheme::LimiterType
00057 {
00058
00059
00060 const surfaceScalarField& faceFlux_;
00061 surfaceScalarField weights_;
00062
00063
00064
00065
00066
00067 multivariateScheme(const multivariateScheme&);
00068
00069
00070 void operator=(const multivariateScheme&);
00071
00072
00073 public:
00074
00075
00076 TypeName("multivariateScheme");
00077
00078
00079
00080
00081
00082 multivariateScheme
00083 (
00084 const fvMesh& mesh,
00085 const typename multivariateSurfaceInterpolationScheme<Type>::
00086 fieldTable& fields,
00087 const surfaceScalarField& faceFlux,
00088 Istream& schemeData
00089 );
00090
00091
00092
00093
00094
00095
00096 class fieldScheme
00097 :
00098 public multivariateSurfaceInterpolationScheme<Type>::
00099 fieldScheme
00100 {
00101
00102
00103 const surfaceScalarField& weights_;
00104
00105 public:
00106
00107
00108
00109 fieldScheme
00110 (
00111 const GeometricField<Type, fvPatchField, volMesh>& field,
00112 const surfaceScalarField& weights
00113 )
00114 :
00115 multivariateSurfaceInterpolationScheme<Type>::
00116 fieldScheme(field),
00117 weights_(weights)
00118 {}
00119
00120
00121
00122
00123
00124 tmp<surfaceScalarField> weights
00125 (
00126 const GeometricField<Type, fvPatchField, volMesh>&
00127 ) const
00128 {
00129 return weights_;
00130 }
00131 };
00132
00133 tmp<surfaceInterpolationScheme<Type> > operator()
00134 (
00135 const GeometricField<Type, fvPatchField, volMesh>& field
00136 ) const
00137 {
00138 return tmp<surfaceInterpolationScheme<Type> >
00139 (
00140 new fieldScheme(field, weights_)
00141 );
00142 }
00143 };
00144
00145
00146
00147
00148 }
00149
00150
00151
00152
00153
00154 #define makeLimitedMultivariateSurfaceInterpolationScheme(SS, LIMITER) \
00155 typedef multivariateScheme<scalar, LimitedScheme<scalar, LIMITER<NVDTVD>, limitFuncs::magSqr> > \
00156 multivariateScheme##LIMITER_; \
00157 defineTemplateTypeNameAndDebugWithName(multivariateScheme##LIMITER_, #SS, 0); \
00158 \
00159 multivariateSurfaceInterpolationScheme<scalar>::addIstreamConstructorToTable \
00160 <multivariateScheme<scalar, LimitedScheme<scalar, LIMITER<NVDTVD>, limitFuncs::magSqr> > > \
00161 addMultivariate##SS##ConstructorToTable_;
00162
00163
00164 #define makeLLimitedMultivariateSurfaceInterpolationScheme(SS, LLIMITER, LIMITER, NVDTVD, LIMFUNC) \
00165 typedef multivariateScheme<scalar, LimitedScheme<scalar, LLIMITER<LIMITER<NVDTVD> >, limitFuncs::LIMFUNC> > \
00166 multivariateScheme##LLIMITER##LIMITER##NVDTVD##LIMFUNC##_; \
00167 defineTemplateTypeNameAndDebugWithName(multivariateScheme##LLIMITER##LIMITER##NVDTVD##LIMFUNC##_, #SS, 0); \
00168 \
00169 multivariateSurfaceInterpolationScheme<scalar>::addIstreamConstructorToTable \
00170 <multivariateScheme<scalar, LimitedScheme<scalar, LLIMITER<LIMITER<NVDTVD> >, limitFuncs::LIMFUNC> > > \
00171 addMultivariate##SS##ConstructorToTable_;
00172
00173
00174
00175 #ifdef NoRepository
00176 # include "multivariateScheme.C"
00177 #endif
00178
00179
00180
00181 #endif
00182
00183