Re: [deal.II] Re: Modifying shape function data in MatrixFree

0}). But from what I learnt, double precision should be at least accurate up to ~10^{-15}. There are five orders of magnitude in between. I probably miss something here in between. Any idea? Regards, David Jean-Paul Pelteret schrieb am Mittwoch, 30. Dezember 2020 um 21:27:30 UTC+1: HI Da

Re: [deal.II] Div-conforming elements in 3D

ted for Nedelec elements. Best, Konrad On Saturday, December 12, 2020 at 8:07:19 PM UTC+1 Konrad Simon wrote: Hi Jean-Paul, On Thursday, December 10, 2020 at 11:39:08 PM UTC+1 Jean-Paul Pelteret wrote: HI Konrad, I have no familiarity

Re: [deal.II] Re: Modifying shape function data in MatrixFree

Hi David, So as I'm sure that you know, if you want to assemble the linear system for quasi-static non-linear (hyper)elasticity with a referential DoFHandler (using no Eulerian mapping to transform the shape functions to the spatial setting) AND without doing a similar transformation by hand

Re: [deal.II] Quasi Static Compressible Finit Strain for Heterogenous Hyperelastic material

he cook_membrane's program and cmake command used to build Trilinos in my system. Thanks! Animesh On Wednesday, November 25, 2020 at 1:08:57 PM UTC+5:30 Jean-Paul Pelteret wrote: Dear Animesh, Well, with such a small incremental displacement its probable that you're either wo

Re: [deal.II] Div-conforming elements in 3D

HI Konrad, I have no familiarity with the H-div elements, so I could be wrong with this suggestion... The Fe_Nedelec element suffered from a similar issue, where adjacent cells had to agree on the edge sense/directionality. This requirement could not be unconditionally guaranteed for that

Re: [deal.II] Condition number and Eigen Values of a Matrix

If it's of any help, our CG and GMRES solvers (and maybe some of the others?) can return the estimated condition number during / after a call to solve(). https://dealii.org/developer/doxygen/deal.II/classSolverCG.html#a52f1fd3d41765ec06d2d0add1a74a37c

Re: [deal.II] Quasi Static Compressible Finit Strain for Heterogenous Hyperelastic material

On Thursday, November 12, 2020 at 1:34:30 PM UTC+5:30 Jean-Paul Pelteret wrote: Hi Animesh, > In the documentation in code gallery <https://dealii.org/developer/doxygen/deal.II/code_gallery_Quasi_static_Finite_strain_Compressible_Elasticity.html>, it is written tha

Re: [deal.II] Use Gmsh to import external grid calculation elasticity problem

Hi Nick, Here's where the documentation for create_source_vector() is https://dealii.org/current/doxygen/deal.II/namespaceVectorTools.html#ac4e531baa71a5183079fd2083aa7d723 . What this

Re: [deal.II] Adaptive Time Increment Algorithm

s willing to share a portion of his code I'm going to try not to re-invent the wheel. James On Thursday, October 15, 2020 at 2:44:01 PM UTC-4 Jean-Paul Pelteret wrote: Dear James, I've got some snippets of code that I'll post below to show you how this can be done within the f

Re: [deal.II] Quasi Static Compressible Finit Strain for Heterogenous Hyperelastic material

Jean, Sorry for late response. Thank you for your detailed answer. Using your suggestions I was able to assign the material ids to different area of my domain of interest. Thanks again! Animesh On Tuesday, October 27, 2020 at 3:05:39 AM UTC+5:30 Jean-P

Re: [deal.II] Quasi Static Compressible Finit Strain for Heterogenous Hyperelastic material

Animesh On Tuesday, October 27, 2020 at 3:05:39 AM UTC+5:30 Jean-Paul Pelteret wrote: Dear Animesh, Let me preface my suggestion by saying that you can introduce some heterogeneity in more than one way. That is, you're not bound to using a material

Re: [deal.II] Re: Error with boost after installing dealii with spack

