Helmholtz solve on GPUs with multicomponent vector

[GOURAB PANIGRAHI]

I am trying to reproduce dealii’s step 64 for the Helmholtz solve on GPUs and I am interested in the source(input) and destination(output) to have more than single component. To elaborate, I am interested in using GPU matrix-free capability of dealii to compute the action of Helmholtz operator on the multicomponent vector. The CUDAWrappers::FEEvaluation has “n_components” as one its template parameters which should make this possible I think but I don’t understand how to pass this multicomponent vectors to functions like vmult and cell_loop of mf_data which currently takes one component input and output vectors.

Thanks

Gourab

[Bruno Turcksin]

Gourab,

Multicomponent vector is currently not supported. Part of the infrastructure is already present as you saw but you cannot use MatrixFree on GPU with multicomponent vectors at the moment.

Best,

Bruno

[GOURAB PANIGRAHI]

Thanks a lot for the info!

Regards

Gourab