[deal.II] Re: Error in GridIn::read_abaqus for a mesh with two blocks

ned in the old version of Cubit that >> you have. >> >> >> >> According to Cubit 15.2 documentation, 'gen_input_file' is an input file >> with the given file name that will be generated when the Sculpt Parallel >> command will be executed. This is a

[deal.II] Re: Development history of deal.II

Very cool :-) Is the tool that you used to make this freely available? On Friday, October 21, 2016 at 1:29:15 AM UTC+2, Wolfgang Bangerth wrote: > > > > All, > a colleague of mine showed me a neat way to visualize the development > history of a code this week. I had to apply this to deal.II, of

[deal.II] Re: How to apply boundary values for a particular point on the boundary instead of the whole boundary surface

Thanks for sharing your solution, Stephen! On Saturday, October 22, 2016 at 1:43:29 PM UTC+2, Stephen DeWitt wrote: > > In case that link dies, here's the snippet of code: > > > // Set constraints to pin the solution if there are no Dirichlet BCs for a > component of a variable in an elliptic equ

Re: [deal.II] Re: Development history of deal.II

Awesome, that looks to be it. Thanks! Cheers, J-P On Saturday, October 22, 2016 at 4:03:28 PM UTC+2, Matthias Maier wrote: > > > > Very cool :-) Is the tool that you used to make this freely available? > > The tool is gource if I am not absolutely mistaken. > > http://gource.io/ > > Best, >

Re: [deal.II] Re: Development history of deal.II

Thanks for pointing that out -- I didn't look in the comments and so I missed that useful command! For anyone who's interested, I had gource added to spack yesterday. On Saturday, October 22, 2016 at 10:19:36 PM UTC+2, Wolfgang Bangerth wrote: > > On 10

[deal.II] Re: Error in GridIn::read_abaqus for a mesh with two blocks

----- > > An error occurred in line <899> of file > > > in function > > void dealii::GridIn<3, 3>::read_abaqus(std::istream &, const bool) > [dim = 3, spacedim = 3] > > The violated condition was: > &g

[deal.II] Re: Error in GridIn::read_abaqus for a mesh with two blocks

t; > On Monday, October 24, 2016 at 11:05:36 AM UTC-4, Jean-Paul Pelteret wrote: >> >> Dear Oded, >> >> The error message from the development version is more useful. >> Unfortunately there's a bug somewhere in the routine that appears from >> time-to-tim

[deal.II] Re: Error in GridIn::read_abaqus for a mesh with two blocks

baqus file format. > > Thanks, > Oded > > On Monday, October 24, 2016 at 4:39:57 PM UTC-4, Jean-Paul Pelteret wrote: >> >> Hi Oded, >> >> Hmm... I didn't have much luck with the script either. That was >> unexpected! I've dug up an older vers

[deal.II] Re: Accessing solution in a vertex given by global index found with find_closest_vertex

Hi Krei, This has been a popular topic as of late. Here's a list of posts that will help you achieve what you're looking to do. Regards, J-P On Tuesday, October 25, 2016 at 4:08:43 PM UTC+2, krei wrote: > > Hi, I

[deal.II] Re: Accessing solution in a vertex given by global index found with find_closest_vertex

ween the "global vertex index" given by > find_closest_vertex and the solution vector index. > > On Tuesday, October 25, 2016 at 5:13:38 PM UTC+3, Jean-Paul Pelteret wrote: >> >> Hi Krei, >> >> This has been a popular topic as of late. Here's a list of po

[deal.II] Re: Error in GridIn::read_abaqus for a mesh with two blocks

ched are the journal file, the UCD file and the resulting screenshot > from ParaView. > > > > Thank you! > > Oded > > > > P.S. > > I suggest you to upload the journal file to the wiki page – it may be > useful to other people. > > On Tuesday, Oct

[deal.II] Re: Modeling Stokes flow in and around a droplet using the approach of step-46

Hi Oded, As you've observed, this class is not yet compatible with the hp::DoFHandler. If you'd like to try to extend it to support this case then we'd be happy accept a patch and to help you in doing so. Regards, J-P On Friday, October 28, 2016 at 11:54:27 PM UTC+2, Oded Yaakobi wrote: > > He

[deal.II] Re: deal.II on youtube and twitter

