I believe all I had to do was upon initializing the MatrixFree operator I needed to include a vector that specified the "selected_row_blocks" which allows you to choose which component to use in the underlying matrix free object by including a vector with the index of the dof_handlers and constriants you choose to use. The n_components function I used still returned that the object had 2 components but the solving went through as if there were only one.
Thanks for your patience, Sean Johnson On Wednesday, July 5, 2023 at 3:03:29 PM UTC-6 Sean Johnson wrote: > Okay I think I have implemented what I could of step-75. I currently have > deal.II compiled without Trilinos and am trying to keep it that way because > it is compiled without Trilinos on the high powered computer we will do > runs on later. That can change but would cause a headache for someone else > so I am just trying to avoid it. > > So I created the operator around a MatrixFree object. It allowed me the > ability to control m() for the matrix. I still have a problem and I went > back to my old version and this was a problem in both of them. Since in the > reinit function of the MatrixFree object is given two DoFHandlers, I have a > system matrix that expects a two component solution vector and rhs. The > solution vector and rhs are scalars with just one component. I am handing > both DoFHandlers over just because the computation of the matrix elements > depend on some variables I have calculated elsewhere with another > DoFHandler. > > Thanks for your time so far. I am still actively working on how to address > this. If you have any suggestions I am all ears. > I will also work on a minimal example that causes the error. > > Thanks for your time again, > Sean Johnson > > On Tuesday, June 27, 2023 at 3:19:10 PM UTC-6 Sean Johnson wrote: > >> I haven't done step-75 but I will look into it today and take a shot at >> it. There is nothing special about the way I did it besides that step-37 >> came earlier so I had already done it and was familiar with that method. >> Thank you for your time and suggestion. I'll give an update probably in a >> day or so after learning and trying to do step-75. >> >> Big picture goal is we are doing a gravity waves simulation where they >> implode and form a black hole. Numerically, it has a few variables that >> have equations for evolving forward in time and then like three that don't >> but do have elliptic equations that if we know all the variables that have >> time equations we can then solve them on each time step. Thus the need for >> a DG DoFHandler for the time equations and a normal CG DoFHandler for the >> elliptic equation variables. >> >> Thanks again so much for your time, >> Sean Johnson >> >> On Tuesday, June 27, 2023 at 2:02:59 PM UTC-6 peterr...@gmail.com wrote: >> >>> 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/68b8eb4f-866c-4105-9c47-7e8dd6ea481cn%40googlegroups.com.