Hi Christian, It sounds like what might be happening is that your project is preferrentially picking up a system boost rather than the one that you've compiled deal.II with. Marc's suggestion of building the development version with that patch (#11024

Re: [deal.II] Quasi Static Compressible Finit Strain for Heterogenous Hyperelastic material

Dear Animesh, Let me preface my suggestion by saying that you can introduce some heterogeneity in more than one way. That is, you're not bound to using a material ID (which is what I'm going to suggest), but you could use some geometric arguments to determine which class or material

Re: [deal.II] Adaptive Time Increment Algorithm

Dear James, I've got some snippets of code that I'll post below to show you how this can be done within the framework of step-44. The basic concept that I implemented is as follows: - The nonlinear solver is called by some function who's job it is to solve for a single timestep.

Re: [deal.II] Application of the non-homogeneous boundary conditions in Finite Elasticity (Compressible case)

Hi Animesh, Here are some answers to your questions: > Essentially, I will have to apply the homogeneous boundary condition at the > zeroth and first Newton's step, and Non-homogeneous BC only at the zeroth > step, as mentioned in the documentation >

Re: [deal.II] Fluid-Structure interaction

Hi Ramiro, I would also suggest taking a look at the publication list to understand a bit better what can be achieved using deal.II (and, hopefully, how to do it). Here, for example, is some work done in vascular biomechanics, and the publication list has several other references to FSI.

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

Hi Behrooz, I must admit that I’ve not been able to follow why you want to store Points when you’re interested in global vertex indices, but I’ll try to give a little help where I can. In case its of any help to you, I believe that you should be able to get the coordinate of the vertex with a

Re: [deal.II] Node duplication in deal.ii

Hi Behrooz, To complement Simon’s suggestion, I can think of two more ways to achieve what you’re wanting but will affect the mesh itself. Whether or the second suggestion is a reasonable approach for you depends on the complexity of the geometry that you’re creating. 1. Use an external mesh

Re: [deal.II] install problems with clang@12.0.0 and Xcode12.0 from spack

Hi Alexander,I am, unfortunately, marginally acquainted with the problem that you’re facing. I’m not 100% sure, but I believe that this is an upstream problem that has to do with how these packages are configured for Apple systems. I too noticed that within a couple of update cycles SuiteSparse,

Re: [deal.II] Eigenvalues and Eigenvectors of BlockSparse Matrix

> AR > On Sunday, September 27, 2020 at 8:36:20 PM UTC+5:30 Animesh Rastogi IIT > Gandhinagar wrote: > Hi Jean, > > Thanks a lot for your detailed answer. Using this idea, I was able to solve > the standard eigenvalue problem. > > Thanks again! > > Animesh

Re: [deal.II] Function apply_boundary_values()

To add to what Wolfgang already said, there’s also the constrained_linear_operator() that may be of use in this scenario. Best, Jean-Paul > On 25 Sep 2020, at 23:29, Wolfgang Bangerth

Re: [deal.II] Issue with unstructured hex mesh

ide-meshing-snappyhexmesh.html> I’ve never used it though, and I have no idea how you would get the mesh in and out of OpenFOAM, and whether one can prevent the creation of meshes with hanging nodes. Nevertheless, maybe its worth looking into. Best, Jean-Paul > On 25 Sep 2020, at 21:58,

Re: [deal.II] Eigenvalues and Eigenvectors of BlockSparse Matrix

ck(u_dof,u_dof)); > const int eigensolver_its = static_cast( > tangent_matrix.block(u_dof, u_dof).m() > * parameters.max_iterations_lin); > SolverControl solver_control(eigensolver_its, 1e-9); > Ar

Re: [deal.II] Issue with unstructured hex mesh

Since you have access to Cubit/Trellis, “that patch” is probably not necessary. Cubit uses the Mesquite library for mesh quality improvement, so if you could try to increase the mesh quality that way. Although I’d have to agree with Wolfgang that it might not help too much in this case, but

Re: [deal.II] Eigenvalues and Eigenvectors of BlockSparse Matrix

Hi Animesh, Although in that code-gallery example the system is assembled into a BlockSparseMatrix, the system actually has only one block. So you can retrieve the underlying SparseMatrix (and Vector for the RHS) via tangent_matrix.block(u_dof, u_dof); system_rhs.block(u_dof); and use them

