[Mapping] Map Vector/Matrix Variables

[ddiezrod]

[Mapping] Map Vector/Matrix Variables

[ddiezrod]

@matekelemen I think I'm going to start with the second option as there is no harm in it, but we can maybe try to think about something else for the future.

[jcotela]

Just a note: the mapper works with variable components (which are Variable<double> internally) so you can already kind-of support other array types. I have already used that in the past to map array_1d<double, 6> variables one component at a time. You just need to loop over components on your side when calling the map function (which is the same the mapper already does for array-3 anyway)

[matekelemen]

I'd go with the second option (extend the interface with array_1d<double,6>).

I think we could also go with a limited support for "dynamic" arrays by extending the interface with bounded_array. Sure, it's a lot of overhead but it's definitely safer (and maybe even faster) than asking the size of each node's Vector (or, god forbid, assuming all nodes have Vectors of the same size).

[ddiezrod]

@jcotela hmm you are right, I think I will just use it like that, thanks for the suggestion!

[ddiezrod]

I'll close this for the moment as @jcotela solution with the workaround of creating an array_1d6 intermediary variable is working fine for me. Thanks all!

[philbucher]

Yep as @jcotela says I am anyway doing it component wise: https://github.com/KratosMultiphysics/Kratos/blob/master/applications/MappingApplication/custom_mappers/interpolative_mapper_base.h#L473-L487

If you want the interface could be extended to array_1d6, I dont see a problem as long as it is fixed size
Dynamic size is different though, not sure if it is worth the pain and hassle ...