← Back to team overview

dolfin team mailing list archive

Re: tensor function/vtk output

 

Kristian Oelgaard wrote:
Quoting Riesen  Patrick <priesen@xxxxxxxxxxxxxxx>:

Hi Kristian

The VTK tensor output in dolfin is still messed up because now the number of
components obtained by

  // Get number of components
  uint dim = 1;
  for (uint i = 0; i < rank; i++)
    dim *= element.value_dimension(i);

returns 16 instead of 4. The ffc code for generating value_dimension() might
also need a hack..., i'm sorry :-(

I was afraid that something like this would happen. Fixing this will require
quite a bit of work I'm afraid and then we might still run into problems.
I think it will be better to simply implement support for tensor elements, but
that will not happen in a foreseeable future.
Saving the components of the tensor individually is your best option for now.

that's ok for me. you can also revert the hack in ffc to keep things straight, if that's preferable. for dolfin, you may want to report it and notify that one needs to output each component individually as a workaround ...

best,
patrick


Kristian

patrick