Re: [deal.II] Making a pair of extracted data and cast into a set

2020-10-13 Thread Behrooz Karami
Dear Jean-Paul, Thanks very much for your explanations. I managed to achieve my objective. You were right, one of the crucial points I had to take into account was the compatibility of data structures. Defining a new operator was another to the point advice, as some of the errors were laid on

[deal.II] Adaptive Time Increment Algorithm

2020-10-13 Thread James Gorman
deal.ii users, I'm learning deal.ii with the goal of using it to non-linear solid mechanics problems. In many of these problems (such as Step-44), a time step is used to move a calculation from the initial condition to the final loaded state. Step-44 uses a uniform time-step, but there are

Re: [deal.II] Re: APPLY BOUNDARY CONDTION ON INTERNAL FACE

2020-10-13 Thread 孙翔
Hi, Daniel, I see. Thank you very much. So, I think our method may be available as long as it can separate the face. Best, Xiang On Tuesday, 13 October 2020 at 14:38:29 UTC-7 d.arnd...@gmail.com wrote: > Xiang, > > deal.II doesn't have functionality to create separate cells sharing a face

Re: [deal.II] Re: APPLY BOUNDARY CONDTION ON INTERNAL FACE

2020-10-13 Thread 孙翔
Hi, Bruno, Thanks a lot. We will have a try based on the example you gave us. Best, Xiang On Tuesday, 13 October 2020 at 10:49:08 UTC-7 bruno.t...@gmail.com wrote: > Xiang, > > Here is an example on how to do it: > >

Re: [deal.II] Re: APPLY BOUNDARY CONDTION ON INTERNAL FACE

2020-10-13 Thread Daniel Arndt
Xiang, deal.II doesn't have functionality to create separate cells sharing a face but only to join them. Hence, whatever input you provide to Triangulation::create_triangulation() must have these faces already. GridIn should work for that if the faces are disjoint in the input file. Best, Daniel

Re: [deal.II] FullMatrix to SparseMatrix

2020-10-13 Thread Wolfgang Bangerth
On 10/13/20 10:28 AM, Nikki Holtzer wrote: The issue that I have when  all matrices are FullMatrix matrices is I can't seem to get the .add function to result in anything other than error. If I keep ONLY both nonlinear_part and linear_part as FullMatrix I receive: error: no matching function

Re: [deal.II] Re: APPLY BOUNDARY CONDTION ON INTERNAL FACE

2020-10-13 Thread Bruno Turcksin
Xiang, Here is an example on how to do it: https://github.com/geodynamics/aspect/blob/master/cookbooks/prescribed_velocity/prescribed_velocity.cc#L297-L303 As you can see, you need to know which dof you need to constrain Best, Bruno Le mar. 13 oct. 2020 à 11:02, 孙翔 a écrit : > > > Hi, Bruno,

[deal.II] Re: building with PETSc and METIS

2020-10-13 Thread Aaditya Lakshmanan
Hi, Although you have set -DDEALII_WITH_MPI=ON it seems from the following lines in your output : -- Found LAPACK -- Use other than Intel MKL implementation of BLAS/LAPACK (consult CMakeFiles/CMakeError.log for further information). -- DEAL_II_WITH_LAPACK successfully set up with external

Re: [deal.II] FullMatrix to SparseMatrix

2020-10-13 Thread Nikki Holtzer
The issue that I have when all matrices are FullMatrix matrices is I can't seem to get the .add function to result in anything other than error. If I keep ONLY both nonlinear_part and linear_part as FullMatrix I receive: error: no matching function for call to

[deal.II] building with PETSc and METIS

2020-10-13 Thread Christos Karliampas
Dear all, i try to build dealii with PETSc and METIS using command, cmake -DDEAL_II_WITH_PETSC=ON -DPETSC_DIR=/home/christos/petsc -DPETSC_ARCH=x86_64 -DDEAL_II_WITH_METIS=ON -DMETIS_DIR=/home/christos/metis-5.1.0 -DDEALII_WITH_MPI=ON -DCMAKE_INSTALL_PREFIX=/usr/local .. having the

Re: [deal.II] Re: APPLY BOUNDARY CONDTION ON INTERNAL FACE

2020-10-13 Thread 孙翔
Hi, Daniel, Thank you very much. So, if I want to disjoint the faces, how should I do using dealII or I need to disjoint them manually? If I manually disjoint them, should I duplicate the joint nodes and assign them with the new number, then assemble the adjacent element with the duplicated

Re: [deal.II] Re: APPLY BOUNDARY CONDTION ON INTERNAL FACE

2020-10-13 Thread 孙翔
Hi, Daniel, Thank you very much. So, if I want to disjoint the faces, how should I do using dealII or I need to disjoint them manually? If I manually disjoint them, should I duplicate the joint nodes and assign them with the new number, then assemble the adjacent element with the duplicated

[deal.II] Re: APPLY BOUNDARY CONDTION ON INTERNAL FACE

2020-10-13 Thread 孙翔
Hi, Bruno, Thank you very much. Does the AffineConstraint force the value at the specified DOF? I need to how the node numbers map to the DOF, right? Best, Xiang On Monday, 12 October 2020 at 06:21:03 UTC-7 bruno.t...@gmail.com wrote: > Xiang, > > You cannot apply a boundary condition on an