[deal.II] Re: Development activity over the past year

2016-07-31 Thread Wolfgang Bangerth
I'm writing my annual report to the NSF and wanted to share a few statistics I thought were interesting: * 40,000 new lines of code in deal.II (excluding testsuite) over the past year * 77 new testcases * 4 new tutorial programs * 1540 emails to dealii@googlegroups.com between January to July a

[deal.II] Development activity over the past year

2016-07-31 Thread Wolfgang Bangerth
All, I'm writing my annual report to the NSF and wanted to share a few statistics I thought were interesting: * 40,000 new lines of code in deal.II (excluding testsuite) over the past year * 77 new testcases * 4 new tutorial programs * 1540 emails to dealii@googlegroups.com between January to

[deal.II] Re: time-dependent boundary conditions for distributed computations

2016-07-31 Thread Daniel Arndt
Marek, You would just set up your ConstraintMatrix anew. This is also the way it is done in step-44. Best, Daniel -- 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 becau

[deal.II] time-dependent boundary conditions for distributed computations

2016-07-31 Thread Marek Čapek
Hello, is it possible to compute in distributed manner like in Step-40 and to impose different boundary condition in each time step. In the step-40 the constraint matrix is set in front constraints.clear

Re: [deal.II] HP with parallel::shared::Triangulation

2016-07-31 Thread Pete Griffin
Thanks, Wolfgang Pete Griffin On Sunday, July 31, 2016 at 2:53:03 PM UTC-4, bangerth wrote: > > On 07/31/2016 12:06 PM, Pete Griffin wrote: > > > > I get an error from an assertion in the call to: > >dof_handler.distribute_dofs (fe_collection); > > The triangulation is a parallel::shared:

Re: [deal.II] HP with parallel::shared::Triangulation

2016-07-31 Thread Wolfgang Bangerth
On 07/31/2016 12:06 PM, Pete Griffin wrote: I get an error from an assertion in the call to: dof_handler.distribute_dofs (fe_collection); The triangulation is a parallel::shared::Triangulation. My C++ is not as good as it needs to be. It appears as if the exception (below) means that HP is n

Re: [deal.II] Re: A question about velocity correction in Step-35

2016-07-31 Thread Wolfgang Bangerth
On 07/31/2016 05:16 AM, Daniel Arndt wrote: I will create a PR next week including some more references to the analysis for that scheme. That would be much appreciated! Thanks! W. -- Wolfgang Bangerth emai

Re: [deal.II] Re: Memory usage using MPI/PETSc on single processor

2016-07-31 Thread Martin Kronbichler
Hi Pete, The difference may be that the NEW version uses a DynamicSparsityPattern while the OLD only guesses on the size with, dof_handler.max_couplings_between_dofs(). Apparently with 3d problems the function overestimates. Presently DynamicSparsityPattern is used in step-8, but, not in ste

[deal.II] HP with parallel::shared::Triangulation

2016-07-31 Thread Pete Griffin
Hi, All. I get an error from an assertion in the call to: dof_handler.distribute_dofs (fe_collection); The triangulation is a parallel::shared::Triangulation. My C++ is not as good as it needs to be. It appears as if the exception (below) means that HP is not implemented for a parallel::shar

Re: [deal.II] Re: Memory usage using MPI/PETSc on single processor

2016-07-31 Thread Pete Griffin
Hello, Wolfgang. I looked at the memory consumption of the system_matrix, the system_rhs and solution vectors. The difference appears, as Timo suggested, in the system_matrix. The difference may be that the NEW version uses a DynamicSparsityPattern while the OLD only guesses on the size with,

Re: [deal.II] Re: A question about velocity correction in Step-35

2016-07-31 Thread Daniel Arndt
Wolfgang, Interesting that this isn't even mentioned in the introduction to that > program, or anywhere else, for that matter. How did you find out? > A major part of my PhD thesis is based on these projection methods. Reading the literature this seems to be the first paper in which the rotatio

Re: [deal.II] Re: A question about velocity correction in Step-35

2016-07-31 Thread Jiaqi ZHANG
Hello Daniel, Thanks for your help again. So the weakly divergence-free velocity is not computed in Step-35. Best, Jiaqi 2016-07-30 6:12 GMT-04:00 Daniel Arndt : > Jiaqi, > > the algorithm in step-35 is the one that is proposed in > "L. Timmermans, P. Minev, and F. Van De Vosse, “An approxima

[deal.II] Re: Floquet periodic conditions and complex valued algebra

2016-07-31 Thread Daniel Garcia
Dear Daniel, Thanks a lot for your answer. I will start to use deal.ii and I will come back with more specific questions. Best, Daniel On Saturday, July 30, 2016 at 2:15:37 PM UTC+2, Daniel Arndt wrote: > > Daniel, > > 1.) You can use complex valued algebra if you split your problem into two >

[deal.II] Re: Using the solution from one problem as a boundary condition in another problem with matching mesh on the boundary

2016-07-31 Thread krei
Thanks for the response. I think this is worth a try. However, as the electric field calculation is entirely decoupled from the nonlinear currents & heating part, which need to be solved by newton's method, then wouldn't it be effective to just once calculate the field and then start the newton