Dear NGSolve team,

We have voxelized diffusion tensor imaging (DTI) data that we want to use in a CoefficientFunction as material properties.

We already used the VoxelCoefficient function to load voxelized MRI data (i.e., only one scalar value per voxel).

However, we could not figure out what the best way is to load tensor data (i.e., one tensor per voxel).

Could you please suggest us a good solution or share your experience?

We thought about loading the tensor component-wise into a VoxelCoefficient and then turning it into a matrix CoefficientFunction.

Thanks for your help and best regards,

Julius