Re: [deal.II] Multi-physics implementation

Hi Ray, I’m happy to see that Simon has already given a great explanation as to how to solve your problem. It’s probably more on point than the code that mentioned that I would direct you towards. So it really was worthwhile asking the question again here :-) Best, Jean-Paul > On 22 Sep

Re: [deal.II] Separated domains in load balanding

Ah, I was just busy writing what Wolfgang said. You can find more info on custom partitions for parallel::shared::Triangulation here: https://dealii.org/current/doxygen/deal.II/classparallel_1_1shared_1_1Triangulation.html#a686a3453dfec098eb64d1510aa1716e1

Re: [deal.II] Can not access the files for Quasi-Static Finite-Strain Compressible Elasticity program

Timo, I think that what Animesh is referring to relates to this issue: https://github.com/dealii/dealii/issues/10295 None of the hyperlinks inside the code gallery documentation are working. J-P > On 10 Sep 2020, at 20:25, Timo Heister wrote: >

Re: [deal.II] extract_dofs() extremely slow in for single-thread run?

disables the reinit()-call, and I have > to call cell->is_locally_owned(), instead of cell->locally_owned(). > Thanks! > > Jean-Paul Pelteret schrieb am Dienstag, 1. September 2020 um 20:28:18 UTC+2: > Hi, > > 30-40% of the runtime for this operation seems a bit excessiv

Re: [deal.II] Re: Thermal Modelling of Additive Manufacturing Process

Dear Pushkar, Indeed, as Bruno said, it’s hard to understand what information you’re looking for. This can certainly be done with deal.II, as it’s a topic of research that my former colleague has been working on for some time:

Re: [deal.II] extract_dofs() extremely slow in for single-thread run?

Hi, 30-40% of the runtime for this operation seems a bit excessive to me. Maybe you should try using the FEValuesExtractors in conjunction with FEValues in order to get the solution at the quadrature points. If you want the solution at the support points, then you can just create a quadrature

Re: [deal.II] How to extract nodes/vertices at the interface of two domains?

Dear Behrooz, I guess that the first thing that you need to be able to do is to traverse the mesh until you can get an iterator and index to the cell face that you’re interested in. If the manifold already has some sort of ID associated with it, then you can just do something like

Re: [deal.II] deal.ii Installation : expand_instantiations Error

ssors(perhaps 4) and if even that doesn't work, then I will try it > in serial. > > Best, > Aaditya > > On Friday, August 21, 2020 at 1:07:54 PM UTC-4 Jean-Paul Pelteret wrote: > Hi Aaditaya, > > I don’t think that I’ve ever seen anyone have this problem before. We

Re: [deal.II] deal.ii Installation : expand_instantiations Error

Hi Aaditaya, I don’t think that I’ve ever seen anyone have this problem before. We build the expand_instantiations binary, and its supposed to be used when expanding files in which we define various combinations of template parameters to instantiate templates within the library. If you’re

Re: [deal.II] Nested reverse and forward-mode Sacado memory leak

Hi Doug,Yes, sadly I have experienced the same but wasn’t able to allot the time to figure out if it was indeed a bug in Sacado (and maybe something that’s already been fixed upstream in a newer version of Trilinos than I was using) or something that we did. I’m really glad that you managed to

Re: [deal.II] ABOUT OUTPUT IN PARALLEL

Hi, The problem is most likely that you’ve got on (static) variable “times_and_names” that is referencing two completely different of output files. This means that for every tilmestep, you’ve now got two entries in this file and Paraview probably doesn’t know what to do with it. You can easily

Re: [deal.II] omp_simd error in step-67 when installing deal.ii 9.2.0 with spack

Hi Ester, There’s a quick workaround for this. You should just edit this line of the deal.II spack package: https://github.com/spack/spack/blob/develop/var/spack/repos/builtin/packages/dealii/package.py#L252

Re: [deal.II] About metallic ductile damage

