[ohiagree]

Hi all,
from the simple case I will present below, I guess you'll understand that I'm new to OF.

I trying to create a solver for the equation:
laplacian(p) = f
where e.g.
f = cos(pi*x)*cos(pi*y).

I'm not quite sure how to create the forcing function. In createFields.H I'm trying something like

Info<< "Reading field pressure\n" << endl;
volScalarField p(
	IOobject(
		"p",
		runTime.timeName(),
		mesh,
		IOobject::MUST_READ,
		IOobject::AUTO_WRITE
	),
	mesh
);

Info<< "Reading field forcing\n" << endl;
volScalarField f(
	IOobject(
	"f",
	runTime.timeName(),
	mesh,
	IOobject::MUST_READ,
	IOobject::AUTO_WRITE
	),
	mesh
);
	
	volScalarField yPos = mesh.C().component(vector::Y);
	volScalarField xPos = mesh.C().component(vector::X); 
	
	forAll(f, counter) 
	{ 
	f[counter] = foam::cos(3.14*xPos[counter])*cos(3.14*yPos[counter]); 
	}

And also, an extract from poissonCustom.C looks something like

#include "fvCFD.H"

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
// Main program:

int main(int argc, char *argv[])
{
    #include "setRootCase.H"
    #include "createTime.H"
	#include "createMesh.H"
	#include "createFields.H"
	// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
	
	solve ( fvm::laplacian(p) == f );
	
	runTime++;
	p.write();
	f.write();

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

    Info<< "\nEnd\n" << endl;
    return 0;
}

This creates the following errors when compiling the solver

./createFields.H:26:17: error: conversion from 'tmp<GeometricField<cmptType, fvPatchField, Foam::volMesh> >' to
      'volScalarField' (aka 'GeometricField<scalar, fvPatchField, Foam::volMesh>') is ambiguous
        volScalarField yPos = mesh.C().component(vector::Y); 
                       ^      ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

Any ideas?

Best Regards,
Bosse