Re: [deal.II] Looking for clarification on a few places on Step-44

his part of the process >>> would be great. * >>> >>> So this is more easy to explain. We specifically choose discontinuous >>> shape function to discretise these fields. As there are no interface/flux >>> contributions, all local element contributions f

Re: [deal.II] Looking for clarification on a few places on Step-44

diagonal, as it is a field that couples with itself. As for > the coupling matrices K_{pJ} and K_{Jp}, they are block-diagonal because we > chose *exactly* the same discretisation for both fields (i.e. the shape > functions and polynomial order match). > > Does that make s

[deal.II] Looking for clarification on a few places on Step-44

Hi, Step-44 has a lot going on and really sets you up to tackle a variety of real world problems. There is a significant jump in complexity between this tutorial and steps 8,17 & 18. I have been reading the references and I am about there but need few clarifications for why things are the

Re: [deal.II] Issue linking GMSH when installing & building dealii 9.4

las support has to be enabled in Gmsh" * > relates to CMake options when building gmsh from source, e.g., > https://gmsh.info/src/gmsh-4.10.5-source.tgz. > I would try if building gmsh yourself with those dependencies/options > resolves the issue. > > Best, > Daniel

[deal.II] Issue linking GMSH when installing & building dealii 9.4

Hi all, I was upgrading the 9.4 and I needed to include GMSH for a project I am working on so I can use the API in my dealii code. I am getting the following error despite getting the all clear from CMAKE. How can I fix this , are there version limits on GMSH? I did notice this line* "netgen,

[deal.II] Reference text for some of the theory with linear systems and block matricies

Hi all, I sort of crashed into a knowledge barrier in my learnings with the library, more specifically with Steps 20-22, which introduce block matrices and crafty manipulations to get the system ready for the CG solver. I am sort of lost with these tutorials and I think it is due to gaps in

[deal.II] Proper use of fe_values.get_function_values to get values from a vector

Hi I am again tinkering with my favorite tutorial step 23 and I am have a question again on how to move forward. What I am attempting to do now is redo this setup but rather than using helper functions to create the system matrix and RHS, I am trying to build the system by looping over the

Re: [deal.II] Wave front not propagating through mesh modification on step 23, is Affine constraints wiping out my solution?

Matt On Sat, 30 Jul 2022 at 22:47, Wolfgang Bangerth wrote: > On 7/30/22 15:45, Matthew Rich wrote: > > > > I am trying a different time integration approach using the generalized > alpha > > method rather then reposing with a velocity variable. I am also > switchin

Re: [deal.II] Wave front not propagating through mesh modification on step 23, is Affine constraints wiping out my solution?

. Maybe this is my problem that my assumption is incorrect? Anyways thanks for the correspondence. Matt On Sun, 31 Jul 2022 at 17:44, Wolfgang Bangerth wrote: > On 7/30/22 21:15, Matthew Rich wrote: > > > > However, I am not sure how to slowly morph the Step 23 tutoria

Re: [deal.II] Wave front not propagating through mesh modification on step 23, is Affine constraints wiping out my solution?

was hoping there was an obvious error in what I am doing since the solution being just zero is pretty stark. I will keep on looking at this. On Sat, 30 Jul 2022 at 22:47, Wolfgang Bangerth wrote: > On 7/30/22 15:45, Matthew Rich wrote: > > > > I am trying a different time integration

[deal.II] Wave front not propagating through mesh modification on step 23, is Affine constraints wiping out my solution?

Hi I am trying a different time integration approach using the generalized alpha method rather then reposing with a velocity variable. I am also switching to the BC handled by affine constraints. I am not sure why the wave wont move through the mesh. I just see the BC at the edge which

[deal.II] Not getting line numbers on stack trace

Hi all, I am not sure if this is a dealii issue or more a gcc issue but I am not getting line numbers on the stack traces, which makes debugging errors hard since I am a raw beginner. So for example I get the following An error

[deal.II] Step 23 - Theta method of time discretization when 2nd order time derivative

Hi My question is in relation to step 23. I am dipping my toes in the water with structural mechanics problems and I would like to avoid the quasistatic simplification. Now the example add velocity to reduce the order of the time derivative and proceeds with the theta method. I wanted to

[deal.II] Re: Program compiles does not finish linking cannot resolve 'void dealii::SparseMatrix::vmult *snip*

ok disregard this message. I found the error. On Thursday, 14 April 2022 at 11:25:08 UTC-4 Matthew Rich wrote: > Hi all, > > Trucking through a problem and hit a snag with linking > > /home/mjrich/dealii/proj/proj.cc:294: error: undefined reference to 'void > dealii::S