Awesome, thanks very much for sharing Praveen! On Wednesday, November 2, 2016 at 2:27:36 PM UTC+1, Praveen C wrote: > > Dear all > > I sometimes upload videos/pictures of solutions obtained with deal.II on > youtube and twitter. I am tagging with #dealii on twitter. > > https://twitter.com/search

[deal.II] Re: deal.II on youtube and twitter

Awesome, thanks very much for sharing Praveen! On Wednesday, November 2, 2016 at 2:27:36 PM UTC+1, Praveen C wrote: > > Dear all > > I sometimes upload videos/pictures of solutions obtained with deal.II on > youtube and twitter. I am tagging with #dealii on twitter. > > https://twitter.com/search

[deal.II] Re: Problems in handling the userdata of Distributed Triangulation

Dear Jack, You'll be happy to know that this functionality has recently been added to the developer version of deal.II via the CellDataStorage and TransferableQuadraturePointData

[deal.II] Re: Help in value_list

Dear Atanu, The RightHandSide class in step-7 derives from the Function class. This in turn defines the value_list

[deal.II] Re: Help in value_list

You're welcome. Since you say that you're a C++ beginner, I'm not sure how much you know about polymorphism: just in case you need some further clarification, here's an example that highlights what's actually going on here with the base and

[deal.II] Re: Help in value_list

On Thursday, November 3, 2016 at 9:42:23 PM UTC+5:30, Jean-Paul Pelteret > wrote: >> >> You're welcome. Since you say that you're a C++ beginner, I'm not sure >> how much you know about polymorphism: just in case you need some further >> clarification, he

Re: [deal.II] Re: Convergence problem arising for large number of DoFs

Dear Hamed, >From the code you've provided you don't seem to impose a dirichlet boundary condition on any "y" degree of freedom. You need to fix at least one degree of freedom in each coordinate direction in order for your problem to have a unique solution. In some cases the CG method can solve

Re: [deal.II] Re: using .vtk file as initial condition

Dear Anup, Probably the simplest way to find examples of the save and load functionality is to look through deal.II's unit tests . Here are a list of tests that are related to the save functionality $ grep --include=\*.{h,cpp,cc} -rnw '.' -e

[deal.II] Re: Question about Step 44 shape function gradient computations

Dear Claire, You are right that this is a total Lagrangian formulation but that doesn't mean that one is restricted to defining the problem in terms of fully referential quantities. One can arrive at the same conclusion from a number of starting points, but ultimately its because we'd chosen

[deal.II] Re: Question about Step 44 shape function gradient computations

> > Le lundi 7 novembre 2016 19:35:18 UTC+1, Jean-Paul Pelteret a écrit : >> >> Dear Claire, >> >> You are right that this is a total Lagrangian formulation but that >> doesn't mean that one is restricted to defining the problem in terms of >> fully refe

[deal.II] Re: Learning Finite Element Method

Dear Divya, Its incredibly difficult to answer your question because the answer is a deeply personal one. Instead of embarking on some long diatribe, I will simply encourage you to do your own research to find a set of resources that work for you. Understandably the lecture videos provided b

Re: [deal.II] Re: Interfacing dealii-8.4.1 with trilinos-12.6.2 - make test failed

Dear Joaquin, It been a while since I've compiled trilinos manually, but if nothings changed since version 11.4 then your should be able to specify which MPI trilinos should be compiled against using the following cmake parameters: -D TPL_ENABLE_MPI:BOOL=ON \ > -D MPI_BASE_DIR:PATH=$DIR_BASE/Op

Re: [deal.II] coordinates of dof (Q1)

Dear Julian, In general its not easy to provide the functionality that you're looking for because not every degree-of-freedom is associated with a support point (e.g. the FE_DGPMonomial element) and let alone a vertex (e.g. an FE_Q of polynomial order 2 has DoFs with support points at face and

[deal.II] Re: Question on Vectortools::integrate_difference.

Dear Jaekwang, The role of the additional quadrature formula is to ensure that your function is not evaluated at "superconvergent" points, which would lead to an underestimation in the error associated with the finite element solution. I have little knowledge of functional analysis, so apart fro

[deal.II] Re: storing sparse matrix in .txt and scan in Matlab

Dear Anup, > 1. if the way of writing the sparse matrix on a .txt file is appropriate > for post-processing in matlab, if not, please suggest > a way; > There are three print functions in the deal.II SparseMatrix. As explained in the documentation

