[mbookin]

Hello!

I run my solver which uses field.oldTime().oldTime().oldTime() function. When I start from initial time the field_0_0_0 updates in objectRegistry and every thing is ok. But then I continue calculation from some time the field_0_0_0 in object registry is not equal to field_0_0_0 its self. In objectRegistry it is equal to initial value read from time directory.

So the question is how to update a field_0_0_0 in its objectRegistry manually?

// create or read T
// mesh is objectRegistry for T
while (runTime.run())
    {
    //  update T
    //  use T.oldTime().oldTime().oldTime()
    Info<<T.oldTime().oldTime().oldTime()<<endl;
    //  and 
    Info<<mesh.lookupObject<volScalarField>("T_0_0_0")<<endl;
    runtime.write()
    }

The code return equal value when I start from the beginning and return different value when I continue calculation. The whole solver is complicated to show it here and I haven't done any simple one to reproduce this problem hoping that someone can answer my question without reproduction. I use OF-2.3.x on Linux based system (SL-6.0)

[jmf]

Hi

You might want to read my answer #2 to a similar problem here:
http://www.cfd-online.com/Forums/ope...dymfoam-c.html

based on the answer #14 in :
http://www.cfd-online.com/Forums/ope...rfield-bc.html

Regards

Jean-Michel

[mbookin]

Thank you for your answer!

The post you gave, as I understood, is about how to write something that is not in registry. You use function set() to update you dictionaries. My field_0_0_0 is in the registry, but is not updated. The problem is that geometricField class has no such function like set() so I cant update it manually or may be I miss something. Did I understand you right?

[jmf]

Dear Maxim


You don't give much details...
Instead of "runtime.write()" would "T.write()" work ?

Jean-Michel

[mbookin]

I made simple testSolver and testCase for someone who wants to help me to handle this problem. PseudoCode as follows, and the code is attached below.

volScalarField T
(
     IOobject
     (
        "T",
         runTime.timeName(),
         mesh,
         IOobject::MUST_READ,
         IOobject::AUTO_WRITE
     ),
     mesh
); 
while (runTime.run())
    {
        runTime++;

        T=runTime.timeIndex();

        Pout<<tab<<"T="<<T[0]<<endl;
        Pout<<tab<<"T.()="<<T.oldTime()[0]<<tab<<"  writeOpt=  "<<T.oldTime().writeOpt()<<endl;
        Pout<<tab<<"T.()()="<<T.oldTime().oldTime()[0]<<tab<<"  writeOpt=  "<<T.oldTime().oldTime().writeOpt()<<endl;
        Pout<<tab<<"T.()()()="<<T.oldTime().oldTime().oldTime()[0]<<tab<<"  writeOpt=  "<<T.oldTime().oldTime().oldTime().writeOpt()<<endl;
        Pout<<tab<<"T.reg_0_0_0="<<T.db().lookupObject<volScalarField>("T_0_0_0")[0]<<endl;
 
        runTime.write();
}

The output is as follows:

Reading field T
 T.().().writeopt=  1
 T.().().().writeopt=  1
 T=0
 T.()=0
 T.()()=0
 T.()()()=0
Time = 1
 T=1
 T.()=0   writeOpt=  0
 T.()()=0   writeOpt=  1
 T.()()()=0   writeOpt=  1
 T.reg_0_0_0=0
Time = 2
 T=2
 T.()=1   writeOpt=  0
 T.()()=0   writeOpt=  0
 T.()()()=0   writeOpt=  1
 T.reg_0_0_0=0
Time = 3
 T=3
 T.()=2   writeOpt=  0
 T.()()=1   writeOpt=  0
 T.()()()=0   writeOpt=  1
 T.reg_0_0_0=0
Time = 4
 T=4
 T.()=3   writeOpt=  0
 T.()()=2   writeOpt=  0
 T.()()()=1   writeOpt=  1
 T.reg_0_0_0=1
Time = 5
 T=5
 T.()=4   writeOpt=  0
 T.()()=3   writeOpt=  0
 T.()()()=2   writeOpt=  1
 T.reg_0_0_0=2
End

T.oldTime().oldTime().writeOpt() is switched to AUTO_WRITE automatically. But T.oldTime().oldTime().oldTime().writeOpt() is not. So T.oldTime().oldTime().oldTime() is not written. That's why it will be updated in a wrong way when you continue. So the problem is how to switch T.oldTime().oldTime().oldTime().writeOpt() to AUTO_WRITE??
I even can write need field by T.db().lookupObject<volScalarField>("T_0_0_0").wri te(); But when I continue it is still not updated in objectRegistry, so the second continue is impossible.

Attachment 52119