Re: [deal.II] Re: element birth using FE_Nothing and FE_Collection in time-dependent problem with adaptive mesh refinement

Just a shot in the dark: Will it work if you set your FE_Nothing elements to be dominating? https://www.dealii.org/developer/doxygen/deal.II/classFE__Nothing.html#a69940158d653b46222943f09c6b25a4c Marc On Thursday, January 18, 2024 at 12:09:02 PM UTC-7 luis@gmail.com wrote: > Hello

[deal.II] Re: How to get a vector graphic (like SVG) in dealii?

On Tuesday, December 12, 2023 at 6:39:31 AM UTC-7 masou...@gmail.com wrote: I also tried to write a Python code for this purpose using vtk and matplotlib libraries. It works almost fine, but when I want to plot the mesh too, there are problems. A while ago I used matplotlib to create vector

[deal.II] Re: question about hybrid MPI and TBB thread parallelism

Hello Tim, > A simple timing tells me the assembly process of system-matrix takes 99% of the whole running time in every newton iteration. Just to make sure: did you compile the deal.II library and your code in Optimized mode/Release mode

[deal.II] Re: question about hybrid MPI and TBB thread parallelism

Hello Tim, > Yet, even though it is universally believed to be superior in terms of convergence properties, it is not widely used because it is often believed to be difficult to implement. One way to address this belief is to provide well-tested, easy to use software that provides this kind of

[deal.II] PhD Intern numerical simulation for additive manufacturing process @ hp, Palo Alto, California

Hello, I just stumbled over this internship posting for PhD students as it explicitly mentions the deal.II library. You can find more details here: https://jobs.hp.com/jobdetails/18492550/phd-intern-numerical-simulation-for-additive-manufacturing-process-palo-alto-ca/ Best, Marc -- The

Re: [deal.II] Re: Regarding tbb error

: DEAL_II_FORCE_BUNDLED_THREADS On Friday, April 28, 2023 at 2:05:36 PM UTC-6 Marc Fehling wrote: > Hello Jau-Uei Chen, > > the configuration did not pick up the TBB version you provided, and > instead uses the version that comes bundled with deal.II (i.e. tbb 2018U2). > However, the linker still tries

Re: [deal.II] Re: Regarding tbb error

sion of dealii 9.4.0. >> With this version I am getting the error as mentioned in the previous >> email. >> >> >> >> >> Thanks & Regards >> Deepika >> >> On Tue, Sep 20, 2022 at 2:48 AM Marc Fehling wrote: >> >>> Hello Deepik

Re: [deal.II] Installation of Deal II hangs

Hello Claudia, have you considered using the prebuilt virtual machine image? A current deal.II release will be readily available for you without compilation. https://www.dealii.org/download.html https://www.math.clemson.edu/~heister/dealvm/ Of course, that doesn't solve the compilation problem.

[deal.II] Re: solution transfer for mesh adaptivity

Hi! Your data lives on quadrature points, correct? And I assume you are using the parallel::distributed::Triangulation since you are checking for locally owned cells. We have a class parallel::distributed::

[deal.II] Re: A problem about cylindrical coordinates.

Hi! The manifold should be attached only to the hull of the cylinder, but in your code snippet you attach it to all cells. This leads to problems when new vertices need to be created on the symmetry axis. After a look in the GridGenerator::cylinder()

[deal.II] Re: Dimension Mismatch

Hello Jack, the error message is quite descriptive and tells you that something is off with the number of components of your Function object: *The violated condition was: exact_solution.n_components == n_components* Check line 200 of your code: ExactSolution () : Function(dim+1)

Re: [deal.II] What's the best strategy to speed up assembly?

Hello Kev, just to make sure: did you compile the deal.II library and your code in Optimized mode/Release mode ? Marc On Thursday, December 8, 2022 at 2:35:22 PM UTC-7 Wolfgang Bangerth wrote: > > Kev, > > > I've done some layman's

[deal.II] Re: How would you document a step/code/module?

Hello Abbas, if you are willing the share your code with the deal.II community, you can submit your code to the deal.II code gallery. The documentation for each of the example programs use the same doxygen formatting as deal.II. https://www.dealii.org/code-gallery.html

Re: [deal.II] Re: 'make test' runs wrong

