[m.maneshi]

Dear all
i
in OpenFOAM/OpenFOAM-1.6/src/lagrangian/molecularDynamics/molecule at line 1158 we have the constrains and thermostats, yet i cant figure out if this member function is working or not, here it is the function

void Foam::moleculeCloud::applyConstraintsAndThermostats
(
    const scalar targetTemperature,
    const scalar measuredTemperature
)
{
    scalar temperatureCorrectionFactor =
        sqrt(targetTemperature/measuredTemperature);

    Info<< "----------------------------------------" << nl
        << "Temperature equilibration" << nl
        << "Target temperature = "
        << targetTemperature << nl
        << "Measured temperature = "
        << measuredTemperature << nl
        << "Temperature correction factor = "
        << temperatureCorrectionFactor << nl
        << "----------------------------------------"
        << endl;

    iterator mol(this->begin());

    for (mol = this->begin(); mol != this->end(); ++mol)
    {
        mol().v() *= temperatureCorrectionFactor;

        mol().pi() *= temperatureCorrectionFactor;
    }
}

must not the function print the two lines every time it is called ??

Thanks

[l_r_mcglashan]

You mean the Info statements? Are you running it in parallel or in serial?

[m.maneshi]

Yes i exactly mean the Info statement, i tried in both serial and parralel, no diffrence.

[l_r_mcglashan]

It won't be called then.

If you run in parallel, I think Info will only output for the master processor (correct me if wrong)? If you want it to output the results for each processor you have to use Pout instead of Info.

[m.maneshi]

So how can i call it in the molecule.c ?