Dear Itachi, To supplement what Bruno has already mentioned, deal.II provides some of the tools what you might find useful to implement constitutive laws or specific governing equations. To get some feeling as to what is possible, I would suggest that you take a look through the list of

Re: [deal.II] Deal.ii installation

Dear Prasad, I’m glad to hear that you solved your problem. I must say, though, that the problem and the solution only seems to me to be tangentially related. Although neither you, I, nor probably anyone else could have guessed from the cryptic and nondescript error message, you were probably

Re: [deal.II] Particle Class Implementation

Dear Victoria, You’re on the right track, but it looks like you got the syntax a bit wrong. To create a particle you need to write something like Particles::Particle particle; You’ll note that the specific constructor that is called is inferred from the arguments given to the instance of the

Re: [deal.II] preconditioner for three block matrix

Dear Zeb, To expand on what Wolfgang said in his previous email, we have a page that provides some hints on how to formulate a question that you can then engage the community with. Please give it a read, and have another go at asking

Re: [deal.II] Resolving with different RHS

Hi Andreas, I concur with Simon. Additionally, if you’re using the Trilinos linear algebra classes, then the TrilinosWrappers::SolverDirect solver will do the same https://dealii.org/developer/doxygen/deal.II/classTrilinosWrappers_1_1SolverDirect.html#a6c0fecc921bc93f2c005e265cfbdd882

Re: [deal.II] problem measure() for quads

Hi Andreas, Thanks for being willing to make a pull request of this. The basic instructions for how to do so are listed in our FAQ, https://github.com/dealii/dealii/wiki/Frequently-Asked-Questions#i-found-a-typo-or-a-bug-and-fixed-it-on-my-machine-how-do-i-get-it-included-in-dealii

Re: [deal.II] vertex and dofcells

Dear Alberto, Quick answer: There’s a entry in the FAQ that may answer your question: https://github.com/dealii/dealii/wiki/Frequently-Asked-Questions#how-do-i-get-the-degree-of-freedom-indices-at-vertices

Re: [deal.II] How to download deal ii 8.2.1 virtual machine image

for sending the link for previous versions. Should I download all the > files or the .tar.gz file is enough to install an 8.2.1 version. > > Thanks & Regards > Ajith > > > On Fri, May 1, 2020 at 2:21 PM Jean-Paul Pelteret <mailto:jppelte...@gmail.com>> wrote:

Re: [deal.II] How to download deal ii 8.2.1 virtual machine image

Dear Ajith, To supplement what has been said before, older versions of deal.II can be found on the GitHub release page and our mirror: https://github.com/dealii/dealii/releases https://dealii.43-1.org/downloads/ It

Re: [deal.II] My machine did't manage to run Deal.II

Dear Wilmar, Indeed, Francis has the right suggestion, although I’m a bit worried that its still not going to help you (unless writing into swap is something that you’d consider). A similar question had been asked a couple of years ago, and the answer that they got at the time was that > 1GB

Re: [deal.II] discrepancy in step-44 material model between implementation and formula.

will take the discussion beyond the scope of > the question I asked earlier. I will ask that later. Thank you once again. > > > Navneet R > > > > On Wed, Mar 18, 2020 at 3:56 AM Jean-Paul Pelteret <mailto:jppelte...@gmail.com>> wrote: > Dear Navneet, &

Re: [deal.II] Periodicity with simultaneous displacement controlled BCs

Dear David, Indeed, Wolfgang is correct. I gather that you’re trying to produce the sort of displacement field that you see in the attached image. This reference > C. Miehe. Computational micro-to-macro transitions for discretized > micro-structures of hetero- geneous materials at finite

Re: [deal.II] Step 44: imported 2d mesh fails to converge to solution

> It might be that one of the constrained boundaries is not traction free... Correction: This should have read "It might be that one of the constrained boundaries is now traction free..." -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see

Re: [deal.II] Step 44: imported 2d mesh fails to converge to solution

Dear James, To me, the most likely point of failure here relates to the application of boundary conditions. It might be that one of the constrained boundaries is not traction free, and so the problem is not free of rigid body modes. Can you please confirm that you apply the same numbering