you install the Intel oneAPI suite? >>>> >>>> Marc >>>> >>>> On Friday, November 11, 2022 at 7:27:23 AM UTC-7 dengm...@gmail.com >>>> wrote: >>>> >>>>> Dear Marc: >>>>> I have tried the '

Re: [deal.II] Re: How to mesh Solid objects in step-54?

the point(...) is not on the manifold". I know what the error is > about and I have confirmed the vertices of the coarse mesh are on the > manifold. However, the error pertains. Are there any other sources that > this error might have happened? > Thanks again. > Himal.

Re: [deal.II] Re: 'make test' runs wrong

b/workflows/linux.yml' at >>> https://github.com/dealii/dealii/pull/14145/files, but I still can not >>> handle this problem. I can not even pass the 'make -j 2' command (the >>> terminal exit automatically when progress reaches 57%). Is it better to >>> cre

Re: [deal.II] Re: How to mesh Solid objects in step-54?

in the file was 2.0. So, I changed it to 3.0 and still I > have got the following error: > [image: Screenshot from 2022-11-11 10-39-03.png] > The surf.vtk file contains following mesh: > [image: Screenshot from 2022-11-11 10-43-01.png] > Thanks! > Himal. > > On Thu, Nov 10,

Re: [deal.II] Re: 'make test' runs wrong

> terminal exit automatically when progress reaches 57%). Is it better to > create a new virtual box environment without intel mkl and mpi lib? > > Best regards > Meng > > Marc Fehling 于2022年11月10日周四 05:15写道： > >> Hello, >> >> a

[deal.II] Re: How to mesh Solid objects in step-54?

Hi Himal, the vertices of your coarse mesh must be located on the geometry provided by the IGES file. It is hard to tell from that amount of information where exactly your error comes from. I advise you to draw a sketch of both your IGES cylinder and your coarse mesh. Usually you should find

[deal.II] Re: 'make test' runs wrong

Hello, at some point, the quicktests failed on our ubuntu CI when using Intel OneAPI in deal.ii 9.4.0, see also https://github.com/dealii/dealii/issues/13477 Do you face the same problem when working on an up-to-date master branch? We fixed it at some point, see

[deal.II] Re: Segmentation fault in Trilinos direct solver

Hello, I can't identify the problem right away, but here are some tips in debugging it: - Does the problem persist if you disable both h- and p-refinement? If so, than an issue might be that your data structures haven't been reinitialized after you updated your discretization. -

[deal.II] Re: Matrix-based p-multigrid preconditioning

Hi! The "Possibilities for extensions" section of step-75 gives you a few hints on how to achieve p-multigrid preconditioning with matrix-based methods: https://www.dealii.org/current/doxygen/deal.II/step_75.html#Solvewithmatrixbasedmethods If you look for an example implementation, the "hpbox"

[deal.II] Re: Dealii and ASPECT installation

Hello Ziqi, which version of deal.II are you using? >From your log, it looks like you are using Intel oneAPI TBB. We started to support oneAPI versions with deal.II 9.4.0. The guide you are mentioning refers to deal.II 9.0. Could you try our latest release? Marc On Monday, September 19, 2022

[deal.II] Re: Regarding tbb error

Hello Deepika, which version of deal.II are you using? >From your screenshots, it looks like you are using Intel oneAPI TBB. We started to support oneAPI versions with deal.II 9.4.0. Could you try our latest release? Marc On Sunday, September 18, 2022 at 6:18:51 AM UTC-6

[deal.II] Re: Add/remove some cells "one by one"

Olivier, we offer weighted load balancing through the p4est library. This allows you to change load balancing dynamically upon a criterion that you choose. However this does not (necessarily) mean that you can do it on a cell-by-cell level. Are you willing to compromise? :-) Marc On Friday,

Re: [deal.II] Matrix-free DG with h-adaptive mesh refinement

Hello! You are calling prepare_for_coarsening_and_refinement() twice which is not allowed. Thanks for reporting this on the mailing list -- we will improve the error message that you encountered. I spotted a few issues in the order of events of your program. Tutorials step-32, step-42,

[deal.II] Re: Different results for varying n.o. of processors when solving Piezoelectric PDE via PETSc

Hi Asmorom, do your results also differ if you run your problem again on the same number of processes, or do you get the same result? If they still differ, there is a nice explanation on why this happens in the FAQ on PETSc:

Re: [deal.II] Periodic bc for dim=1

