[bmikuz]

Hi!

I need a simple operator which works on a tensor and calculates another tensor with squared components (a_{ij} -> (a_{ij})^2). Is there already an operator, which does that in OF?

I've tried an outer product (gradU*gradU):

tmp<volTensorField> S2(const tmp<volTensorField>& gradU) const
        {
        return (gradU * gradU);
        }

but it seems not to be defined on tensors since compiling produces:

error: no match for ‘operator*’ in ‘gradU * gradU’

Than I tried also sqr(gradU), but compiling says

error: no matching function for call to ‘sqr(const Foam::tmp<Foam::GeometricField<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh> >&)’

Any help would be appreciated. THX!

Blaz

[bmikuz]

I may not be clear enough. I need to create a new volTensorField, which tensor's elements are squared tensor elements of gradU field (gradU is also a volTensorField typedef):

PHP Code:

gradU[cellI] =  ( a, b, c, d, e, f, g, h, i ) 
newTensor[cellI] = ( a^2, b^2, c^2, d^2, e^2, f^2, g^2, h^2, i^2 ) 


Now, I realized that there is no existent operation on volTensorField, which could directly produce such a result.

So, my idea is to define a new operation among volTensorField objects, which would directly produce required new volTensorField. Is this easy to implement or not? What do you suggest?

Thanks!
Blaz

[gu1]

Originally Posted by bmikuz

I may not be clear enough. I need to create a new volTensorField, which tensor's elements are squared tensor elements of gradU field (gradU is also a volTensorField typedef):

PHP Code:

gradU[cellI] =  ( a, b, c, d, e, f, g, h, i ) 
newTensor[cellI] = ( a^2, b^2, c^2, d^2, e^2, f^2, g^2, h^2, i^2 ) 


Now, I realized that there is no existent operation on volTensorField, which could directly produce such a result.

So, my idea is to define a new operation among volTensorField objects, which would directly produce required new volTensorField. Is this easy to implement or not? What do you suggest?

Thanks!
Blaz

Is this possible to do?

[zhangyan]

cmptMultiply(gradU,gradU)?