[deal.II] Program compiles does not finish linking cannot resolve 'void dealii::SparseMatrix::vmult *snip*

Hi all, Trucking through a problem and hit a snag with linking /home/mjrich/dealii/proj/proj.cc:294: error: undefined reference to 'void dealii::SparseMatrix::vmult, dealii::PackagedOperation > >(dealii::Vector&, dealii::PackagedOperation > const&) const' collect2: error: ld returned 1 exit

[deal.II] How to use create_mass_matrix tools with FESystem

Hi all, I am trying to use the helper tools to assemble a system for a vector valued problem. I am looking to have 2 or 3 displacement DoFs. So I use FESystem fe; and initialize with fe(FE_Q (1), dim) So my code

[deal.II] Implementing a fiber reinforcement in a biphasic model

Hi all, I am trying to add some anisotropy to my model via fibers. These fibers would only be active in tension and have different orientations depending on spatial location in the model. My plan of attack was to simply store and extra value at the quadrature points that had a direction

Re: [deal.II] Re: Compile Error: no type named ‘active_cell_iterator’

If some one would guide me on how to submit the patch, I am willing to give it a go! On Monday, 12 July 2021 at 09:54:32 UTC-4 Wolfgang Bangerth wrote: > On 7/12/21 7:28 AM, Matthew Rich wrote: > > > > At one time did dataout have an active_cell_iterator property. The code &g

[deal.II] Re: Compile Error: no type named ‘active_cell_iterator’

OK, At one time did dataout have an active_cell_iterator property. The code gallery stuff may not get the same upkeep as the tutorials and maybe the library updated and broke this piece of the code. I see both typename DoFHandler::active_cell_iterator and DataOut::active_cell_iterator,

[deal.II] Compile Error: no type named ‘active_cell_iterator’

Hi all, Once again perplexed by the c++ compiler. I am getting the following error when trying to compile a program from the code gallery. prog1.cc:In instantiation of ‘typename dealii::DataOut::cell_iterator Prog1::FilteredDataOut::first_cell() [with int dim = 3; DH = dealii::DoFHandler<3,

Re: [deal.II] New location for sacado.hpp

lag it does not work On Thursday, 8 July 2021 at 13:40:52 UTC-4 Matthew Rich wrote: > So I thought I was in good shape. I had cmake run with the proper flags > and everything and built deal.ii > > But. > > I tried to build step-33 and got the following > > -- Using

Re: [deal.II] New location for sacado.hpp

54 > > I don’t think that there’s an Ubuntu package for that yet, but it’s not > too hard to install by hand. You can take a look at their GitHub page for > build instructions. At the bottom of the page there is their "Recommended > options to build”, which you might want to take not

Re: [deal.II] New location for sacado.hpp

ow for > sure that there’s still something wrong. From these lines > > https://github.com/dealii/dealii/blob/master/cmake/configure/configure_2_trilinos.cmake#L58-L68 > there should also be some messages that indicate exactly which modules it > could not find. > > I hope that this helps

Re: [deal.II] New location for sacado.hpp

So I tried running cmake again to get trilinos picked up. It does not seem to be detected. I even try supplying the directory with cmake -/usr/local ./ -DTRILINOS_DIR=/usr/include/trilinos But it still does not get detected. Should I build trilinos from source as well? I attached another

Re: [deal.II] New location for sacado.hpp

mplemented > inside the library that make using Sacado a bit more tractable. > > Best, > Jean-Paul > > > On 7. Jul 2021, at 17:41, Matthew Rich wrote: > > Hi all, > > I am a bit confused on step-33 which uses the differentiation of sacado. > I installed the la

[deal.II] New location for sacado.hpp

Hi all, I am a bit confused on step-33 which uses the differentiation of sacado. I installed the latest version. 9.3 which has trilinios in it, and therefore sacado. Browsing that directory I see other header files for the package but not sacado.hpp The tutorial has the line #include Was this

Re: [deal.II] Install and running deal.II with latest debian package

> cd step-1 > > cmake . > > make > > make run > > > > Best, > > Matthias > > > > > > > > On Wed, Jan 13, 2021, at 23:22 CST, Matthew Rich > wrote: > > > >> Hi, > >> > >> I am new to deal.II and I

[deal.II] Install and running deal.II with latest debian package

Hi, I am new to deal.II and I wanted to mess with the tutorials. I am running the mint LDE distro of linux and just installed the package. per the website. On step one when you are just creating a mesh, I am getting an error regrading Sacado.h which appears to be a solver package external to