Re: [deal.II] Adding the manifold ID to a DataOutFaces

value. Maybe someone else on the mailing list knows more about this class, and could provide more hints? > On 10 Apr 2020, at 00:45, Jean-Paul Pelteret wrote: > > Hi Bruno, > > I haven’t used this class before, and admittedly haven’t gone digging too > deep for examples of w

Re: [deal.II] Adding the manifold ID to a DataOutFaces

Hi Bruno, I haven’t used this class before, and admittedly haven’t gone digging too deep for examples of where its used. But here are a few examples (one tutorial and a couple from the test suite) that may be able to guide you in the right direction. At first glance, its seems to work

Re: [deal.II] Difficulty setting manifold with Opencascade using gmsh mesh + salome STEP (or IGES) - Step-54

he functionnality I will add it to the >> FAQ. Thanks for adding this section. >> Best >> Bruno >> >> >> On Friday, 3 April 2020 16:11:05 UTC-4, Jean-Paul Pelteret wrote: >> Dear Bruno, >> >> I guess that it would be best to mention this in th

Re: [deal.II] Local tutorial installation missing some examples

Hi Xuefeng, The list of tutorials that are contained in the latest release (currently 9.1.1, which you’re using) are shown on this page: https://www.dealii.org/current/doxygen/deal.II/Tutorial.html The tutorials that your find

Re: [deal.II] ABOUT ASSEMBLING RIGHT HAND SIDE

Hi Xiang, You’ve also not mentioned whether you’re comparing vectors before or after you’ve distributed the constraints. Perhaps you get the same result if you call constraints.distribute() on both of your vectors? Best, Jean-Paul > On 07 Apr 2020, at 22:22, Wolfgang Bangerth wrote: > > On

Re: [deal.II] Difficulty setting manifold with Opencascade using gmsh mesh + salome STEP (or IGES) - Step-54

ily available. Clearly I would not have found that out > without Luca's help. > Best > Bruno > > > On Friday, 3 April 2020 09:43:33 UTC-4, Jean-Paul Pelteret wrote: > Dear Bruno, > > I’m glad that you managed to solve this issue with Luca’s aid. I think that >

Re: [deal.II] Difficulty setting manifold with Opencascade using gmsh mesh + salome STEP (or IGES) - Step-54

Dear Bruno, I’m glad that you managed to solve this issue with Luca’s aid. I think that the information and workflow that Luca gave was really valuable. Would you care to add this to the Wiki page? Best, Jean-Paul > On 03 Apr 2020, at 06:03, Bruno Blais wrote: > > Ciao Luca, > It works

Re: [deal.II] discrepancy in step-44 material model between implementation and formula.

case of > three variable formulation and are independent hence , as . > > Regards, > Navneet R > > > > > On Fri, Mar 6, 2020 at 2:47 AM Jean-Paul Pelteret <mailto:jppelte...@gmail.com>> wrote: > Dear Naveet, > > I think that I see where thi

Re: [deal.II] Installation with spack - no access to adol-c

Hi Konrad, Yes, we need some version of ADOL-C greater than 2.6.3, as we upstreamed some required patches to the project. I installed deal.II using Spack a few weeks ago, and didn’t run into any difficulty... but I see that the Gitlab project is indeed locked down, which is a bit surprising.

Re: [deal.II] Macro periodic boundary condition for a elasticity equation

Hi Lukas, If I’ve interpreted your emails correctly, you’ve mentioned two distinct things here, namely (1) the imposition of periodic displacement fluctuation field (thereby, in essence, prescribing a macroscopic average deformation gradient, with stress unknown) or (2) the imposition of a

Re: [deal.II] Instantiation problem for Utilities::MPI::sum (const ArrayView< const T > , const MPI_Comm _communicator, const ArrayView< T > )

Dear Feimi, On first glance, it looks like the object you send in as the first argument to this function has the wrong type. The function expects an ArrayView< const int >, but you are passing it an ArrayView< int >. That’s the first thing that I spotted, but hopefully its just this that’s an

