Re: [deal.II] Loading solution in step 6

Ah, right, Thanks for the correction Wolfgang! > On 27 Feb 2019, at 06:58, Wolfgang Bangerth wrote: > > On 2/26/19 7:54 AM, Jean-Paul Pelteret wrote: >> >> Do you mean a 6-sided hexagon? Well, this is the sort of geometry that one >> might want to use a mesh genera

Re: [deal.II] Loading solution in step 6

Dear John, > Another simple question is that how to obtain the triangulation of a > standard hexagonal in 2d? Do you mean a 6-sided hexagon? Well, this is the sort of geometry that one might want to use a mesh generation tool (e.g. GMSH, Cubit) for since you can’t trivially create a

Re: [deal.II] Position of a given node

I wrote a small wiki entry for this this morning. You can see it here: https://github.com/dealii/dealii/wiki/Frequently-Asked-Questions#how-do-i-get-the-degree-of-freedom-indices-at-vertices

Re: [deal.II] Loading solution in step 6

Dear John, You don’t say how you write this "solution-…txt” file, so its difficult to say exactly what the best course of action is. Do you know about the block_write() and block_read()

Re: [deal.II] How to get the coordinates of a given degree of freedom on an edge?

Dear Kien, The problem that you’re facing (in the example that you gave us) is not one to do with deal.II itself, but rather C++ and how you’ve structured your code. When I build your code in debug mode, I see this warning: $ make Scanning dependencies of target problem [ 50%] Building CXX

Re: [deal.II] Parallel Nonlinear Poisson code in DEALII

Dear Boron, As a new user, I’d like to welcome to the forum and deal.II! Am I correct that you only see this issue when executing your program in parallel? I think that the problem results from the way that you initialise the solution vector. That is, this line: >

Re: [deal.II] How to compute convergence rate of L2 norm of error without exact solution? And how to compute convergence rate of in different golobally refinement?

Hi Chucui, I can only offer you a very brief reply right now: You might want to look into the ConvergenceTable class to help you with (1). The introduction to the class directs you to some tutorials that indicate its

Re: [deal.II] Trouble getting UMFPACK to work

Dear Darius, I just thought I’d ask whether you were able to solve your problem? Best, Jean-Paul > On 16 Feb 2019, at 05:22, Darius Barreto wrote: > > Dear Jean-Paul, > > Thanks for the reply. > > > > On Thursday, February 14, 2019 at 9:20:39 PM UTC+5:30,

Re: [deal.II] Trouble getting UMFPACK to work

alii/wiki/deal.II-in-Spack>). I hope that this clarify what you need to do. Best, Jean-Paul > On 14 Feb 2019, at 08:16, Jean-Paul Pelteret wrote: > > Dear Darius, > > Your installation of deal.II has not been compiled with UMFPACK supplied as > an external library. The bund

Re: [deal.II] Trouble getting UMFPACK to work

Dear Darius, Your installation of deal.II has not been compiled with UMFPACK supplied as an external library. The bundled version has also not been built, likely because you do not have the required dependencies installed. I believe that you need to have LAPACK installed for it to work, so you

Re: [deal.II] Automatic Differentiation of the Finite Element Jacobian

> The derivatives with respect to the > mesh points will be needed for optimization purposes. I’m not familiar with the algorithms used in shape optimisation, but one *possible* work-around if you need derivatives wrt. the mesh points on the finite element cell level would be to do the

Re: [deal.II] Automatic Differentiation of the Finite Element Jacobian

Dear Doug, > Looking at the current implementation in fe_values, mapping_q_generic, and > local integrators, I see that it looks relatively easy to template the value > type since it is already partially templated. > > Is there a reason why those functions were not templates? No, its likely

Re: [deal.II] Warning messages about text-based stub files not syncing with library files when compiling on Mac.