Hello Praveen, the parallel::distributed::Triangulation class is not instantiated for dim=1, only for dim>1. The reason behind it is that p4est does not allow to create 1D domains. With any other Triangulation class you should be fine. Marc On Tuesday, February 8, 2022 at 1:01:02 AM UTC-7

[deal.II] Re: Question about "make_periodicity_constraints"

Hello there! There are many tests that demonstrate how to make periodicity constraints with scalar values finite elements. For example, have a look at tests/bits/periodicity_01.cc . Let us know if that helped! Marc

[deal.II] Re: Use of DoFRenumbering

Hi I've found application examples for parallel::distributed::Triangulation in our tests. See for example https://github.com/dealii/dealii/blob/master/tests/mpi/renumber_cuthill_mckee.cc As you specifically mentioned step-40, I think this is what you were asking for? Best, Marc On Monday,

[deal.II] Re: Problems installing dealii-9.3.2 (AVX)

Hi Daniel, to me it looks like deal.II has been configured fine. During the cmake configuration phase, we try to compile small code snippets to see if your system is actually capable of certain features, i.e., vectorization (AVX) support. This is nothing you need to worry about. Best, Marc

[deal.II] Re: Obtain degrees of freedom starting from the cell

You can convert your Triangulation::cell_iterators into a DoFHandler::cell_iterators. An example is shown here . Marc On Tuesday, November 30, 2021 at 3:30:49 PM

[deal.II] Re: Obtain degrees of freedom starting from the cell

Hi Bob, you should be able to create a `DoFAccessor` iterator object with your index in some way, and then call this particular function. Or you can iterate over all active cells, and then check for their index and store the corresponding dof indices in some container. In general, I would

[deal.II] Re: Obtain degrees of freedom starting from the cell

Hello Bob, I think the function DoFAccessor::get_dof_indices() is the one you are looking for. Its use is demonstrated in multiple tutorial programs, e.g. step-4

[deal.II] Re: Coarse Cell IDs Distributed Triangulation

Hi Jonathan, there is a unique way to identify cells, even in parallel::distributed::Triangulation objects! For this, have a look at the CellId class. You can create an object of this class from a cell accessor using this

[deal.II] Re: Mesh export, high order

Hello, yes, you can use high-order output using the `DataOut` class after setting the corresponding `DataOutFlags` flag *write_higher_order_cells* . The tutorials step-11,

[deal.II] Re: SolutionTransfer Error

t them another time with every call of read_checkpoint or write_checkpoint, which might cause the error in your case. Best, Marc On Wednesday, October 27, 2021 at 12:49:29 AM UTC-6 Marc Fehling wrote: > Hello Toni, > > this error occurs when you load a triangulation object, but do not >

[deal.II] Re: SolutionTransfer Error

Hello Toni, this error occurs when you load a triangulation object, but do not de-serialize all the data associated with it, before you attempt to write another checkpoint. This error tells you that you should check your de-serialization process. Please make sure all data has been

Re: [deal.II] Building p4est with no system-wide zlib

Hi Michał, maybe it will help to invoke `configure --help` on the source of your p4est tarball. Hopefully there will be instructions on how to connect custom zlib libraries -- I couldn't find anything in their README :( Best, Marc On Thursday, September 16, 2021 at 6:26:23 PM UTC-6 Wolfgang

[deal.II] Re: PETSC not found by configuration cmake

Hello Tom! If you have a closer look at the error message, you'll see which files our configuration is looking for: -- FIND_LIBRARY(PETSC_LIBRARY NAMES petsc libpetsc HINTS / PATH_SUFFIXES lib lib64 lib) -- FIND_PATH(PETSC_INCLUDE_DIR_ARCH petscconf.h HINTS / PATH_SUFFIXES petsc include

[deal.II] Re: Subdivided cylinder and boundary ids

By the way, this function has just been introduced in April with #12003 <https://github.com/dealii/dealii/pull/12003>. Suggestions for improvement on this new feature are welcome! Best, Marc On Friday, June 25, 2021 at 12:37:31 PM UTC-6 Marc Fehling wrote: > Hi Alex, > > >

[deal.II] Re: Subdivided cylinder and boundary ids

Hi Alex, > I scrolled through grid_generator.cc and found that subdivided_cylinder has a hardcoded epsilon of 1e-5 for determining whether a face is at the boundary. Of course, with a system of my size, all faces will end up being considered to be at the boundary. > > Because it not only sets