Re: [deal.II] problem with user_index()

The user_index is a placeholder for you, the user, to store some data. By default, all user_indices are set to zero until you set a value with cell->set_user_index() . On Sunday, November

[deal.II] Re: Very basic question on setting boundary id

Hi Jaekwang, With what you've shown I'm going to guess that your problem might simply be related to numerical / round-off error associated with floating point calculations. if ( cell->face(f)->center()[0]==1 ) // The LHS value is a computed double, > and the RHS value is an integer which will

[deal.II] Re: Very basic question on setting boundary id

> > } > > } > > > > <https://lh3.googleusercontent.com/-vxXjt_6106w/WCjEWZXkl9I/A9M/HtP_S_FZR3wqpaBSHpe7SjgNG9tCnYX5gCLcB/s1600/Screen%2BShot%2B2016-11-13%2Bat%2B1.50.56%2BPM.png> > > > > 2016년 11월 13일 일요일 오후

[deal.II] Re: mesh quality during deformations using MappingQEulerian

Hi Thomas, Trilinos' Mesquite package provide mesh optimisation tools. I'd be happy to see this integrated into deal.II, so if you decide that this would help you solve your problem then I'd be willing to help you to do so. The challenge here would be t

[deal.II] Re: mesh quality during deformations using MappingQEulerian

Hi Tom, Yes, you're essentially correct. With Mesquite I believe one would essentially redefine a triangulation in whatever format it wants (I think you provide nodal coordinate and connectivities), tell it which nodes remain fixed (i.e. provide some constraints) and then simply specify some a

[deal.II] Re: mesh quality during deformations using MappingQEulerian

Hi Tom, An alternative to implemented a complex mesh update strategy is to define an initial mesh that, when mapped via your Euler vector, produces a higher quality mesh. For example, you could remove

[deal.II] Re: Periodic boundary condition in Meshworker.

Sudarshan, Does MeshWorker take a constraints matrix? I suspect that it does . If so then yes, you should be able to implement periodic boundary conditions. I

[deal.II] Re: a small issue with step44

True, the condition at the bottom should probably read AssertThrow > > > (newton_iteration < parameters.max_iterations_NR, > ExcMessage > > ("No convergence in nonlinear solver!")); Thanks, I'll make a note to fix this. On Thursday, November 17, 2016 at 1:22:54 AM UTC+1, Anup Basak wrote:

[deal.II] Re: mesh quality during deformations using MappingQEulerian

> > Yes, this is what I am struggling to understand, you have put it very > simply. > Great! I'm glad that gave you a bit more insight. I only hope you didn't take that as a patronising comment :-) > I can see how this is trivial for linear FE_Q's, and these three links > should make this

[deal.II] Re: How to change material ids during iteration.

Dear Kyusik, Cell iterators have the function set_material_id , so you can change the material ID of any cell at any time. If by change the boundary of 2 materials you mean physically

[deal.II] Re: How to change material ids during iteration.

> > First of all, Thank you very much for your kind reply. > You're welcome :-) > So, is it better for me to study step-28 than step-18? > > Or, both examples would help for me to solve my problem(change of boundary > between two different materials during iteration)? > I'm not 100% certain

[deal.II] Re: mesh quality during deformations using MappingQEulerian

Hi Thomas, Yes, I agree in that the online documentation (like in many other Trilinos packages) is not up to scratch. Have you had a look at the user guide

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

Dear Arjun, The ordering of the extracted data should align with how the components are internally unrolled. Both the Tensor and SymmetricTensor classes have two functions to assist with this type of operation, namely unrolled_to_component_indices

Re: [deal.II] Accesing cell from cell index

Dear Kyusik, Two points: 1. Please create a new thread for a new question. Even though yours might be related to this one, this question has been asked and answered and its definitely tangential to the discussion here. 2. You have asked a question wanting general commentary on code that you hav

Re: [deal.II] Re: mesh quality during deformations using MappingQEulerian

Hi Tom, Its great that you're getting somewhere with this! I have taken some good advice from Oded (it turns out we work one building > over from each other!), > What a happy coincidence :-) > I'm trying hard to maintain the MappingQEulerian approach in my current > model setup, and I want

[deal.II] Re: constraining solution in a part of the domain

