Re: [deal.II] Does MatrixFree support multiple different types of Quadrature formula?

2025-02-18 Thread 'Martin Kronbichler' via deal.II User Group
Dear Chengjiang, This is a good question. This case is supported, but one needs to work slightly harder (and I admit we do not show this case). You should do as follows: std::vector> quad;     quad.push_back(QGauss<1>(fe_degree + 1));     quad.push_back(QGaussLobatto<1>(fe_degree + 2));     m

Re: [deal.II] Iterating over two DoFhandlers for a thermo-elastic problem

2025-01-05 Thread &#x27;Martin Kronbichler&#x27; via deal.II User Group
the coupling where information from the two vectors needs to be accessed simultaneously. Subramanya. सुब्रह्मण्य . On Sun, Jan 5, 2025 at 1:51 AM 'Martin Kronbichler' via deal.II User Group wrote: Dear Subramanya, The actual cells associated with the integer indices 'ce

Re: [deal.II] Iterating over two DoFhandlers for a thermo-elastic problem

2025-01-04 Thread &#x27;Martin Kronbichler&#x27; via deal.II User Group
Dear Subramanya, The actual cells associated with the integer indices 'cell' (that are batches in case of SIMD vectorization) are not necessarily the same if you create two separate MatrixFree objects and initialize each with the respective DoFHandler for the two fields. If you intend to use c