Re: [deal.II] time dependent navier stokes (gallery) does not work

Hi Giang, > I did the git bisect (14 iterations) and it reveals that the commit 0c59a54d7385a8e862fa114201e91c0f3c53 gives the problem. For convenience, I linked to the commit with that specific hash. This

[deal.II] Re: Convergence rate when the analytical solution is unknown

is, this may require lots of memory. Marc On Friday, June 4, 2021 at 2:56:52 PM UTC-6 Marc Fehling wrote: > Hello Naren! > > you can use the `FEFieldFunction` > <https://www.dealii.org/current/doxygen/deal.II/classFunctions_1_1FEFieldFunction.html> > > class for this pu

Re: [deal.II] Re: Unable to Build Dealii on WSL Ubuntu 18.04

wrote: > Hello, > > Error after "make install". > > Regards, > Sanket Kulkarni > > On Fri, May 21, 2021 at 7:21 PM Marc Fehling wrote: > >> Hi Sanket, >> >> Yes, you need to install g++. >> >> For ubuntu, I would e

Re: [deal.II] Re: Unable to Build Dealii on WSL Ubuntu 18.04

KULKARNI wrote: Hello, I am attaching the screenshot. Do i need to install g++ On Fri, May 21, 2021 at 7:00 PM Marc Fehling <mailto:mafehl...@gmail.com>> wrote: Hi Sanket, what do you get if you type in `g++` in your console? Marc On 5/21/21 4:58 PM, SANKET KULKA

[deal.II] Re: Unable to Build Dealii on WSL Ubuntu 18.04

Hello, from your error message it appears to me that you have an `openmpi` wrapper installed on your system which doesn't know anything about the `gcc` compilers. Are you sure you have them installed? Marc On Friday, May 21, 2021 at 12:15:56 PM UTC-6 sanketku...@gmail.com wrote: > Hello, >

[deal.II] Re: initial velocity

Hello, you provided the `initial_values` and `initial_velocity` functions as lambda functions in your `Coefficients` class. You create an object of this class in the main function and then store a `SmartPointer` to it in an `OfflineData` object, which in turn is stored as a `SmartPointer` in a

Re: [deal.II] Re: Problem related to make_hanging_nodes_constraints in the hp cases

Hello Chong, you've set `dominate = true` for your `FE_Nothing` element. This means that your solution will be constrained to be continuous on interfaces with other elements. Since `FE_Nothing` represents the zero space, all DoFs on interfaces will be constrained to zero. You can disable this

[deal.II] Re: Removing cells in parallel distributed triangulation

Hello Reza, `parallel::distributed::Triangulation` objects store the entire coarse mesh, see also the module documentation . This allows for a quick execution of adaptive refinement. So from my understanding, you need to

Re: [deal.II] On mesh refinement & solution transfer with Raviart-Thomas

Hi Charlie, to 2) Yes, this update would be necessary. It should be okay to overwrite the .output file with the actual output of the test after you applied 1). to 3) I now wonder to where to place it as well! I'm surprised that we have lots of tests for SolutionTransfer in the hp context or the

Re: [deal.II] On mesh refinement & solution transfer with Raviart-Thomas

Hi Charlie! It looks like you found the cause for the issue! The transferred solution now looks like one would expect. It seems like this is a general problem with the restriction matrices in the Raviart-Thomas elements. Would you mind writing a patch for the deal.II library so that your

Re: [deal.II] How to refine mesh in different field

Hi Chen, step-46 gives you a good example on how to perform grid refinement when coupling different kinds of equations in different parts of the domain. I would highly suggest reading this

Re: [deal.II] Problem in the installing dealii-9.2.0

s Fehling, which version of Petsc is compatible with dealii? > > Best regards, > > Em sáb, 27 de mar de 2021 às 8:50 PM, Marc Fehling > escreveu: > >> Hello! >> >> > I wonder whether the PETSc interfaces have changed again. >> >> Sigh, I guess e

Re: [deal.II] Problem in the installing dealii-9.2.0

Hello! > I wonder whether the PETSc interfaces have changed again. Sigh, I guess exactly this happened, see https://github.com/petsc/petsc/commit/569ea7c476928c4ccab5493c49720b88ce3320f4#diff-10d5cd204c8b9d29ef9bc54965f06c2ac039d4cfcdaab6a01ce027a345b9d232 I'll prepare a fix. @Emmanuel: To

