[GerhardHolzinger]

Hi,

I want to create a post-processing utility, which traverses all time steps and compares the data with the data of all other time steps. So I end up with two nested loops traversing all time steps.

So, in pseudo code, it looks like this:

    // go over all time steps
    forAll(timeDirs, timeI)
    {
        runTime.setTime(timeDirs[timeI], timeI);
        
        // go over all time steps
        forAll(timeDirs, timeJ)
        {
            secTime.setTime(timeDirs[timeJ], timeJ);
            
            if (timeI == timeJ)
            {
                // easy
                recurrenceMatrixAlpha[timeI][timeJ] = 1;
            }
            else
            {
                // do stuff
            }
        }
    }

In the "// do stuff" section, I read the fields at the two different times. Right now, I call the constructor, like this:

volScalarField alpha1
(
    IOobject
    (
        "alpha.air",
        runTime.timeName(),
        mesh,
        IOobject::MUST_READ,
        IOobject::NO_WRITE
    ),
    mesh
);
volScalarField alpha2
(
    IOobject
    (
        "alpha.air",
        secTime.timeName(),
        mesh,
        IOobject::MUST_READ,
        IOobject::NO_WRITE
    ),
    mesh
);

I assume, that this might be done faster.

Is there a way to just read the field at the current time, like read()?