[deal.II] Convective Heat Flux Boundary Conditions in deall ii?

2021-09-14 Thread Morris Jowas
Hello everyone, I am trying to solve a simple thermal problem but finding the documentation helpfully detailed yet overwhelming. My problem is on heat conduction within a cubic domain with internal heat generation. I have studied the transient heat transfer problem in step 26, in

[deal.II] What I learned from 20 years of leading open source projects

2021-09-14 Thread Wolfgang Bangerth
All: in case you wonder what is going on in the heads of the principal developers of deal.II: I will be giving a webinar tomorrow on this topic: https://bssw.io/events/webinar-what-i-learned-from-20-years-of-leading-open-source-projects I believe that you have to register (it's free!) but

[deal.II] deal.II Newsletter #182

2021-09-14 Thread 'Rene Gassmoeller' via deal.II User Group
Hello everyone! This is deal.II newsletter #182. It automatically reports recently merged features and discussions about the deal.II finite element library. ## Below you find a list of recently proposed or merged features: #12758: Rename VECTOR to VectorType (proposed by peterrum; merged)

Re: [deal.II] Re: step-29 & step-40

2021-09-14 Thread Hermes Sampedro
Thank you very much. Regards, Hermes El mar, 14 sept 2021 a las 16:28, Bruno Turcksin () escribió: > Hermes, > > You probably don't need to use reinit() in your case. You want to use > reinit() if the structure of the matrix changes with frequency. For > example, if you change the degree of

Re: [deal.II] MatrixTools::apply_boundary_values for BlockSparseMatrix> and BlockVector> entries

2021-09-14 Thread Daniel Arndt
Mariia, You need to make sure that the std::map also uses values of type std::complex. The compile error you shared shows that you use doubles. Best, Daniel Am Di., 14. Sept. 2021 um 07:30 Uhr schrieb Мария Бронзова < masianic...@gmail.com>: > Hi everyone, > > I ran into the following

Re: [deal.II] Re: step-29 & step-40

2021-09-14 Thread Bruno Turcksin
Hermes, You probably don't need to use reinit() in your case. You want to use reinit() if the structure of the matrix changes with frequency. For example, if you change the degree of finite elements you need to call reinit(). If you want to change the values but the structure matrix is

Re: [deal.II] Re: step-29 & step-40

2021-09-14 Thread Hermes Sampedro
Dear Bruno thank you very much for pointing step-17, I could solve the problem. I would like to ask the last question. I am computing step-29 in parallel for different frequencies. I have a loop for each of the frequencies as follows: make_grid(); setup_system();

Re: [deal.II] Re: step-29 & step-40

2021-09-14 Thread Bruno Turcksin
Hermes, Le mar. 14 sept. 2021 à 05:19, Hermes Sampedro a écrit : > > Should I use dealii::PETScWrappers::MPI::SparseMatrix system_matrix instead? > If so could ou please help me to with the reinit() function? I do not fully > understand how to call it. That's right, you need the matrix to be

[deal.II] MatrixTools::apply_boundary_values for BlockSparseMatrix> and BlockVector> entries

2021-09-14 Thread Мария Бронзова
Hi everyone, I ran into the following difficulty: I have a vector-valued problem with complex numbers as entries for my BlockSparseMatrix system_matrix and BlockVector system_rhs and solution vectors and I am trying to apply Dirichlet boundary conditions, using the

Re: [deal.II] Re: step-29 & step-40

2021-09-14 Thread Hermes Sampedro
Dear Bruno, thank you for your answer it helped me to solve the problem. I had some issues with Trilinos library so I decided to implement step-29 in parallel (similar to step-40) using dealii::PETSxWrappers. I get an error on the solve() function on the red line: *template* <*int* dim> *void*