Re: [deal.II] Re: Transfer vector of solutions

Hello Karthik, it is perfectly reasonable to treat refinement for the initial mesh separately. I noticed that both your refine and coarsen fractions always add up to 100%. This is not a requirement! You can adjust both fractions independently until you are okay fine with the results. Marc

Re: [deal.II] Re: Transfer vector of solutions

inement at time zero makes > sense but as the solution decays I was hoping the mesh would coarsen (but > it refines further). I am clearly doing something wrong. I need some help > in fixing this issue. > > > Thank you! > > > Karthi. > > On Mon, Jan 25,

[deal.II] Re: deal.ii installation on NERSC Cori

Hello I am not familiar with the details of the NESRC Cori machine. In its documentation, I found the following manual . I hope this helps. I can only speak from my experience on HPC machines, and we had

[deal.II] Re: Transfer vector of solutions

Hi Karthi, if you work on the same DoFHandler, one SolutionTransfer object is sufficient. There are already member functions that take a container of solutions like yours as a parameter. Have a look at SolutionTransfer::prepare_for_coarsening_and_refinement

Re: [deal.II] Re: Parallel distributed hp solution transfer with FE_nothing

to use the `p::d::CellDataTransfer` class for your use case as described in the last message. Marc On Thursday, December 31, 2020 at 6:02:00 PM UTC-7 Marc Fehling wrote: > Hi Kaushik, > > Yes, this is possible by changing a cell from FE_Nothing to FE_Q using > p-refinement. &g

Re: [deal.II] Re: Parallel distributed hp solution transfer with FE_nothing

all time steps. I was > hoping to save some computation time, by only forming a system consists of > cells that are in the "active" layers only. > > Please let me if this makes sense? Is there any other method in deal.ii > that can simulation such a process? > Thank you

Re: [deal.II] Re: Parallel distributed hp solution transfer with FE_nothing

The FiniteElementDomination logic in the codim=0 case would indeed make up a cheap a priori check in this context. In case a FE_Nothing has been configured to dominate, the solution should be continuous on the interface if I understood correctly, i.e., zero on the face. I will write a few

Re: [deal.II] Re: Parallel distributed hp solution transfer with FE_nothing

once. I don't know if we have such a test for the general SolutionTransfer class. I will check that. Marc On Monday, December 28, 2020 at 1:39:33 PM UTC-7 Wolfgang Bangerth wrote: > On 12/27/20 8:48 PM, Marc Fehling wrote: > > > > 2) I did not know you were trying to interpol

Re: [deal.II] Re: Parallel distributed hp solution transfer with FE_nothing

_Q > elements after p-refinement? Are those always set to zeros after the > refinement? > > [image: image.png] > Thank you, > Kaushik > > On Wed, Dec 23, 2020 at 5:35 PM Marc Fehling wrote: > >> Hi Kaushik, >> >> Be careful on what you are doing here: Y

Re: [deal.II] Re: Parallel distributed hp solution transfer with FE_nothing

nk you Mark. >>> I am using the dealii lib that I got from apt-get from >>> deal.ii-9.2.0-backports. >>> I used PETSc and the abort was on even 1 cpus. I tried 2, 3, 6 cpus and >>> all aborted similarly. >>> >>> I will get the latest master b

Re: [deal.II] Re: Parallel distributed hp solution transfer with FE_nothing

>From your stacktrace I can see you are using PETSc and deal.II 9.2.0 which already incorporates the specified patch. Would you try to build the actual master branch anyways? On Wednesday, December 9, 2020 at 2:11:59 PM UTC-7 Marc Fehling wrote: > Hi Kaushik, > > I am unable to re

Re: [deal.II] Re: Parallel distributed hp solution transfer with FE_nothing