Dear Anup, You'll probably just have to manually add entries to your ConstraintMatrix using add_line and set_inhomogeneity

[deal.II] Re: deal.ii and Reduced Basis method

Dear Giulia, Judging by the lack of reply it seems that no-one knows of (or has implemented) an open-source library to specifically do model order reduction with deal.II. However, at least one of my colleagues uses deal.II to perform ROM via proper orthogonal decomposition, so deal.II definitel

[deal.II] Re: Apply Neumann Boundary condition

Dear Benhour, As I explained to someone else the day before yesterday, the point of this forum is not to provide a code review service without you having done due diligence on your s

Re: [deal.II] Re: mesh quality during deformations using MappingQEulerian

Hi Tom, > J-P, you (and the entire dealii community) are a paragon of encouragement, >> and I appreciate that very much. >> > Thank you, your kind words do mean a lot! > here's what I came up with yesterday > That's a handy piece of code to keep around. Thanks for sharing it! This works

[deal.II] Re: Apply Neumann Boundary condition

Dear Benhour, Thank you for taking my point heart and providing more specifics about your actual problem. For the horizontal axis, I have considered the zero displacement along > y-axis because of symmetry. Is it correct? Should I fix the center of the > circle as well? > If you don't provid

[deal.II] Re: Distorted cell when it is crated by gmsh

Dear Kyusik, The error message that deal.II has outputted explains what the problem is: > However, this is an internal line not located on the boundary. You cannot > assign a boundary indicator to it. You cannot give a boundary indicator value to internal lines/faces. If appropriate, you can

Re: [deal.II] Re: Apply Neumann Boundary condition

Dear Benhour, > These are 9 concepts, in 9 lines of text, all without any kind of > elaboration. > In other words, we have no real idea what you are doing, and equally > important, we have no idea what is actually wrong. When debugging problems > like yours, you need to learn to reduce thing

Re: [deal.II] Re: Modeling Stokes flow in and around a droplet using the approach of step-46

Hi Oded, Is this the recommended way to address point a), or is there another way > which is better? > I don't work primarily in fluid-mechanics, so this is a bit out of my field of knowledge (and comfort zone) and I'm going to defer to other's wisdom here. I suppose you're uncomfortable that

Re: [deal.II] Re: Modeling Stokes flow in and around a droplet using the approach of step-46

Dear Oded, I believe that Wolfgang has answered your question in much more detail and much more precision than I would have, so I think that there's no need for me to elaborate on further on what I said. Regards, J-P On Monday, December 5, 2016 at 5:47:17 PM UTC+1, Oded Yaakobi wrote: > > Hi W

[deal.II] Re: GitHub: make developers team public and add @dealii to ping all

No way to move it, but I can delete it if you want? On Tuesday, December 6, 2016 at 5:41:32 PM UTC+1, Denis Davydov wrote: > > wrong forum, sorry, i meant to post it here > https://groups.google.com/forum/?fromgroups#!forum/dealii-developers > If there is a way to move the post, feel free. > > On

[deal.II] Re: GitHub: make developers team public and add @dealii to ping all

Oh, actually Wolfgang could probably move it - I don't have write access to the developer forum. On Tuesday, December 6, 2016 at 5:47:36 PM UTC+1, Jean-Paul Pelteret wrote: > > No way to move it, but I can delete it if you want? > > On Tuesday, December 6, 2016 at 5:41:32 PM UTC

Re: [deal.II] Re: GitHub: make developers team public and add @dealii to ping all

gt; On 12/06/2016 09:49 AM, Jean-Paul Pelteret wrote: > > Oh, actually Wolfgang could probably move it - I don't have write access > to > > the developer forum. > > Nope, can't do -- completely separate mailing lists... > > -- > --

Re: [deal.II] Re: GitHub: make developers team public and add @dealii to ping all

