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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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

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

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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] 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] Stabilized FEM implementation with bubble function

Hi Lixing, Another tutorial that might be of interest to you to look at is step-44 . In that tutorial, two discontinuous fields are condensed out in one of two ways (there's a switch to choose which method is applied). The first

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] 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] 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] 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] 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] 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] 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] 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] ADOL-C not found through SPACK installation of deal.ii 9.2

Dear Behrooz, The problem here is that your Trilinos installation includes the SEACASChaco package. We'd discovered that there's a common global symbol that both ADOL-C and Trilinos-SEACASChaco emit, which causes some very strange issues when using ADOL-C. Naturally, we had to make sure that

Re: [deal.II] Re: issue in installing dealii@9.2 on ubuntu through spack

Hi Alberto, I concur with Bruno. ADOL-C changed the location of their project and repository, so Spack's master branch must be directing you towards the old site. Clearly it doesn't have an updated package for deal.II either, because its trying to deduce the location of the deal.II tarball on

Re: [deal.II] ADOL-C not found through SPACK installation of deal.ii 9.2

find the 'spack-build.out' file as attached. This is what Spack mentioned as a build log for details. Hope to be the correct one. Thank you very much, Behrooz On Wed, Jan 13, 2021 at 8:14 PM Jean-Paul Pelteret mailto:jppelte...@gmail.com>> wrote: Hi Behrooz, Great, I'

Re: [deal.II] ADOL-C not found through SPACK installation of deal.ii 9.2

any way around? Thank you very much, Behrooz On Friday, January 8, 2021 at 7:22:45 PM UTC+1 Jean-Paul Pelteret wrote: Dear Behrooz, The problem here is that your Trilinos installation includes the SEACASChaco package. We'd discovered that there's a common global symbol tha

Re: [deal.II] ADOL-C not found through SPACK installation of deal.ii 9.2

Hi Behrooz, When you install deal.II using Spack, it is completely isolated from the system environment. That way, you can install multiple versions of the same package without them conflicting with one another. Before you can use deal.II though spack, you populate all of the environment

Re: [deal.II] Re: issue in installing dealii@9.2 on ubuntu through spack

<http://m4lab.unibs.it/faculty.html> On Thu, Jan 14, 2021 at 7:57 PM Jean-Paul Pelteret mailto:jppelte...@gmail.com>> wrote: Hi Alberto, I concur with Bruno. ADOL-C changed the location of their project and repository, so Spack's master branch must be directing you to

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] Transferring data between two triangulations (dof handlers)

Hi Yu, Have you tried the VectorTools::interpolate_to_different_mesh() functions? It sounds to me like they perform the operation that you're trying to do. They require that the two

Re: [deal.II] Linear operators - TrilinosWrappers

op_S_inv = inverse_operator(opa_S, solver_S, preconditioner_S > ); > > I am not sure what is happening because all the operators are declared as > "LA::MPI::Vector". > If you have any suggestion, would be greatly appreciated. > > Thank you so much, > >

Re: [deal.II] Linear operators - TrilinosWrappers

ver_S, preconditioner_S ); > //(Case I could invert) > > (Please excuse me for the extension of the code, but the solver function > where I have the issue is too short. In the case you consider it necessary, I > could code a shorter one) > > Thank you so much, > >

Re: [deal.II] Deal.II Make error || Make release

Hi David and Sumanth, There is a new Mac package for the soon to be announced 9.3 release that is built for BigSur. Can you try that, and see if it works for you? If you have any troubles then we’d appreciate the feedback.

Re: [deal.II] Question about deal.II flexibility

Dear Abdelrahman, To add to some further details to what Konrad has already provided: - the step-44 tutorial implements one of many incompressible elasticity formulations. There is also the Quasi-Static Finite-Strain

Re: [deal.II] Questions on Step-18

velocity gradient is the vorticity > which is twice the angular velocity - hence I think you need a 1/2. Happy to > be corrected on this. > > Best, > Andrew > > > On 24 Jun 2021, at 21:03, Jean-Paul Pelteret <mailto:jppelte...@gmail.com>> wrote:

Re: [deal.II] Questions on Step-18

Hi Michael, I cannot comment on the first question, but might be able to assist a bit with the second. But may I first ask, what precisely are you trying to achieve with this extension? As interesting as it is, in the past I had found step-18 to deviate too significantly from the “classical”

Re: [deal.II] Different get_function_gradients

Hi Michael, > Does the (1) do all the work that the last two do? I.e., (1) gives the whole > chunk of gradient, and the last two extract the scalar and vector parts from > the chunk, respectively. > They are very closely related. (1) has no knowledge on the nature of the field(s) for which the

Re: [deal.II] Deal.II compilation error issued at example steps phase

Dear Saad, Well, after briefly browsing the Ginkgo repository online, its a bit of a mystery to me. The mask looks to be declared by this call to this piece of

Re: [deal.II] Marking Boundary Nodes in Triangulation

Hi Chen, Is the DoFTools::extract_boundary_dofs() function perhaps what you’re looking for? Best, Jean-Paul > On 23. Jun 2021, at 10:08, Chen R wrote: > > Hi all, > > So I am trying to

Re: [deal.II] Linear operators - TrilinosWrappers

Hi Feilpe, Firstly, I agree that the documentation is very light on details on how to use the linear operators with Trilinos linear algebra types. We can definitely improve this, and any contributions to that effect would be greatly appreciated! Right now, I can direct you to a few of the

<    1   2   3   4   5   6   >