Re: [deal.II] Different shape representations with manifolds on the same triangulation

er" > triangulation and > apply different mappings (alpha_m) to the cell iterator as we loop. > > Maybe this is not currently possible? > > Thanks in advance, > > Juan Carlos Araújo, > > > On Wednesday, 15 January 2020 22:19:46 UTC+1, Jean-Paul Pelter

Re: [deal.II] Using solution on from one FE problem as a boundary condition for another

Hey Ernesto, I’ve used the approach outlined by Wolfgang and Daniel in my own work (multi physics problems, single DoFHandler with BC’s being transferred between a primary and secondary problem solved on subsets of the domain). Overall, it works really well so I’d back up their suggestion and

Re: [deal.II] Facilities for computing the "goodness" of a mesh (i.e. something like a mesh-quality indicator)

Dear Krishna, There is also this long-standing pull request https://github.com/dealii/dealii/pull/8196 that hooks into the Verdict library to do mesh quality assessment. One day, hopefully in the not too distant future, I’ll get back to integrating

Re: [deal.II] SOS - UNABLE to complete coding assignment 1 of Coursera - problems with UMFPACK - sparsity issues

Hi Jordi, Indeed, it is as Wolfgang says. Unfortunately none of the developers know this course, or how to what extent deal.II is used in it. It was good timing that I was able to help you out during the vacation period (and you’re most welcome for that, if it did assist you), but time is

Re: [deal.II] SOS - UNABLE to complete coding assignment 1 of Coursera - problems with UMFPACK - sparsity issues

f > the Jacobian in the basis gradients but, for 1d and 2d this has not to be > affecting due to the same spacing between nodes. Can you somehow "decript" > your answer a little bit, I am really struggling with this assignment XD > > Kind regards and thanks again! > &g

Re: [deal.II] SOS - UNABLE to complete coding assignment 1 of Coursera - problems with UMFPACK - sparsity issues

Dear Jordi, Firstly, welcome to the forum. Well, as this is a course assignment I think that its fair to assume that some course mentor should be available to guide you toward the solution to the problem. You say that you’ve had no luck on the course site forum — that’s a bit poor. Anyway, I

Re: [deal.II] Solve for eigenvalues for small matrix in serial

Hi Pai, If you’ve got LAPACK installed then the LAPACKFullMatrix class is able to do the job for you. It has several methods implemented to compute eigenvalues, and you’d then use the eigenvalue()

Re: [deal.II] Instantiation of Vector< Sacado::Fad::DFad >

H Doug and Wolfgang, At the time when we introduced the AD framework I took a stab at quickly adding support for to our Vector class. To summarise, I hit the same issue and couldn’t find an elegant solution past to get past it (so hopefully this is the only outstanding problem to add support

Re: [deal.II] Number of Lines in the Deal.ii Library

-- > Language files blankcomment > code > --- > C++ 4842 114261 103709 > 399925 > C/C++ Header84 3254 3014

Re: [deal.II] Number of Lines in the Deal.ii Library

Hi Jason, Sorry, replying to your message has been on my TODO list for quite some time. Thanks for sharing! Those numbers are certainly thought provoking. It would appear that we (being the collective community) have been putting in a few minutes here and there to write this FEM library :-D

Re: [deal.II] Different solution with different software

e too hard on yourself. We’ve all made (and will continue to make) these sorts of mistakes from time to time :-) Best, Jean-Paul > On Thursday, July 25, 2019 at 7:28:13 PM UTC+2, Jean-Paul Pelteret wrote: > Dear Felix, > > I don’t claim to have read through your comprehensive document

Re: [deal.II] Different solution with different software

Dear Felix, I don’t claim to have read through your comprehensive documentation of your issue in any detail, but one observation of your code is that you have a couple of suspect terms in your local matrix assembly. In particular, the terms of this form > 1/2/phiscale*(u_j_y-v_j_x) do not

Re: [deal.II] How can I get the global coordinate of Integrated Point

