[cfd@kgp]

Dear Foam experts',

here is a coded function to calculate a derived variable using a existiong variable IO object variable sMS1:alpha1. The following quoted text was added in the controlDict,

TIntegral
{
functionObjectLibs ("libutilityFunctionObjects.so");
type coded;
redirectType integral;
outputControl timeStep;
code
#{
const volScalarField& fl = mesh().lookupObject<volScalarField>("sMS1:alpha1") ;
// volScalarField magU(mag(U));

scalar volIntegral = 0;
scalar volIntegralV = 0;

forAll (fl, cellI)
{
volIntegral += fl[cellI]*mesh().V()[cellI];
volIntegralV += 1.0*mesh().V()[cellI];
}
volIntegral=volIntegral/volIntegralV;
Info<<"ANSWER: " << volIntegral << endl;
//foo.write();

#};
}

In serial running of the solver it performs and does the job. But when executed in parallel the output is always zero.

Somebody kindly help.

Thanks and regards,
cfd@kgp

[Tobi]

Dear cfd@kgp,

I am not 100% sure about the error but you can check (output) the size of your lookedUp field in serial and parallel. If both are identical, the lookUp command might be different. However, I think you should also use gSum instead of your primitive summing because gSum is the summation function for parallel stuff.

Good luck

[cfd@kgp]

Thanks Tobi!

can you please instuct me further on "check (output) the size of your lookedUp field in serial and parallel"

gsum do not work.

Using dynamicCode for functionObject TIntegral at line 104 in ".TIntegral"
Creating new library in "dynamicCode/integral/platforms/linuxGccDPInt32Opt/lib/libintegral_622d34218b4809de0bc1877891d4b79ccf99a3 25.so"
Invoking "wmake -s libso /home/ojas/MY_works_VLP/OpenFoam_learning/practice/beta_-1/dynamicCode/integral"
wmakeLnInclude: linking include files to ./lnInclude
Making dependency list for source file FilterFunctionObjectTemplate.C
Making dependency list for source file functionObjectTemplate.C
.TIntegral: In member function ‘virtual void Foam::integralFunctionObject::write()’:
.TIntegral:126:30: error: ‘gsum’ was not declared in this scope
make: *** [Make/linuxGccDPInt32Opt/functionObjectTemplate.o] Error 1


--> FOAM FATAL IO ERROR:
Failed wmake "dynamicCode/integral/platforms/linuxGccDPInt32Opt/lib/libintegral_622d34218b4809de0bc1877891d4b79ccf99a3 25.so"


file: .TIntegral from line 104 to line 109.

From function codedBase::createLibrary(..)
in file db/dynamicLibrary/codedBase/codedBase.C at line 213.

FOAM exiting

but simple "sum" operation does the job.

But problem is partially solved.

If the expresion could be written in Swak4Foam then we get a clean output file, with the method of coded function one has to extract the output from LOG file.

I found some important hint for using the varibles with colen is to define a alias. refer #13 https://www.cfd-online.com/Forums/op...variables.html

please have a look

GFlAvg
{
aliases {
sMS1:alpha1 myfl;
// thermo:cp myCp;
}
type swakExpression;
valueType cellZone;
zoneName solid; // or whatever is your zoneName
accumulations (
average
);

expression "myfl";
verbose true;
// outputControlMode timeStep;
// outputInterval 1;
}

but " aliases { sMS1:alpha1 myfl; }" is not recognised with OF3.X and the corresponding version of Swak4Foam

Regards,

[Tobi]

Hi,

if sum is working, gSum should work too because this is the function for parallel stuff. If you use sum, it should be not possible to run the stuff in parallel.

Sent from my HTC One mini using CFD Online Forum mobile app

[cfd@kgp]

Hi Tobi,

Your argument seem logical Tobi...

But sum gives correct answer in parallel too....... I think sum is the compiler based function and thus mpi is already tuned for compiler functions, it should implicitly applies the reduce command.

Regards.

