[krishna_kant_IITH]

Hello Foamers

I am on my quest of learning openfoam trying to solve laplace equation on my own without using laplacian function.

I am trying the following

#include "fvCFD.H"
#include "simpleControl.H"

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
scalar calculateTemperature(scalar x,int id);

int main(int argc, char *argv[])
{
#include "setRootCase.H"
#include "createTime.H"
#include "createMesh.H"

simpleControl simple(mesh);

#include "createFields.H"
//#include "initContinuityErrs.H"

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

Info<< "\nStarting time loop\n" << endl;
scalar x[1]={0};
scalar y[1]={0};
scalar z[1]={0};

while (simple.loop())
{
Info<< "Time = " << runTime.timeName() << nl << endl;

while (simple.correctNonOrthogonal())
{
gradT=fvc::grad(T);
for (label cellid= 0; cellid <mesh.C().size(); cellid++)
{
for (label faceid= 0; faceid <mesh.owner().size(); faceid++)
{
x[0]=mesh.Sf()[faceid].x();
y[0]=mesh.Sf()[faceid].y();
z[0]=mesh.Sf()[faceid].z();
T[cellid]=T[cellid]-runTime.deltaTValue()*1.00*(gradT.component(0)*x[0]+gradT.component(1)*y[0]+gradT.component(2)*z[0]);
}
}
}

runTime.write();

Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
<< " ClockTime = " << runTime.elapsedClockTime() << " s"
<< nl << endl;
}

Info<< "End\n" << endl;

return 0;
}

I am getting this error :
error: cannot convert ‘Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> >’ to ‘double’ in assignment
for the line T[cellid] ...