Dear Jacob, You should have a look at step-8 , for example. In assemble_system() some local material properties are computed by evaluating some functions at all cell quadrature points: FEValues fe_values(fe, quadrature_formula,

Re: [deal.II] configuration with LAPACK/BLAS fails

Hi Ahmad, Although I don’t have an answer for you, I just want to say that what you’re going here looks somewhat reasonable to me. Its unclear to me why this test would fail: /usr/bin/cc -DLAPACK_SYMBOL_CHECK -pthread-rdynamic CMakeFiles/cmTC_3cdf4.dir/src.c.o -o cmTC_3cdf4 -lgfortran

Re: [deal.II] libraries for line search

Dear Dhananjay, We have a family of line search algorithms implemented in the LineMinimization namespace (https://github.com/dealii/dealii/blob/master/include/deal.II/optimization/line_minimization.h

Re: [deal.II] how to get the inverse of a sparse matrix?

ons to global block. I feel that this approach may not > work as I am dealing with continuous fields. > > I would be glad to know about your thoughts regarding this. > > Thanks, > Krishanu > > On Thursday, 6 June 2019 14:39:16 UTC-5, Jean-Paul Pelteret wrote: > Hi K

Re: [deal.II] how to get the inverse of a sparse matrix?

Hi Krishna, One of my answers in this thread provides the solution you’re looking for, namely the suggestion of the LinearOperator class and its associated helper functions. Best, Jean-Paul > On 06 Jun 2019, at 19:41, Krishanu Sen wrote: > > Hello, > > I have a similar problem where

Re: [deal.II] Interface boundary attribute of partitioned meshes

Dear Tom, You’ve come to the right place to ask your question but, occasionally, it takes some time for the right person to come along to be able to give a meaningful answer to it. I, for example, have no familiarity with Additive Schwarz methods, so I have no insights to share with you on

Re: [deal.II] Problem Solving BlockSparsematrix system blockwise

Dear Gabriel, It looks to me like you forgot to dereference the shared pointer for the preconditioner. So, solver_gmres.solve( system_matrix.block(0,0), distributed_solution.block(0), rhs.block(0), *Mp_preconditioner); // <— here You probably have the same problem for the other call.

Re: [deal.II] Adding TrilinosWrapper::BlockSparseMatrix

Dear Gabriel, There is a function to do this operation https://dealii.org/9.0.0/doxygen/deal.II/classBlockMatrixBase.html#a9ee9987ece33a052d2b7a44034a8a01a So you should write this as

Re: [deal.II] Controllins fill-in, abs. tol and rel. tol or smoother and coarse solver of TrilinosWrappers::PreconditionAMG

_ml.set("smoother: ifpack absolute threshold",ilu_atol); > parameter_ml.set("coarse: ifpack level-of-fill",ilu_fill); > parameter_ml.set("coarse: ifpack absolute threshold",ilu_atol); > preconditioner.initialize(system_matrix,parameter_ml); > > >

Re: [deal.II] Controllins fill-in, abs. tol and rel. tol or smoother and coarse solver of TrilinosWrappers::PreconditionAMG

e. I have found the MULUE documentation, but it seems to be more tailored > around their XML interface? > > > Thank you so much, > Bruno > > > > On Tuesday, 12 March 2019 05:06:03 UTC-4, Jean-Paul Pelteret wrote: > Dear Bruno, > > I believe that my prop

Re: [deal.II] Installation problem

Dear Wenguo, It looks like you’re linking deal.II, built as a shared library, against LAPACK, built as a static library. As the error message suggests, you’ll need to build LAPACK with the -fPIC flag in order to do this. As suggested here

Re: [deal.II] Accuracy of Dirichlet condition for p in step-20

Dear Jane, > const Functions::ZeroFunction no_tang_bcs; This should be const Functions::ZeroFunction no_tang_bcs(n_components); Best, Jean-Paul -- 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 ---

Re: [deal.II] Finding the y-coordinate of the boundary during the assembly

There’s a function to compute the bounding box in GridTools: https://www.dealii.org/current/doxygen/deal.II/namespaceGridTools.html#ae1ec55abefa31cf001fd29d8d4d993f1 > On 16 Mar 2019, at

Re: [deal.II] Interfacing deal.ii for solution, pre and post-processing

Dear Muhammad, You can look at the documentation to the GridIn namespace for some ideas as to which software supports the generation of mesh files that deal.II can read, and the and GridOut

Re: [deal.II] Finding the y-coordinate of the boundary during the assembly

Dear Jane, There is the BoundingBox class that could tell you the overall extents of the domain. So if your geometry is rectangular and not rotated then you could work out the difference between the y-coordinate of the top

Re: [deal.II] what is the different between VectorTools::point_value and fe_values.get_function_values()?

Dear Pham, When you use the FEValuesExtractors in conjunction with the FEValues class using the operator[] like you’ve done here > fe_values[vector_re] then what is returned is an FEValuesViews type object. Since your vector_re extractor is of type FEValuesExtractors::vector, the returned

Re: [deal.II] Controllins fill-in, abs. tol and rel. tol or smoother and coarse solver of TrilinosWrappers::PreconditionAMG

Dear Bruno, I believe that my proposed changes are working now, so you could always pull that branch from my repository and use it immediately. That might be the most convenient way for you to move forward. Thanks for being willing to share your findings wrt. AMG parameters. I think that

Re: [deal.II] Controllins fill-in, abs. tol and rel. tol or smoother and coarse solver of TrilinosWrappers::PreconditionAMG

Dear Bruno, Indeed, it appears that one isn't able to change the settings passed to the AMG preconditioner other than those set by the AdditionalData object. Your only current alternative is to create the Teuchos::ParameterList yourself, which you can send in to the call to initialize()

Re: [deal.II] Error in the code

Dear Muhammad, Can you confirm something: does the directory /media/muhammad/Data/Deal/dealii-8.3.0/ point to your source directory for deal.II, or the location of the installed library? Also, what does your CMakeLists.txt look like? If you are invoking autopilot like this

Re: [deal.II] Regarding billet proof jacket

Dear Nitin, Your question is exceptionally broad, so its not possible to give you any meaningful answer or feedback on it. Perhaps you should take a looks at our guidelines for asking questions on the forum which are located here https://groups.google.com/forum/#!topic/dealii/GRZMUTLIm2I

Re: [deal.II] Shared memory parallelism of Trilinos

Dear Bruno, > Thank you for the rapid and very detailed answer. This makes this community > so great. Bruno T. has given you far more insight into the inner workings of the Trilinos linear algebra packages than I’d have been able to (thanks :-) ), so I’m really glad that we were able to help

Re: [deal.II] Error installing complex_bessel library

Dear Muhammad, This is not really a deal.II related question, but I’ll try to help you anyway. CMake is telling you that you don’t have a Fortran compiler installed, so you should install one. You should be able to do this with “sudo apt install gfortran “. I hope that this help you. Best,

Re: [deal.II] Shared memory parallelism of Trilinos

Dear Bruno, > The suggested compilation options for Trilinos do not suggest the use of > OpenMP and the flag is not enabled by default. Can you please indicate where you read this? I’m looking at our documentation for interfacing to external libraries, i.e.

Re: [deal.II] Compile error with function compute_bounding_box on GCC 8.5

Hi Reza, Thanks for letting us know about this problem. Yes, its best to report these types of issues directly on the GitHub repository :-) I’ve opened an issue based on the details that you given. You can track it here: https://github.com/dealii/dealii/issues/7771

Re: [deal.II] Can't find affineconstraints related header file

Dear Chen, In deal.II version 9.0, affine constraints are applied using the ConstraintMatrix class (deal.II/lac/constraint_matrix.h). This class and header have been deprecated in the developer version, so if you’re using the developer version the ConstraintMatrix is typedef’d to

Re: [deal.II] Tutorial step-35 fails with my inp

Dear Eugeny, Have you tried to increase the “max_iterations” in the parameter file? Basically, with the way the tutorial program is defined the solver is currently limited to use the number of iterations defined by the variable “vel_max_it”. This seems to be prescribed using a parameter file,

<    1   2   3   4   5   6   >