OK. The problem are these lines: https://github.com/dealii/dealii/blob/8d68b3192ed1c8b928dfce03ca1d615201859904/include/deal.II/matrix_free/operators.h#L1225-L1235. The number of rows are the sum of the number of DoFs of both attached DoFHandlers...
Personally, I would not use these matrix-free operators but create the operator myself around a MatrixFree object. Just like in step-75. Would that be possible? Out of curiosity, what is your goal? Best, PM On Tuesday, June 27, 2023 at 7:38:48 PM UTC+2 sljoh...@gmail.com wrote: > Here is step-37 with minimal changes. All that was changed was in the > creation of a DG DoFHandler and the fe_system and constrains and what not > to go along with it. Then of the functions only the LapalceProblem > setup_system() function was altered and nothing else. It was just altered > to also be reinit with the dof_handler_DG. > > Obviously the LaplaceOperator in my actual code actually uses variables > that depend on the DG DoFHandler but that was left out of this example. > > Thank you again so much for your time and for all you guys do for this > library. > > Thanks, > Sean Johnson > > On Tuesday, June 27, 2023 at 10:58:53 AM UTC-6 Sean Johnson wrote: > >> Yes of course. I can get one with minimal changes made to step-37. Full >> code do you want the entirety of the changed step-37 code or just the >> LaplaceProblem class or just the setup_system() function? >> >> On Tuesday, June 27, 2023 at 10:53:04 AM UTC-6 peterr...@gmail.com wrote: >> >>> I'll try to take a look at the issue tomorrow. Could you attach the full >>> code. I have a feeling what might be an issue. >>> >>> Peter >>> >>> On Tue, Jun 27, 2023, 18:18 Sean Johnson <sljoh...@gmail.com> wrote: >>> >>>> Okay. Thank you. >>>> >>>> On Tuesday, June 27, 2023 at 10:02:16 AM UTC-6 Wolfgang Bangerth wrote: >>>> >>>>> On 6/27/23 09:53, Sean Johnson wrote: >>>>> > I keep seeing other questions getting at least one response. Is >>>>> there >>>>> > something else that I should provide in addition to the question or >>>>> that is >>>>> > missing or is there somewhere this question has already been >>>>> addressed that I >>>>> > cannot find? I am really all ears >>>>> >>>>> Sean, it's all about who has the expertise to answer questions and who >>>>> has the >>>>> time. I don't know anything about matrix-free stuff, so your questions >>>>> require >>>>> someone else to have the time and expertise. >>>>> >>>>> Best >>>>> W. >>>>> >>>>> -- >>>>> ------------------------------------------------------------------------ >>>>> >>>>> Wolfgang Bangerth email: bang...@colostate.edu >>>>> www: http://www.math.colostate.edu/~bangerth/ >>>>> >>>>> >>>>> -- >>>> The deal.II project is located at http://www.dealii.org/ >>>> For mailing list/forum options, see >>>> https://groups.google.com/d/forum/dealii?hl=en >>>> --- >>>> You received this message because you are subscribed to the Google >>>> Groups "deal.II User Group" group. >>>> To unsubscribe from this group and stop receiving emails from it, send >>>> an email to dealii+un...@googlegroups.com. >>>> To view this discussion on the web visit >>>> https://groups.google.com/d/msgid/dealii/dab3cc77-bd64-445a-82f5-da55b00a812en%40googlegroups.com >>>> >>>> <https://groups.google.com/d/msgid/dealii/dab3cc77-bd64-445a-82f5-da55b00a812en%40googlegroups.com?utm_medium=email&utm_source=footer> >>>> . >>>> >>> -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/a8ec30df-c93d-46c6-a971-cbbe3e52b299n%40googlegroups.com.
