cmptMag(g.value()) and mag(g.value())

sharonyue · March 6, 2016, 04:26

Hi guys,

In readgh.H, we have:

Code:

dimensionedScalar ghRef
    (
        mag(g.value()) > SMALL
      ? g & (cmptMag(g.value())/mag(g.value()))*hRef
      : dimensionedScalar("ghRef", g.dimensions()*dimLength, 0)
    );

g.value() is a scalar. So in my opinion, cmptMag(g.value()) and mag(g.value()) is the same. Why do we need this operation?

Code:

inline Scalar cmptMag(const Scalar s)
  277 {
  278     return mag(s);
  279 }

alexeym · March 6, 2016, 09:12

Hi,

In general g is uniformDimensionedVectorField, so g.value() is vector (http://foam.sourceforge.net/docs/cpp/a02728.html).

March 6, 2016, 04:26	cmptMag(g.value()) and mag(g.value())	#1
sharonyue Senior Member Dongyue Li Join Date: Jun 2012 Location: Beijing, China Posts: 838 Rep Power: 17	Hi guys, In readgh.H, we have: Code: dimensionedScalar ghRef ( mag(g.value()) > SMALL ? g & (cmptMag(g.value())/mag(g.value()))hRef : dimensionedScalar("ghRef", g.dimensions()dimLength, 0) ); g.value() is a scalar. So in my opinion, cmptMag(g.value()) and mag(g.value()) is the same. Why do we need this operation? Code: inline Scalar cmptMag(const Scalar s) 277 { 278 return mag(s); 279 } __________________ My OpenFOAM algorithm website: http://dyfluid.com By far the largest Chinese CFD-based forum: http://www.cfd-china.com/category/6/openfoam We provide lots of clusters to Chinese customers, and we are considering to do business overseas: http://dyfluid.com/DMCmodel.html

March 6, 2016, 09:12		#2
alexeym Senior Member Alexey Matveichev Join Date: Aug 2011 Location: Nancy, France Posts: 1,930 Rep Power: 38	Hi, In general g is uniformDimensionedVectorField, so g.value() is vector (http://foam.sourceforge.net/docs/cpp/a02728.html). sharonyue, wwzhao and Milica like this.