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

Dear Everyone, During installation of deal.ii 9.2 through SPACK I receive the following error on ADOL-C: 1 error found in build log: 549-- Could not find a sufficient ADOL-C installation: Possible symbol clash between the ADOL-C library and Trilinos' SEACASChaco detecte

Re: [deal.II] Stabilized FEM implementation with bubble function

Hi W, You are right. I am not bounded by the memory at this moment. The only potential issue I can think of at this stage is the computation of the error norm. Technically, if I feed the slots in global RHS corresponding to bubble support as zero, the solution at the bubble support should

Re: [deal.II] AlignedVector set at index

Zachary, Ideally I would like to create an AlignedVector with a begin and end and then push_back only to those indices within those bounds. Is there a way to accomplish this? for (auto x=begin; x!=end; ++x) my_vec.push_back (*x); I assume the push_back function works similar to

Re: [deal.II] initializing values using random numbers (of doubles in an interval)

On 1/6/21 9:12 AM, Karthi wrote: Question (1) Should I use VectorTools::project or VectorTools::interpolate for initializing the variables? In practice, it does not matter very much, but interpolate() is much cheaper than project(). Both the initial conditions are not giving the desired

Re: [deal.II] PETSC SparsMatrix initialization error

On 1/6/21 8:07 AM, Zachary 42! wrote: Here is my code again with hopefully better comments. Looking at the loop structure for the dynamic_sparsity_pattern should make it clear. I think the 3rd argument in the constructor for the dynamic sparsity pattern is wrong if I want these dimensions.

Re: [deal.II] Stabilized FEM implementation with bubble function

On 1/7/21 6:58 PM, Lixing Zhu wrote: I implemented the local elimination as you suggested and it works. However, this ends up with many extra "redundant" DOFs in the global system. For example, if I am using N*N linear elements for a scalar 2D problem with structured mesh, instead of (N+1)^2

Re: [deal.II] Re: Slepc Solver Dimension error

Hi Zachary, It is difficult for many of us to find the time to debug problems from our users like this, I hope you understand. May I suggest that you look into this problem in more detail by yourself first? Have you used a debugger to see where the error occurred exactly? Do you understand what

Re: [deal.II] Stabilized FEM implementation with bubble function

Dear Jean-Paul, Thanks for your suggestion. I'll try to test with these two approaches. Regards, Lixing On Friday, January 8, 2021 at 4:36:44 AM UTC+8 Jean-Paul Pelteret wrote: > Hi Lixing, > > Another tutorial that might be of interest to you to look at is step-44 >

Re: [deal.II] Stabilized FEM implementation with bubble function

Dear W. Many thanks for your suggestion. I guess I have to implement a new class for my own bubbles at some point. A follow-up question on (1). I implemented the local elimination as you suggested and it works. However, this ends up with many extra "redundant" DOFs in the global system. For

[deal.II] Re: Slepc Solver Dimension error

Sorry, forgot the file! -- 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

[deal.II] Slepc Solver Dimension error

Hi everyone, I get a “sum of local lengths 32 does not equal global length 128” PETSC error in the following test file. It dies at the eigensolver.solve step. I tried to use the test and tutorial programs but to no avail. Still wrapping my head around the logic so please excuse me. Thank

Re: [deal.II] Stabilized FEM implementation with bubble function

1. I looked into the step-51 of the tutorial. It does illustrate a paradigm of segregating local DOFs and global DOFs. If I utilize this paradigm, the workflow would solve the local DOFs first (virtual node of the bubble function), which is a block matrix since bubble support from adjacent

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

[deal.II] deal.II Newsletter #148

Hello everyone! This is deal.II newsletter #148. It automatically reports recently merged features and discussions about the deal.II finite element library. ## Below you find a list of recently proposed or merged features: #11492: Make one variable 'const'. (proposed by bangerth)