[cfd@kgp]

Originally Posted by Tobi

Hi,

if sum is working, gSum should work too because this is the function for parallel stuff. If you use sum, it should be not possible to run the stuff in parallel.

Sent from my HTC One mini using CFD Online Forum mobile app

Hi Tobi!

I need your suggestion on this unique problem

the coded function which runs in OF3.01 and prints the value of "ANSWER" fails to do so in OF4.0. Here is the extract from the controlDict from one of the above directory

Also, please note that there are no errors during comiling the dynamic code in OF4.0 .

Dummy
{
functionObjectLibs ("libutilityFunctionObjects.so");
type coded;
redirectType dummy;
outputControl timeStep;
code
#{
const volScalarField& fl = mesh().lookupObject<volScalarField>("T") ;
// volScalarField magU(mag(U));

scalar volIntegral = 0;
scalar volIntegralV = 0;

forAll (fl, cellI)
{
volIntegral += fl[cellI]*mesh().V()[cellI];
volIntegralV += 1.0*mesh().V()[cellI];
}
volIntegral=volIntegral/volIntegralV;
Info<<"ANSWER: " << volIntegral << endl;
//foo.write();

#};

This is the output during dynamic compiling of the coded function,

Reading "/home/ojas/OPENFOAM_4.0_18_november/withGMSH/Ra10e8/L1000/system/controlDict.functions.Dummy" from line 135 to line 140
Using deprecated 'outputControl'
Please use 'writeControl' with 'writeInterval'
Using dynamicCode for functionObject Dummy at line 135 in "/home/ojas/OPENFOAM_4.0_18_november/withGMSH/Ra10e8/L1000/system/controlDict.functions.Dummy"
Creating new library in "dynamicCode/dummy/platforms/linux64GccDPInt32Opt/lib/libdummy_a2d2bd98ff2a92c53bf305564b91683acd9e80aa. so"
Invoking "wmake -s libso /home/ojas/OPENFOAM_4.0_18_november/withGMSH/Ra10e8/L1000/dynamicCode/dummy"
wmakeLnInclude: linking include files to ./lnInclude
Making dependency list for source file functionObjectTemplate.C
'/home/ojas/OPENFOAM_4.0_18_november/withGMSH/Ra10e8/L1000/dynamicCode/dummy/../platforms/linux64GccDPInt32Opt/lib/libdummy_a2d2bd98ff2a92c53bf305564b91683acd9e80aa. so' is up to date.
--> FOAM Warning :
From function void Foam::timeControl::read(const Foam::dictionary&)
in file db/functionObjects/timeControl/timeControl.C at line 89
Reading "/home/ojas/OPENFOAM_4.0_18_november/withGMSH/Ra10e8/L1000/system/controlDict.functions.sh" from line 160 to line 170
Using deprecated 'outputControl'
Please use 'writeControl' with 'writeInterval'

[Tobi]

Most of the warnings you got are related to the functions dictionary in controlDict and are not related to your coded boundary condition (if I got it right)

Sent from my HTC One mini using CFD Online Forum mobile app

[cfd@kgp]

Originally Posted by Tobi

Most of the warnings you got are related to the functions dictionary in controlDict and are not related to your coded boundary condition (if I got it right)

Sent from my HTC One mini using CFD Online Forum mobile app

Thanks, but do you have any suggestion regarding how to make it run in OF4?

[kodexys]

cfd@kgp, have you discovered a fix for the function dictionary warnings? I'm not sure how to be rid of them.

[cfd@kgp]

Originally Posted by kodexys

cfd@kgp, have you discovered a fix for the function dictionary warnings? I'm not sure how to be rid of them.

Hi kodexys,

the fix is very simple just replace "( blah blah " by "{ blah blah }"

Hope this will help you. All the best!

Thanks

[elomri]

Originally Posted by cfd@kgp

the fix is very simple just replace "( blah blah " by "{ blah blah }"

could you please specify what do you mean here ? i'm facing the same warning problem