:TriaIterator > >, > std::allocator 2>::CellStatus, dealii::TriaIterator > > > > > const&, std::vector > > (dealii::TriaIterator >, > dealii::Triangulation<2, 2>::CellStatus)>, > std::allocator > > (dealii::TriaIterator >, > dealii::Triangulatio

[deal.II] Re: Periodic boundary conditions : Error using GridTools::collect_periodic_facepairs

Hi Aaditya, on first look your implementation looks good to me. Does the same error occur when you are using a standard `Triangulation` object instead of a `parallel::distributed::Triangulation`? As far as I know, the direction parameter does not matter for scalar fields (see also step-45).

Re: [deal.II] Re: Parallel distributed hp solution transfer with FE_nothing

Hi Kaushik, the `p::d::SolutionTransfer` class should be able to deal with `FENothing` elements in your example. The tricky cases are when you're coarsening a `FENothing` element with others as Bruno already pointed out (h-coarsening), or if you change a `FENothing` element to a different

Re: [deal.II] Re: DEAL.II INSTALLATION ERROR

folder, i.e., /path/to/petsc/include Best, Marc On Thursday, November 26, 2020 at 7:47:50 PM UTC-7 pushkar...@gmail.com wrote: > Yes I did follow the above instructions but I am still facing the same > issue . > > On Fri, Nov 27, 2020 at 3:53 AM Marc Fehling wrote: >

[deal.II] Re: DEAL.II INSTALLATION ERROR

Hi Pushkar! It appears the PETSc has been found during the configuration of deal.II with cmake, but you can not find the header files of the PETSc libraries during compilation. Did you follow all instructions on how to interface deal.II to PETSc on this particular guide

[deal.II] Re: Refinement on a parallel distributed triangulation with periodic BC

Hi Maurice! On Monday, November 16, 2020 at 8:09:14 AM UTC-7 maurice@googlemail.com wrote: > Looking at the doc of `collect_periodic_faces` (This function will collect > periodic face pairs on the coarsest mesh level of the given mesh (a > Triangulation >

Re: [deal.II] outer product of two vectors

Please have a look at this particular test which showcases how an outer product can be achieved with deal.II! https://github.com/dealii/dealii/blob/master/tests/full_matrix/full_matrix_57.cc Hope this helps! Marc Wolfgang Bangerth schrieb am Dienstag, 6. Oktober 2020 um 18:31:47 UTC-6: > On

[deal.II] Re: hp fem error assigning Fourier

Hi Ishan! You are correct: We opted for a more versatile approach in transforming solutions into Fourier or Legendre series with deal.II 9.2. Glad you figured it out! Marc -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see

Re: [deal.II] Re: Deal.ii installation

Hi Prasad! My guess now is the following: You have PETSc and Trilinos installed on your device, which deal.II finds, but complains that PETSc has been installed with a different MPI configuration as is available for deal.II. -- Include

Re: [deal.II] Re: Deal.ii installation

Prasad! If you look closely into your CMakeError.log file, you'll find that there are multiple tests failing. This is not a bad thing: deal.II figures out your system configuration this way and enables/disables certain features. However in your case, it seems that there is a mandatory test

Re: [deal.II] Re: Deal.ii installation

Hi Prasad! On Tuesday, April 28, 2020 at 10:54:54 AM UTC+2, Prasad Adhav wrote: > > I am using cmake version 3.10.2 > With version 3.10.2, it is less likely that my suggestion will fix your problem. Did cmake actually finished the configuration, as Wolfgang pointed out? Going through the

[deal.II] Re: Deal.ii installation

Hi Prasad! Thank you for providing the logs! It seems like this here is the cause of it: /usr/bin/ld: cannot find -lpthreads It seems like this is again related to the `-lpthread` problem. Just out of curiosity: Which version of cmake and which OS are you using? They may have updated CMake

[deal.II] Re: "libpthreads"?

Hi Victor! This has been fixed upstream with pull request #9117 . This issue only occurs with CMake version >= 3.16. Just compile from the master branch and you're good! Best, Marc On Wednesday, April 22, 2020 at 3:16:40 PM UTC+2, Victor Eijkhout

Re: [deal.II] Mesh refinement and the ability to transfer the data to the quadrature points of the new mesh on parallel::shared::triangulation.

Hi Alberto! If I understood you correctly, you transfer quadrature point data, with the `SolutionTransfer` class which is meant to transfer finite element approximations. A different class dedicated to the transfer of quadrature point data already exists: It is called

[deal.II] Re: Error in make_hanging_node_constraints() while using parallel::distributed::Triangulation with hp::DoFHandler

Hi Chaitanya, we turned your minimum working example into a test case for the deal.II library #9555 . Thank you for providing your code! Would you mind to give it a look, since we reduced and changed a few parts of it? Best, Marc -- The deal.II

[deal.II] Re: Error in make_hanging_node_constraints() while using parallel::distributed::Triangulation with hp::DoFHandler

Hi Chaitanya, This should've been fixed in #8365 , which is not included in deal.II-9.1.1. Compiling the most recent version of deal.II from the master branch should do the trick. Marc -- The deal.II project is located at http://www.dealii.org/

[deal.II] Re: Error installing with Candi

Hi David! On Saturday, October 12, 2019 at 4:34:22 AM UTC+2, David Ryan wrote: > > I'm trying to get deal.ii installed using candi on my Mac running macOS > Mojave. > > Everything seems to work up till the deal.ii compiling where it tells me > that it can't find the lapack libraries. > I've

Re: [deal.II] installation fails with intel/19.0.5

On Monday, September 30, 2019 at 11:01:45 PM UTC+2, Victor Eijkhout wrote: > > > > On Sep 30, 2019, at 3:23 PM, Marc Fehling > wrote: > > Victor, have you tried disabling C++17 support? Maybe that'll do the > trick... > > > cmake option please? > >

Re: [deal.II] installation fails with intel/19.0.5

On Friday, September 27, 2019 at 11:24:12 PM UTC+2, Wolfgang Bangerth wrote: > > Didn't we recently merge a patch where ICC reported that it understands > C++17, but doesn't in fact support this attribute? Does that ring a bell > for anyone? > Intel published a list of all C++17 features

Re: [deal.II] Parallel distributed hp solution transfer

Hi Doug! On Tuesday, August 27, 2019 at 3:41:11 AM UTC+2, Doug wrote: > > Thank you very much for the quick fix! Looking forward to pull this once > it goes through all the checks. > The patch has been merged. Let me know if this fix does the trick for you. We introduced a test named

Re: [deal.II] Parallel distributed hp solution transfer

Hi Doug! Hi Wolfgang! On Sunday, August 25, 2019 at 3:25:06 AM UTC+2, Wolfgang Bangerth wrote: > > On 8/23/19 6:32 PM, Marc Fehling wrote: > > > > Your scenario indeed revealed a bug: Currently, we set and send > > `active_fe_indices` based on the refinement fl

Re: [deal.II] Parallel distributed hp solution transfer

Hi Doug! Your scenario indeed revealed a bug: Currently, we set and send `active_fe_indices` based on the refinement flags on the Triangulation object. However, p4est has the last word on deciding which cells will be refined -- and in your specific scenario p4est makes use of it. I came up

Re: [deal.II] Parallel distributed hp solution transfer

Hi Doug! On Wednesday, August 21, 2019 at 4:00:49 AM UTC+2, Doug wrote: > > 8134: void dealii::parallel::distributed::SolutionTransfer VectorType, DoFHandlerType>::unpack_callback(const typename > dealii::parallel::distributed::Triangulation space_dimension>::cell_iterator&, typename >

Re: [deal.II] Parallel distributed hp solution transfer

Hi Doug, when dealing with distributed meshes, ownership of cells change and we may not know which finite element lives on cells that the process got recently assigned to. Thus, we need to transfer each cell's `active_fe_index`, which we do automatically during coarsening and refinement.

Re: [deal.II] Re: Heat equation (step-26): Negative values with small time step

Hi Wolfgang, On Monday, December 18, 2017 at 1:45:15 AM UTC+1, Wolfgang Bangerth wrote: > > I think your observation of negative values is an interesting one (and > surprising one, for many). Would you be interested in writing a couple of > paragraphs about time step choice for the introduction

[deal.II] Re: Heat equation (step-26): Negative values with small time step

Hi Bruno, I only heard about applying flux limiters on advection/convection problems, but not on diffusion-related ones. This conforms with what I recently found in literature, but I may skipped something crucial. The equation of interest is the heat equation:

[deal.II] Heat equation (step-26): Negative values with small time step

Dear deal.ii community! I stumbled over some interesting behavior of the heat equation from step-26. If I reduce the time step to a smaller value, let's say to 1e-6, I observe negative values for the solution near the sources (where gradients are large), which I would not expect. I guess it is

[deal.II] Convergence rate of solution scheme for incompressible Navier-Stokes equations

are changing over time. Why is the convergence rate in space inconsistent? Am I missing some crucial point? Best regards, Marc Fehling -- 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 receive