small常数

Samuel-Tu

自己定义了一个类，里面有个除法，由于担心分母太小，因此在除法之前想判别分母是否小于一个很小的数，在OpenFOAM的源代码中有如下例子：

if ((d1 - d3) < rootVSmall)

或者

     epsilonMin_
     (
         dimensioned<scalar>::lookupOrAddToDict
         (
             "epsilonMin",
             RASDict_,
             kMin_.dimensions()/dimTime,
             small
         )
     ),

我自己写成

if(d<small)或者if(d<rootVSmall)

就报错未定义small或者rootVSmall。请问如何引用这两个小量呢？

李东岳

直接用就行

Samuel-Tu

@东岳东岳老师，你是说直接除以这个数而不判断吗?这个例子是在算一个点到一个线段的距离，如果分母很小，说明线段两个端点过近，所以直接算点到端点的距离就可以。所以我想判断一下分母是否小于small。不过编译错误，提醒我未定义small。我后来直接用了10e-10了。不过不知道默认的small是多大。。

李东岳

@Samuel-Tu 我记得是1e-16，是大写的SMALL，我的意思你直接用SMALL，不需要写什么include<small>之类

Samuel-Tu

@东岳感谢，大写SMALL成功了。我还输出了一下SMALL和GREAT。输出了下是1e-15和1e+15。

星星星星晴

@Samuel-Tu ROOTVSMALL

Samuel-Tu

@星星星星晴谢谢，大写的SMALL编译成功了。ROOTVSMALL还没试过。

thegame

https://www.cfd-online.com/Forums/openfoam-solving/120990-small-great-rootvsmall-what.html
里面写了small常数还有这些以及大小

static const doubleScalar doubleScalarGREAT = 1.0e+15;
static const doubleScalar doubleScalarVGREAT = 1.0e+300;
static const doubleScalar doubleScalarROOTVGREAT = 1.0e+150;
static const doubleScalar doubleScalarSMALL = 1.0e-15;
static const doubleScalar doubleScalarVSMALL = 1.0e-300;
static const doubleScalar doubleScalarROOTVSMALL = 1.0e-150;

以及在OpenFOAM-7数值改变了(https://github.com/OpenFOAM/OpenFOAM-7/commit/c902c7a396c43ff733b33294f79acab03d7b92d8#diff-118c1e97b96036fb1f9d7f1ee5912713)

测试了一下OpenFOAM-7的VSMALL结果：

CFD中文网

small常数