[stevenvanharen]

Hi does anyone have an idea how to create a volTensorField from volScalarField.

For instance I have volScalarField a and want to create the following volTensorField:

random example:

[a , 0 ,0,
a, 0 ,0,
0, 0, a]

How do I program this?

[boger]

You could always set each value explicitly; e.g.

    forAll(T,cellI)
    {
      T[cellI].xx() = a[cellI];
    }

but in your case, what about something like this:

    tensor select(1,0,0,1,0,0,0,0,1);
    T = a*select;

or more directly

T = a*tensor(1,0,0,1,0,0,0,0,1);

In all cases, a is a volScalarField and T a volTensorField.

[stevenvanharen]

Hi David,

no I solved it like this:

volTensorField T
(
IOobject
(
"Transform",
runTime.timeName(),
mesh
//IOobject::NO_READ
//IOobject::NO_WRITE
),
mesh,
tensor(0,0,0,0,0,0,0,0,0)
);
T=tensor(0,0,0,0,0,0,1,0,0)+(tensor(0,-1,0,0,0,0,0,0,0)*sin(theta))+(tensor(0,0,1,0,0,0,0 ,0,0)*cos(theta))+(tensor(0,0,0,0,-1,0,0,0,0)*cos(theta))+(tensor(0,0,0,0,0,-1,0,0,0)*sin(theta));
T.write();

Where theta is a volScalarField.

So there really is no constructor which can handle volScalarFields as an argument I guess.

P.S. How do I post code like you? Couldn't find how to do it but it looks more nice.

[olesen]

Originally Posted by stevenvanharen

...
P.S. How do I post code like you? Couldn't find how to do it but it looks more nice.

If you are posting with the WYSIWYG forum editor, use the '#' from the toolbar. This wraps "[ CODE ]" ... "[/ CODE ]" tags around the selected text, but of course without the spaces that I've need to add in order to post it.

[stevenvanharen]

Ok, thanks!

Cool