Maybe you also don't have rights on the other group, ha ha :-) On Tuesday, December 6, 2016 at 6:26:27 PM UTC+1, Wolfgang Bangerth wrote: > > On 12/06/2016 10:14 AM, Jean-Paul Pelteret wrote: > > Don't you have the little "move" button at the top of the post (se

Re: [deal.II] general question

Dear Benhour, > > 2-Deviatoric strain > > dev_T_elastic = deviator(T_elastic); > > > > dev_T_elastic = T_elastic - (1/dim) * trace(T_elastic) * (Tensor<2, > > dim>(StandardTensors::I))/3.0; > > The expression '1/dim' is computed in integer arithmetic and evaluates to > zero. You want to w

Re: [deal.II] deal.ii workflow?

Dear Michael, It could definitely be overwhelming if you dive straight into an advanced application using deal.II, and is certainly not advised! Your options for information resources are listed in this post . A general recommendatio

[deal.II] Re: general question

> > For 2D (plane strain, plane stress and axisymmetric) and 3D the 1.0/3.0 > should be used (not 1/dim). > Good point! Special consideration definitely must be made for the axisymmetric / plane stress / plane strain cases. Cheers, J-P -- The deal.II project is located at http://www.deali

[deal.II] Re: Tensor output via ConditionalOStream pcout

Dear S. A. Mohseni, >From the class documentation , I figure that you should be able reach the equivalent output with (something like) the following steps: if (pcout.is_active()) { system_rhs.print(pcout.get_stream());

[deal.II] Re: Tensor output via ConditionalOStream pcout

Yes, well TrilinosWrappers::SparseMatrix only has a print() function. I took from what you wrote about that you were using dealii:: SparseMatrix. As could be expected, you sho

Re: [deal.II] square pyramid boundary

Dear Loylick, For these more complex grids that are not defined in GridGenerator, you have a few options. You can: 1. Use other GridGenerator functions in conjunction with the GridGenerator::merge_triangulations

Re: [deal.II] square pyramid boundary

Dear Loylick, Great, I'm glad you've managed to create a grid using the merge function. The way that you've specified the boundary id's is typically fine. The other function sets the boundary id's of both the faces and edges, which is unnecessary for most cases. Regards, Jean-Paul On Tuesday,

[deal.II] Re: BlockSparseMatrix and BlockVector with a nested hierarchy of blocks

Dear Oded, I think that the BlockMatrixArray class might offer the functionality that you're looking for? If so then I'm not quite sure what the equivalent class for the BlockVector would be. If there is none, then you h

[deal.II] Re: Tensor output via ConditionalOStream pcout

Dear S.A. Mohseni, You need to look carefully at the documentation for the specific SparseMatrix class. Its unclear to me whether LA::MPI::SparseMatrix represents a TrilinosWrappers::SparseMatrix

[deal.II] Re: ifem

Dear Joaquin, The warning appears because you are calling a function that is marked as deprecated, so although it remains in the library for now it will be removed in the near future. There are two a

[deal.II] Re: Tensor output via ConditionalOStream pcout

Dear Seyed, Thanks for sharing your print function. For future reference, the FullMatrix class has two print

Re: [deal.II] Re: ifem

Dear Joaquín, >From your previous posts I assume that all of PeFT, Pe[qs] and F[qs] represent rank 2 tensors. If this is correct then no, the contraction indices that you've specified are not correct. They're out of the range [0,1] that are valid for rank-2 tensors, and this would presumably tr

[deal.II] Re: deal.II cmake stuck at "Performing Test DEAL_II_HAVE_USABLE_FLAGS_DEBUG"

Dear Praveen, Unfortunately I don't have any suggestions about CMake hanging. I've just reinstalled deal.II through Spack on both my Apple and Linux machines and I didn't run into the problem you're experiencing. Denis has some experience with installing deal.II though spack on a number of diff

Re: [deal.II] Re: Problem with installation

Dear Sudarshan, I see that you've posted this question already in as a new topic . Since you've made that post, lets rather keep keep the discussion there. Regards, Jean-Paul On Friday, December 16, 2016 at 1:27:24 PM UTC+1, Sudarshan

[deal.II] Re: make test

Dear Sudarshan, Could you please post the other log files, namely "detailed.log" and "summary.log". If the experience of previous posters with MPI problem is anything to go by then you've either (1) got a broken MPI installation, or (2) have compiled one or more of deal.II's dependencies using a

Re: [deal.II] Re: ifem

Dear Joaquin, As a supplement to Luca's message, with every release of deal.II there is the possibility that we introduce some incompatibilities that break existing code. This is what must have happened here - I assume that you've upgraded from version 8.2.1 or earlier to something newer, and t

[deal.II] Re: Computation of B-operator

There are additional checks, warning and error messages that are processed in debug mode. So they're specifically in place to help detect problems like you're having. For this reason it is *highly* recommended to do your development in debug mode. This is mentioned numerous times in the documen

[deal.II] Re: Computation of B-operator

Its no problem. I'm glad you found the bug! On Monday, December 19, 2016 at 10:49:55 AM UTC+1, mohsengineer...@gmail.com wrote: > > Got it ;) Thanks to Jean-Paul Pelteret I found the problem. I tried to get > the jacobians without being inside the cell loop, silly and embarrassing

[deal.II] Re: internal compiler error: recipe for target 'source/numerics/CMakeFiles/obj_numerics.debug.dir/all

Dear ExcitingDDT, It may be possible that you've run out of system resources on the VM during compilation. How much RAM have you allocated to this VM image, and are you compiling with multiple threads? Regards, Jean-Paul On Monday, December 19, 2016 at 2:17:03 PM UTC+1, exciting...@gmail.com

Re: [deal.II] Re: BlockSparseMatrix and BlockVector with a nested hierarchy of blocks

Hi Oded, Have you managed to get anywhere with this? I have no particularly creative suggestions for this issue. The only thing I can think of is to derive a new class from SparseMatrix, which would contain your BlockMatrixArray as a member variable, and work with that. There is also the BlockL

[deal.II] Re: deal.II FE structure

Dear Seyed, You may also be interested in looking at the module on the interplay between finite elements and mappings , and the description to the mapping class. Note that th

[deal.II] Re: Evaluating the gradient of the solution at the centroid of the element

Firstly, a tangential observation: For the mean dilatation element (Q1-P0-P0), the dilation and pressure variables can be condensed out on an element level and a modified one-field formulation remains. For this element the "mean dilatation" means that one must compute \bar{\theta} = 1/V_{e} *

[deal.II] Re: Evaluating the gradient of the solution at the centroid of the element

efficient than the dilatation element because for each element I just > need to calculate one quantity instead of solving a linear system to remove > the static condensation from the displacement system. > > On Thursday, December 22, 2016 at 11:22:40 AM UTC-5, Jean-Paul Pelteret

[deal.II] Re: deal.II FE structure

One option would be to create a single cell triangulation that matches the reference cell's geometry. Then you can get all of this information easily without having to dig into the inner working of the FEValues and Mapping classes. J-P On Friday, December 23, 2016 at 11:29:41 AM UTC+1, mohsen

[deal.II] Re: cout stiffness matrix in step-3

Dear Michael, There are a number of "print" functions in the sparse matrix classes. See the answer to this question , as well as some of the surrounding discussion for more details. Kin

[deal.II] Re: Relationships between material parameters

Dear Seyed, There is a table in wikipedia that details how to convert between fundamental elastic material parameters. Due to the number of permutations that are available, I would be hesitant to add this to deal.II itself. However,

Re: [deal.II] Re: Relationships between material parameters

Constitutive_Parameters::Bulk_modulus::from_lambda_and_mu(double lambda, double mu) This can get untidy quite quickly, but I can’t think of a much better alternative… Best, J-P > On 10 Jan 2017, at 14:07, Wolfgang Bangerth wrote: > > On 01/10/2017 04:21 AM, Jean-Paul Pelteret wrote: >> >

Re: [deal.II] Re: Relationships between material parameters

AssertThrow(kappa > 0.0, ExcMessage("Non-positive bulk modulus!")); return kappa; } On Tuesday, January 10, 2017 at 2:45:13 PM UTC+1, Wolfgang Bangerth wrote: > > On 01/10/2017 06:18 AM, Jean-Paul Pelteret wrote: > > Do you hav

Re: [deal.II] Re: Relationships between material parameters

I am most swayed towards Wolfgang's suggestion that takes in some identifying enumeration type and the associated material parameters. I think its probably the most transparent and easy to maintain option that we've come up with thus far. Seyed, would you like to try to implement this? Regards

[deal.II] Re: Discontiuities in material properties

Dear Bjorne, There are quite a number of examples of Laplace and Poisson problems in deal.II's tutorials . The picture that you've shown could well be produced by a derivative of step-6

[deal.II] Re: Relationships between material parameters

Dear Seyed, It would be great to have your contribution to deal.II! When you have the time to implement this, open up an issue on the GitHub page and we could discuss this further. We can always refer back to this conversation at that point in time

[deal.II] Re: interpolate_to_different_mesh produces strange results

Hi Krzysztof, I don't have a direct comment on the result that you've shown here, but I think it would be useful if you mentioned what type and order of FE you are using here. One thing that you should be aware of is that Paraview does not interpolate data in the same way as deal.II does, so t

Re: [deal.II] Re: interpolate_to_different_mesh produces strange results

will give a better (although still not exact) visualisation of the approximated solution field. J-P On Thursday, January 12, 2017 at 9:43:46 PM UTC+1, Wolfgang Bangerth wrote: > > On 01/12/2017 06:45 AM, Jean-Paul Pelteret wrote: > > > > That said, it does seem a bit str

Re: [deal.II] Examples makefile

Dear Victor, CMake has some facilities to add to the link line. Here are a couple of stackoverflow posts t

[deal.II] Re: Elastic energy computation

Dear Seyed, There are two get_strain() functions in step-18. The first computes the symmetric part of the shape function gradients, and the second computes the small strain tensor. To store the values of the en

Re: [deal.II] Re: Elastic energy computation

Dear Seyed, These two things are very different, and cannot be compared with one another. The “strain” in the assembly routine stems from the weak form, which is presented in the introduction to the tutorial. If one considers a variational approach, is coincides with the variation of the strain

[deal.II] Re: How to solve a linear system with two RHS and a same matrix?

Dear Pasha, At first glance, it would appear that there is a further ConstraintMatrix::distribute_local_to_global function designed with this in mind. It takes in the local_matrix but d

Re: [deal.II] Re: How to solve a linear system with two RHS and a same matrix?

ch iteration. > In other hand, the local matrix is not available for the second RHS during > matrix assemble process. > > Regards > Pasha > > On Tuesday, January 17, 2017 at 12:12:06 PM UTC+3:30, Jean-Paul Pelteret > wrote: > Dear Pasha, >

[deal.II] Re: Computation of shape function gradients at the reference cell

Dear Tulio, This is the type of information that you can get directly from the chosen FiniteElement, e.g. this function provides the shape function gradient on the unit cell. You could also

Re: [deal.II] Re: How to solve a linear system with two RHS and a same matrix?

Pretty much exactly as I said. You first need to store the total displacement vector, but this would be the displacement from the spatial to the material configuration, accumulated by subtracting your incremental update from this global vector. You then create an Eularian map

Re: [deal.II] Re: How to solve a linear system with two RHS and a same matrix?

values); > > It works when I run the code with one CPU but I got the following error when > I run it with more than one processor > > [1]PETSC ERROR: Caught signal number 11 SEGV: Segmentation Violation, > probably memory access out of range > > Regards, > Pasha > &