as A University > B.S. Nuclear Engineering > ––– > ––– > Cell:(682) 667-3646 > School Email: zach.ha...@tamu.edu <mailto:zach.ha...@tamu.edu> > Personal Email: zhardy1...@gmail.com <mailto:zhardy1...@gmail.com> > > -- > > > On Thu,

Re: [deal.II] Warning messages about text-based stub files not syncing with library files when compiling on Mac.

h.ha...@tamu.edu> > Personal Email: zhardy1...@gmail.com <mailto:zhardy1...@gmail.com> > > -- > > > On Wed, Jan 30, 2019 at 11:57 PM Jean-Paul Pelteret <mailto:jppelte...@gmail.com>> wrote: > Dear Zachary, > > I googled this problem and it seems

Re: [deal.II] Warning messages about text-based stub files not syncing with library files when compiling on Mac.

Dear Zachary, I googled this problem and it seems that, dependent on what version of Xcode and MacOS you’re running, there might be one of two solutions. Here’s a similar post on the Apple forums that outlines the two possibilities . The first

Re: [deal.II] command line compilation error on a linux machine using spack installation of deal.ii dev

Dear Nicola, This is pretty strange, since you say that the same code compiles on other machines. From what I can tell from a brief search of this error message, it typically appears when there is a stray brace or comma or the like. Are you completely sure that source files on this problematic

Re: [deal.II] output vector solution to vtk file

Dear Pham, Welcome to deal.II! Your observation is related to this frequently asked question . Basically, deal.II outputs the

Re: [deal.II] MPI processes: use only a part of the available processes for dealii

Dear Daniel, If appropriate, you could also consider using MPI_Comm_split() to create a new communicator based off of a subset of the world / master communicator, and pass that along to all classes / functions in your FEM code.

Re: [deal.II] How to implement GridGenerator::channel_with_cylinder?

glad to contribute and > will really appreciate your help. > > Nick. > > > On Wednesday, January 16, 2019 at 11:52:09 AM UTC-5, Jean-Paul Pelteret wrote: > Dear Nick, > > Indeed, it would be quite straightforward to extend the > GridGenerator::plate_with_a_hole()

Re: [deal.II] apply Pouseuille BC at multiple inlets with different velocity magnitudes

