[deal.II] The type of particle properties

2019-09-11 Thread Yidong ZHAO
Hello all, >From the descriptions of Particle( https://www.dealii.org/current/doxygen/deal.II/classParticles_1_1Particle.html) and PropertyPool( https://www.dealii.org/current/doxygen/deal.II/classParticles_1_1PropertyPool.html), I know the properties of a particle should be ArrayView type. I

Re: [deal.II] Function parser for tensor function

2019-09-11 Thread Wolfgang Bangerth
On 9/6/19 1:36 AM, Konrad wrote: > > I was wondering if there is a simple way to parse input (e.g., from a > parameter file) and use it similarly to the FunctionParser class (which > works for vector valued functions) just for a TensorFunction. Nothing I would know already does that. But the Fu

[deal.II] Re: tutorial step-36 failed with PETSc error number 76 while calling a SLEPc function

2019-09-11 Thread Pragalv Karki
Dear Prof. Wolfgang, Thank you very much for your reply and insight. I will change the PETSc/SLEPc versions to 3.7.4/3.7.3 and see if that will solve the issue. On Wednesday, September 11, 2019 at 10:39:08 AM UTC-7, Pragalv Karki wrote: > > Dear deal.II community, >

Re: [deal.II] step-21 and 43

2019-09-11 Thread Wolfgang Bangerth
> I read step-21 and 43 tutorials which are closer to my task. > > I'd like to ask some guidelines on how to implement a water/oil with > polymer concentration in that model. We still have two phases, but in > the water phase there will be two components, aqueous and polymeric. > > Is there a

Re: [deal.II] tutorial step-36 failed with PETSc error number 76 while calling a SLEPc function

2019-09-11 Thread Wolfgang Bangerth
Pragalv, I've tried to run step-36 on my workstation and it works just fine. Are you running in parallel? You might also want to try a different version of PETSc/SLEPc -- I am using 3.7.4/3.7.3 for the two. I also tried to find out what error 76 actually means. I *think* that it is this error

Re: [deal.II] Question about (parallel) linear solver

2019-09-11 Thread Wolfgang Bangerth
On 9/11/19 11:03 AM, Konrad wrote: > > That indeed depends on the geometry of the domain and it is not an > obvious result. Example: Think of a 3D hyper-shell with an inner and > outer radius. There is a one-dimensional space of non-trivial radial > vector fields (with respect to the origin) th

[deal.II] Re: On using project_boundary_values_div_conforming for (RT-DG) in parallel case

2019-09-11 Thread Konrad
Hi Charlie, I can compile and run it. Could there be a problem with your deal.ii setup? I am using v9.1.1. Best, Konrad -- 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 messa

Re: [deal.II] Requirements for citing deal.II

2019-09-11 Thread 'Maxi Miller' via deal.II User Group
Ok, will do that, thanks! Am Mittwoch, 11. September 2019 16:52:51 UTC+2 schrieb Wolfgang Bangerth: > > On 9/11/19 3:27 AM, 'Maxi Miller' via deal.II User Group wrote: > > In one of my side projects I include the library deal.II, but use it > > only for reading in parameter files and printing ou

[deal.II] tutorial step-36 failed with PETSc error number 76 while calling a SLEPc function

2019-09-11 Thread Pragalv Karki
Dear deal.II community, I am new to deal.II and I had posted a question about this same issue yesterday. Perhaps I did not ask it correctly. I have been trying to simulate the wave equation on a particular kind of surface (I have attached a picture of that surface and the mesh

Re: [deal.II] Question about (parallel) linear solver

2019-09-11 Thread Konrad
Hi Wolfgang, > > Using approach 1.) I get timing results like that: > > (using MPI with Petsc, AMG preconditioning for *A* on one machine with > > four cores, no preconditioning for the Schur complement *S = > > -B*A^{-1}*B^T + C* from 1.) > > > > Running using PETSc. > > Number of active

Re: [deal.II] Question about (parallel) linear solver

2019-09-11 Thread Wolfgang Bangerth
On 9/11/19 7:04 AM, Konrad wrote: > > The system that needs to be solved formally reads > *(P)   A sigma + B^T u = 0 > * > *        B sigma + C u = f* > > *A* is sort of a mass matrix so it is positive (corresponds to > ** in *H(curl)*) and *C* is positive semi-definite > (corresponds to ** in

Re: [deal.II] Stabilization for Q1Q0 Stokes by Bochev [2006]

2019-09-11 Thread Richard Schussnig
Am Dienstag, 10. September 2019 13:26:56 UTC+2 schrieb Bruno Blais: > > I second Wolfgang comment on the fact that Q1Q1 is not difficult to > implement. You can also scale it to arbitrary Qn-Qn elements if you are > interested in higher order. > We have implemented such an approach in our code

Re: [deal.II] Stabilization for Q1Q0 Stokes by Bochev [2006]

2019-09-11 Thread Richard Schussnig
Am Dienstag, 10. September 2019 13:26:56 UTC+2 schrieb Bruno Blais: > > I second Wolfgang comment on the fact that Q1Q1 is not difficult to > implement. You can also scale it to arbitrary Qn-Qn elements if you are > interested in higher order. > We have implemented such an approach in our code ba

Re: [deal.II] Stabilization for Q1Q0 Stokes by Bochev [2006]

2019-09-11 Thread Richard Schussnig
Dear Wolfgang, Thank you very much for looking taking the time to answer my question, I really appreciate it - such things do not happen to often on the internet! ; ) Please find below my response per answer, I hope it is not too unstructured that way! Am Dienstag, 10. September 2019 06:03:27 U

[deal.II] Re: On using project_boundary_values_div_conforming for (RT-DG) in parallel case

2019-09-11 Thread Charlie Z
> > The minimal code. > -- 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

Re: [deal.II] Requirements for citing deal.II

2019-09-11 Thread Wolfgang Bangerth
On 9/11/19 3:27 AM, 'Maxi Miller' via deal.II User Group wrote: > In one of my side projects I include the library deal.II, but use it > only for reading in parameter files and printing out data in the HDF5 > format. I included it mostly because it is in use in most of my main > projects, and th

[deal.II] On using project_boundary_values_div_conforming for (RT-DG) in parallel case

2019-09-11 Thread Charlie Z
Hi there, I'm trying to make step-20 parallel through trilinos (by following step-40/55), but haven't succeeded. In setup_system() function, I set owned_partitioning, relevant_partitioning and constraints in the following way: dof_handler.distribute_dofs (fe); DoFRenumbering::component_wise(d

[deal.II] Re: On the usage of Utilities::MPI

2019-09-11 Thread Konrad
> > However, when I run the command "mpirun -np 5 ./MPI_TEST", I get the > following message, > > Hello world from process 0 of 1 > Hello world from process 0 of 1 > Hello world from process 0 of 1 > Hello world from process 0 of 1 > Hello world from process 0 of 1 > > It seems that 5 cpu cores a

[deal.II] Re: On the usage of Utilities::MPI

2019-09-11 Thread Bruno Turcksin
Yang, Are sure that you compiled deal.II with MPI support? Also can you post what you see when you do make VERBOSE=1 Best, Bruno On Wednesday, September 11, 2019 at 8:44:29 AM UTC-4, Yang Liu wrote: > > Dear Deal.II developers, > > Greetings! > > I am a new user of Deal.II. Currently I would

[deal.II] Question about (parallel) linear solver

2019-09-11 Thread Konrad
Dear deal.ii community, Apologies for so many questions lately... but here is another one. This one is a bit technical and aims at good solvers. In a simplified form you can formulate it this way: I implemented a PDE solver with a for the 3D vector Laplacian *curl(curl u) - grad(div u) = f*. I

[deal.II] On the usage of Utilities::MPI

2019-09-11 Thread Yang Liu
Dear Deal.II developers, Greetings! I am a new user of Deal.II. Currently I would like to test the hello world program on Utilities::MPI. A minimum working example is attached to this email. However, when I run the command "mpirun -np 5 ./MPI_TEST", I get the following message, Hello world fro

[deal.II] Requirements for citing deal.II

2019-09-11 Thread 'Maxi Miller' via deal.II User Group
In one of my side projects I include the library deal.II, but use it only for reading in parameter files and printing out data in the HDF5 format. I included it mostly because it is in use in most of my main projects, and thus it already was available. Furthermore, being able to reuse parameter