Hello Tariq,
thanks for the clarification. Which rank has the tensor-value of your
function, 1 or 2?
Best Regards,
Markus
Am 16.05.10 16:17, schrieb tarbaig:
Hi Markus,
I took my time to think about what I want to do so here is my problem
and my Idea.
Of course you are right to say that a tensor in dealii as object living
at a point. The term tensor-field would be better. In deal I think the
most common tensor field is that constituted by the "Gradient" of a
Vectorvalued Function.
If I want to solve an PDE depending on a tensor-field using e.g. newton
iterations and do not want to compute the linearization by hand , but
rather using sacado, I will have to do something similar to stepp 33 but
with tensors. This gets ugly very fast, so I thought it would be nice to
have something like a class derivable_tensor that can be used for this
purpose.
Does anyone have an idea how to implement something like that ?
The sacado documentation is a little bit sparse and I am not sure how to
implement methods that do not destroy the property of being derivable.
The ugly C way to generate something that does the trick would be
preprocessor magic by macros but of course this is undesirable.
Greetings tariq
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