Dear Laila, I have two suggestions that may help you here: 1/ You could add the boundary ID as an argument to the PoiseuilleBoundaryValues class: template class PoiseuilleBoundaryValues : public Function { public: PoiseuilleBoundaryValues (const typename types::boundary_id b_id)

Re: [deal.II] How to implement GridGenerator::channel_with_cylinder?

Dear Nick, Indeed, it would be quite straightforward to extend the GridGenerator::plate_with_a_hole() function to one that has a filled inclusion using the steps that David has listed.

Re: [deal.II] Cooperation proposition

Dear Chinedu, Please be aware that the original post was deleted because it was unsolicited spam. Best, Jean-Paul On Wednesday, December 5, 2018 at 8:05:02 AM UTC+1, Chinedu Nwaigwe wrote: > > Yes, I am interested. > -- The deal.II project is located at http://www.dealii.org/ For mailing

Re: [deal.II] cmake linker errors using spack view

Dear Praveen, This appears to be a documented issue - please see the “Library conflicts” section of our CMake documentation . There may not be anything you can do about this. I build deal.II on some machines in exactly the same manner as

Re: [deal.II] Interpolated Boundary Conditions from Distributed Solution

Dear Kyle, I’m not too familiar with FEFieldFunction, but I’m aware of the suggested usage within the parallel::distributed context. I’ve taken a quick look into the source code for FEFieldFunction (in particular the compute_point_locations()

Re: [deal.II] Ordering of entries in FEValuesExtractors::SymmetricTensor

a vector base ,we need nomalize the > entry to be one. for the case of symmetric tensor , if the deviator part is > 1 , after normalization , it will become 1/sqrt(2). > so , what is the real implementation for the symmetric tensor 's base? > > 在 2016年11月29日星期二 UTC+8下午2:35:44，Jean

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

to know how to get A. > Or what do you think of solving this equation? > > best， > FU > > > 在 2018年11月15日星期四 UTC+8下午10:45:24，Jean-Paul Pelteret写道： > Hi, > > In general there is, to the best of my knowledge, no easy answer to your > question. Can you explain wh

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

> A is the sparse matrix that I want to get in the end, where N and M are > sparse matrices. > > What else can you do to solve this problem? > > Thank you > > FU在 2018年11月14日星期三 UTC+8下午2:32:03，Jean-Paul Pelteret写道： > Hi, > > To compute th

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

Hi, To compute the sparse matrix inverse M^{-1} you would have to use a direct solver such as SparseDirectUMFPACK . If this operation (M^{-1})N appears in the context of matrix-vector products, i.e. you’re actually

Re: [deal.II] How to restart a time dependent code from its breakpoint?

ss and headfile for parallel one. Shall I > rewrite my code into a parallel one? Thank you very much! > > Best, > Chucui > > 在 2018年11月5日星期一 UTC+8下午5:01:06，Jean-Paul Pelteret写道： > Dear Chucui, > > I don’t use this functionality myself, so I can only point you to the

Re: [deal.II] Modelling Anistropic materials for electrostatic simulations

Dear Philip, Let me offer an alternative to what Wolfgang has suggested. I don’t think that it is necessary to use a mixed formulation to introduce material anisotropy into your problem. If you go back to the derivation in our previous discussion, there was this line: (6) div [\epsilon_{0}

Re: [deal.II] How to restart a time dependent code from its breakpoint?

ental point of C++? Does it only work at > some projects such as the ASPECT library > <https://www.google.com/url?q=https%3A%2F%2Faspect.geodynamics.org%2F=D=1=AFQjCNEohyzhbZ5uZr1UJ6cdCXdI3mAj1w> > you mentioned above? Thank you very much! > > Best, > Chucui > > 在 2

Re: [deal.II] How to restart a time dependent code from its breakpoint?

Dear Chucui, The Triangulation , DoFHandler and Vector

Re: [deal.II] Question on step-33

Dear Mathieu, As Maxi already pointed out, step-33 is indeed not written to take advantage of any MPI capabilities. Perhaps somewhat unintuitively the call to Utilities::MPI::MPI_InitFinalize mpi_initialization(argc, argv, ::numbers::invalid_unsigned_int); in the main() function configures

Re: [deal.II] function interpolation (projection) to PETScWrappers::MPI::BlockVector mixed up after dof re-numbering

Dear Rochan, This sounds like a bug to me, but I must admit that I have no good guesses as to where the problem might be located. Would you be able to come up with a minimal working example that illustrates the issue so that we can try to determine the source of the problem? Also, what version

Re: [deal.II] Error in installing deal.ii with spack on CentOs7

@unibs.it <mailto:alberto.salvad...@unibs.it> > web-pages: > http://m4lab.unibs.it/faculty.html <http://m4lab.unibs.it/faculty.html> > http://dicata.ing.unibs.it/salvadori <http://dicata.ing.unibs.it/salvadori> > > > On Mon, Oct 22, 2018 at 8:43 AM

Re: [deal.II] Error in installing deal.ii with spack on CentOs7

Hi Alberto, As best I recall, their is a problem with a hidden clash in the dependencies between the latest PETSc and Trilinos versions. (I thought that there was an issue/PR on github/spack, but I couldn’t find it right now.) The former wants a new version of SuiteSparse while the latter does

Re: [deal.II] Configuration with Adol-C fails with wrong version

FpNF5rM/JLZgRu8sAQAJ> > there is link on development 2.6.4 version, and as I realized DEAL.II boost > conflicts with adol-c boost so the last one I should be set up wit flag > --without-boost > > With regards, > Oleg Kmechak > > On Monday, 8 October 2018 14:26:13 UTC+2, Jea

Re: [deal.II] Configuration with Adol-C fails with wrong version

Dear Oleg, Please see this recently closed issue on our GitHub repository: https://github.com/dealii/dealii/issues/7293 There have been some patches that have been recently submitted upstream to fix some bugs in the most recent release. You should

Re: [deal.II] Boundary condition from nodal traction values

Dear Qin, If I understand this statement > I have a vector named traction[dof_handler.n_dofs()], which gives the > various traction on each dof correctly, you have a pre-assembled global traction vector (of size n_dofs) that you wish to add to the system. Why is it that you want to

Re: [deal.II] Can GridIn read 3D quad mesh?

jppelteret/dealii/blob/4d6c06d4ebae1b342861d2cf937f0ce3bb331d5b/source/grid/grid_in.cc#L960> > > This scenario would currently raise an error. Would you think that would be > easy? I would be happy to provide sample files if that would be helpful. > > Best, > Shawn > >

Re: [deal.II] Repost Is the approach for the electrostatic Bi-linear form correct?

tion and I came across > this class: > > https://dealii.org/9.0.0/doxygen/deal.II/classDataPostprocessorVector.html > <https://dealii.org/9.0.0/doxygen/deal.II/classDataPostprocessorVector.html> > > Is this what you we referring to? > > On Wednesday, September 5,

Re: [deal.II] Can GridIn read 3D quad mesh?

> That's so fast! Well, it was an easy fix so we got lucky there. But the test case you provided was really useful, so we appreciate you providing that. > Thank you Jean-Paul :) You’re quite welcome. Best, Jean-Paul > On Thursday, September 20, 2018 at 1:20:46 AM UTC-7, Jean-Paul

Re: [deal.II] Update of Ubuntu

s is largely unrelated, since any (MPI dependent) deal.II program is likely to fail until you rebuilt deal.II and its dependencies. I hope that you get this issue solved! Best, Jean-Paul > On Thursday, September 20, 2018 at 3:02:08 PM UTC+2, Jean-Paul Pelteret wrote: > Dear Yaakov, >

Re: [deal.II] Update of Ubuntu

Dear Yaakov, > Recently, I have updated Ubuntu Where does your deal.II installation come from? If you’ve compiled it yourself, then have you rebuilt deal.II and its entries dependency list? With a distribution update often comes an update in the compiler as well as other dependencies. >

Re: [deal.II] Animation with gnuplot

A quick google search came up with this result: http://www.gnuplotting.org/tag/animation/ Does this answer your question? > On 20 Sep 2018, at 12:00, chandra sekhar yaswanth devarakonda > wrote: > > Can anyone tell me how to create a movie with a

Re: [deal.II] Can GridIn read 3D quad mesh?

that the ABAQUS > equivalent is "node set". Please let me know if you'd like me to try that! > > Best, > Shawn > > On Wed, Sep 19, 2018 at 10:00 PM Jean-Paul Pelteret <mailto:jppelte...@gmail.com>> wrote: >> Hi Shawn, > > >> Thank you for y

Re: [deal.II] Can GridIn read 3D quad mesh?

know about ABAQUS itself (I would have thought it would be possible). I’ll try to modify the input file you supplied, which should be possible since the geometry is simple enough. Best, Jean-Paul > > On Tue, Sep 18, 2018 at 10:57 PM Jean-Paul Pelteret <mailto:jppelte...@gmai

Re: [deal.II] Can GridIn read 3D quad mesh?

Dear Shawn, The GridIn::read_abaqus() function is currently not equipped to deal with the codimension 1 case. I might be able to find some time to look at this in the next few days though. In addition to the mesh file that you’ve already given us, may you please construct and give us a fully

Re: [deal.II] Compilation error in dealii development version (9.1.0-pre, shortrev c6b7876) using gcc/7.2.0 compiler and cuda/9.2

Hi Sambit, Your new error looks like the same one reported (https://github.com/dealii/dealii/issues/7204 ) yesterday so I don’t think that using another version of CMake will fix it. I’m not quite sure what the answer to the problem is, since the

Re: [deal.II] Abbreviation of deal-II

Dear Chandra, The answer to your question is summarised on our webpage (see https://dealii.org/about.html ): deal.II — a name that originally meant to indicate that it is the successor to the Differential Equations Analysis Library — is a C++ program library

Re: [deal.II] Compilation problem with Development sources

Dear Michał, I think that calling git checkout is insufficient. You probably want to call git pull (or git fetch && git rebase origin/master) to fetch and merge the most current state of the remote repository. Best, Jean-Paul > On 17 Sep 2018, at 15:38, Michał Wichrowski wrote: > > I've

Re: [deal.II] hermite element

Dear Chandra, There’s a module that lists the finite element spaces and associated finite elements that have been implemented in deal.II. It looks like only a discontinuous Hermite FE

Re: [deal.II] Repost Is the approach for the electrostatic Bi-linear form correct?

Dear Philip, It looks to me like you’re mixing things up a little bit. In all of the equations that you presented there is a mixture of governing equations, assumed relationships and the consequence of combining them. From what I can tell you’re actually trying to use the same equation more

Re: [deal.II] implementing traction neumann boundary

Hi Shaqir, So, I have a short question to clarify something: By commenting out the calls to VectorTools::interpolate_boundary_values() and MatrixTools::apply_boundary_values(), are you implying that you do not impose equation 4 (the Dirichlet boundary conditions) anywhere on the boundary of

Re: [deal.II] Design question about multiscale FEM with deal.ii

Hi Kendra, > Thank you, Jean-Paul, for the hint. Unfortunately this is not really what I > am looking for. I found a solution (that I don't like) by iterating through > two meshes on different levels simultaneously. Once I find a good solution I > will post it here and link to the code. I

Re: [deal.II] how to remove rows from system matrix

Hi, > I started with the code published in > https://journals.ub.uni-heidelberg.de/index.php/ans/article/view/22317 and > adopted it for my problem. Since you had removed Step-47 tutorial it is hard > to figure out how to use FE_Enriched class. FYI. It looks like step-47, which was removed in

Re: [deal.II] Design question about multiscale FEM with deal.ii

Hi Konrad, > Do you think it would be possible to (ab)use the multigrid classes for that? > As far as I understand there are dofhandlers on each level of the mesh. Can I > get dofhandlers for all active subcells of say coarse cell number n in level > 5? The DoFHandler::begin_active()

Re: [deal.II] step-8

Dear Chandra, For me, one of the most instructive textbooks on linear elasticity (and that fundamentals of FEM) is this one: @Book{Hughes2000a, title = {The Finite Element Method: Linear Static and Dynamic Finite Element Analysis}, publisher = {Dover Publications Inc.}, year =

Re: [deal.II] projecting values of the solution on the boundary

ing going > wrong. > > In fact one could also move the nodes a bit to generate "almost" zero > volumes. In some cases it can be done easily, but in general this approach is > unfeasible for complex meshes. > > I appreciate your comments. > Best > > Alber

Re: [deal.II] projecting values of the solution on the boundary

ion worth having on github rather than here. Best regards, J-P > > Il giorno lunedì 13 agosto 2018 21:47:29 UTC+2, Wolfgang Bangerth ha scritto: > On 08/08/2018 02:03 AM, Jean-Paul Pelteret wrote: > > Hi Alberto, > > > > I have dealt with a similar problem where

Re: [deal.II] make run error: "dyld: Library not loaded"

Dear Zhibo, I guess that this could be fixed in Candi (? ping Uwe, Timo). As of the El Capitan release (I think), its not possible to use DYLD_LIBRARY_PATH anymore if SIP (System Integrity Protection) is enabled. For an immediate solution I would suggest trying to use install_name_tool to add

Re: [deal.II] implemention of periodic boundary condition

Dear Chucui, Thanks for the explanation - clearly I was wrong with regards to exactly how much complexity could be stripped from this example. Anyway, I think that Daniel has made the critical observation as to where the fault may lie. Hopefully his suggestion will sort this out for you.

Re: [deal.II] How to use intellisense to parse dealii code using visual studio remotely?

Hi there, I’m really happy to hear that you managed to work this out for yourself, and thanks a lot for posting your solution! I’m sure that will be very useful to our other Visual Studio users in the future. Best, Jean-Paul > On 08 Aug 2018, at 17:07, 程迪 wrote: > > > Hi everyone, > > I

Re: [deal.II] implemention of periodic boundary condition

Dear Chucui, Thanks for giving this a go. This is better than before, but I think that you still misinterpret what a minimal working example is. This is a 1000 line piece of code, but I estimate only 100 lines of code are required to reproduce the problem. If the problem arises in the first

Re: [deal.II] UMAT subroutine in dealii.

Dear Mahesh, I guess this is less of a “deal.II” question and more of a C++ question. It is possible to call a Fortran function from C++ code - if you do a google search for “c++ calling fortran routine” you get numerous helpful links showing how it would be done in concept, e.g.

Re: [deal.II] implemention of periodic boundary condition

Dear Chucui, I think that what Daniel was asking for was the absolute minimum size test case that reproduces the error in building the sparsity pattern. I interpreted (and he’s welcome to correct me) that wants to help you find out why you hit this error: An error occurred in line <1764> of

Re: [deal.II] projecting values of the solution on the boundary

Hi Alberto, I have dealt with a similar problem where I needed to transmit solution values between two different problems that shared a common interface. If I remember correctly, the way that I did this was to precompute the positions at which I would need the solutions from problem 1 for

Re: [deal.II] How to give dof per each node?

Dear Chandra, Its a little unclear to me as to what you mean when you say > to give dof per each element/node Do you mean that you assign DoFs to the mesh? If so, then what Rakesh has sketched out for you gives an outline of the steps needed to do this are. This is mirrored in the second

Re: [deal.II] UMAT subroutine in dealii.

Dear Mahesh, If you browse through both the tutorials and the code-gallery contributions , you will notice that deal.II is completely agnostic about which finite element

Re: [deal.II] implemention of periodic boundary condition

Hi Chucui, I’ll try to answer some of your questions quickly: > 1. Why the source code of function "GridTools::collect_periodic_faces" use > "cell_iterator" rather than "active_cell_iterator”? Because your mesh should be periodic from the coarsest representation all the way down to the finest

Re: [deal.II] implemention of periodic boundary condition

Dear Chucui, Since you haven’t provided a minimal working example, I could only speculate what might be the issue. From the incomplete example that you’ve shown here, I would guess that this might be the problem: > On 05 Aug 2018, at 09:28, chucui1...@gmail.com wrote: > >

Re: [deal.II] Modifying the mesh topology in Triangulation object

Dear Rakesh, I believe that the move_mesh() function in step-18 demonstrates exactly this procedure. Does this help? Best, Jean-Paul > On 02 Aug 2018, at 20:51, Rakesh wrote: > > Hello All, > > I have

Re: [deal.II] Getting RHS values at nodes with DBC

; > > On Wednesday, August 1, 2018 at 10:08:12 AM UTC-4, Jean-Paul Pelteret wrote: > Hi Rajat, > > Could you perhaps assemble a second RHS vector with an empty constraint > matrix, and then extract the entries in this vector that are constrained by > the

Re: [deal.II] Nedelec Elements for unstructured adaptively refined meshes

> On 10 May 2018, at 22:41, Bruno Turcksin wrote: > > If I remember correctly, J-P couldn't make Nedelec work for his problem a few > months ago. Yes, that’s correct. I was working on a magneto-static problem in 3d, and wasn’t able to achieve sensible results. That

Re: [deal.II] Deprecated .h file compressed_sparsity_pattern.h

Dear Henry, The DynamicSparsityPattern family of classes has superseded many of the older “Compressed” sparsity patterns that were deprecated in deal.II 8.3 (or thereabouts). The header file that you now need to

Re: [deal.II] set manifold for merged boundaries

Hi Weixiong, Great, I’m glad that you managed to fix your problem! Best, Jean-Paul > On 29 Mar 2018, at 23:47, Weixiong Zheng <zwxne2...@gmail.com> wrote: > > Both idea work fine now. Thank you so much! > > 在 2018年3月29日星期四 UTC-7上午11:31:23，Jean-Paul Pelteret写道： > Hi Wei

Re: [deal.II] set manifold for merged boundaries

> #7 9 libdyld.dylib 0x7fffe1c02255 start + 1: 9 > libdyld.dylib 0x7fffe1c02255 start > > > make[3]: *** [CMakeFiles/run] Abort trap: 6 > make[2]: *** [CMakeFiles/run.dir

Re: [deal.II] set manifold for merged boundaries

Hi Weixiong, The issue that you’re seeing here is actually related to a discussion that we’re currently having regarding the behaviour of boundaries and manifolds (and will hopefully be fixed soon in the developer version). You’re hitting the issue discussed here

Re: [deal.II] Any easy way to output cellwise average solution

Dear Weixiong, There’s an interpolation function that should do this for you: https://dealii.org/8.5.1/doxygen/deal.II/namespaceVectorTools.html#a5e3af70a47cedfaf361cf5c621e94e3d The

Re: [deal.II] FE w/ Gauss Point DoFs

Dear Jonathan, Is FE_DGQArbitraryNodes perhaps what you’re looking for? From the documentation: > The primary purpose of this class is to provide an element for which the mass > matrix can be made diagonal by choosing

Re: [deal.II] Force and PK1

Dear David, > I studied the step- 44 Tutorial Program and read already a lot about how to > apply a Neumann contribution. But I have still two questions: > In my case, I like to apply a force (not a force per unit reference area) as > a Neumann contribution on the RHS. But for the Integral, Ni

[deal.II] Re: Solution value with respective vertex coordinate

Dear David, This is quite a frequently asked question (well, some permutation of it), and so I refer you to the list of answers for this search query on the forum instead of repeating the discussion here. This

Re: [deal.II] Tutorial fail to run in my enviroment and i need help ....

To comment further on what Bruno’s said, it looks like you probably need to rebuild and reinstall deal.II again. This will ensure that the latest version of these libraries are detected during the configuration stage. Best, Jean-Paul > On 05 Jan 2018, at 16:14, Bruno Turcksin

Re: [deal.II] "Someone is tinkering with a vacant arena slot" Error

Dear Mustafa, You haven’t provided us with very much information to help you here. It would appear that you are hitting this assertion within the TBB library but

Re: [deal.II] installation process on Ubuntu

ckout a different version/branch. Please let me > know if this is the case. > I guess I could also try with a newer Trilinos version, what version should > work? > > El viernes, 17 de noviembre de 2017, 11:28:31 (UTC-5), Jean-Paul Pelteret > escribió: > Dear Juan Carlos, &

[deal.II] Re: installation process on Ubuntu

Dear Juan Carlos, Can you please tell me what the commit hash of the deal.II respository that you're working from is? I did briefly introduce a bug that showed similar symptoms but I thought that I fixed it. It might also have to do with your version of Trilinos, which appears to be older than

Re: [deal.II] Creating periodic Boundary Conditions after doing a uniform subdivision

Dear Phani, Great, thanks for being willing to contribute a patch! I’ve created a new article on our wiki that should hopefully explain these steps in sufficient detail. Best, Jean-Paul > On 11 Nov 2017, at 08:15, Phani Motamarri

Re: [deal.II] Creating periodic Boundary Conditions after doing a uniform subdivision

Dear Phani, The problem, as the error message suggests, is that you’re calling “triangulation.add_periodicity()” after you’ve performed some (global) refinement. You can only call this function on the coarsest mesh, and therefore you typically only need to do so once after you’ve created the

Re: [deal.II] Re: VectorTools::interpolate fails with MappingCartesian

Dear Praveen, I’ve opened up an issue on the GitHub repository for you: https://github.com/dealii/dealii/issues/5402 Kind regards, Jean-Paul > On 05 Nov 2017, at 11:07, Praveen C wrote: > > I forgot to attach the error I get.

Re: [deal.II] CellDataStorage with mesh refinement

nt? data_quadrature is the quadrature > rule I'm using in my model and projection_fe has to be of an order high > enough to interpolate this data correctly, am I at least right there? > > Best regards, > Frederik > PS: Thanks, I forgot GeometryInfo at this point > > Am Freitag

Re: [deal.II] Compilation error during the use of CellDataStorage datastructure

Hi Phani, Hmm… interesting. I think that you may be using the class in a way that we did not expect - I think that we intended it to store a struct/class rather than a primitive data type. I think that this should be an easy fix, namely that this assert (and perhaps others) need to be extended

Re: [deal.II] CellDataStorage with mesh refinement

Hi Frederik, I’m being wildly presumptuous here (so feel free to say if I’ve presumed incorrectly), but unless you have a special set of data that needs to be considered then the ContinuousQuadratureDataTransfer

[deal.II] Re: extract the mesh information

Dear Peng, The first few tutorial programs demonstrate how to query basic information related to the mesh. The documentation for the triangulation class and

Re: [deal.II] Using values from the time step before while solving a heat equation fails

the heat equation? The boundary values for the update is > zero, and the boundary values for the solution is always a given value (for > constant boundary conditions). > > Am Mittwoch, 18. Oktober 2017 18:16:26 UTC+2 schrieb Jean-Paul Pelteret: > boundary_constraint

Re: [deal.II] Re: Using values from the time step before while solving a heat equation fails

boundary_constraints.distribute(distributed_solution); boundary_constraints.distribute(distributed_old_solution); This looks suspicious to me. Don't the Dirichlet constraints change between time-steps? On Wednesday, October 18, 2017 at 4:16:08 PM UTC+2, Maxi Miller wrote: > > Removing the

Re: [deal.II] Reflect filter in ParaView

Dear Chenchen, With the way that we output data, Paraview cuts each quadrilateral into two triangles for the purpose of visualisation. So what you see is a linear interpolation of the data, and not a bilinear interpolation (which I’m guessing you were expecting if you used FE_Q elements).

Re: [deal.II] EigenValue Problem

Dear Nitish, I’m not familiar with this tutorial and I don’t have much to suggest, but I have one basic question: When increasing the polynomial order of the finite element, do you increase the corresponding quadrature order as well? Regards, Jean-Paul > On 19 Sep 2017, at 14:53, Nitish Anand

[deal.II] Re: Design of output of computationally expensive parameter function

Dear Dominik, I believe that your question is directly along the lines of this recent discussion: https://groups.google.com/forum/#!topic/dealii/FW9Jy5_CUGY Does the information contained in that post help you? Regards, Jean-Paul On Thursday, September 14, 2017 at 10:54:01 AM UTC+2, Dominik

[deal.II] Re: Object still in use by other objects

Presumably in relation to this ? On Tuesday, September 12,

Re: [deal.II] SynchronousIterators compilation errors

Judging by Wolfgang’s comment, I’m probably mistaken on this point! > On 11 Sep 2017, at 18:33, Jean-Paul Pelteret <jppelte...@gmail.com> wrote: > > Dear Claire, > > In the next version of deal.II we will be requiring the compilers be C++11 > compatible. We therefor

<    1   2   3   4   5   6   >