Re: [deal.II] Re: How to solve a linear system with two RHS and a same matrix?

Dear Pasha I think that there are two problems with how you've created your Eulerian mapping: dealii::TrilinosWrappers::MPI::BlockVector > rel_solution2(partition_relevant); > Here you've provided the partitioning, but not an MPI communicator. You probably only need to use the locally own

Re: [deal.II] Re: How to solve a linear system with two RHS and a same matrix?

kVector5printERSojbb+0x126) > > [0x7fd462513f60] > [pasha:03624] [21] > ./FEM(_ZN3FEMILi2EE26assemble_large_deformationEb+0x2327) [0x560017] However, maybe this is not the case and some else has some other ideas. On Wednesday, January 18, 2017 at 10:37:48 AM UTC+1, Jean-Paul Pelteret wrote:

Re: [deal.II] Re: How to solve a linear system with two RHS and a same matrix?

ndler, rel_solution2.block(0)); > > > Thank you > Pasha > > On Wednesday, January 18, 2017 at 1:07:48 PM UTC+3:30, Jean-Paul Pelteret > wrote: > Dear Pasha > > I think that there are two problems with how you've created your Eulerian > mapping: > > d

Re: [deal.II] Re: How to solve a linear system with two RHS and a same matrix?

Great, thanks a lot. I've added it to the issue on GitHub that I opened this morning. So clearly the problem is elsewhere. Have you tried passing a normal MPI::Vector to the Eulerian map instead of the MPI::BlockVector that you had used in the pr

<    1   2   3   4   5   6   7   8   >