Re: [deal.II] fixing one component of solution to the same value

2017-04-12 Thread Wolfgang Bangerth
On 04/12/2017 05:03 PM, RAJAT ARORA wrote: I am not sure why such an error is occurring. This shows that sparsity pattern was generated that such a position (2249, 6) will not be filled, so no space was allocated while declaring the Petsc sparse matrix. However, constraints.distribute()

Re: [deal.II] fixing one component of solution to the same value

2017-04-12 Thread RAJAT ARORA
Hello Professor, I did what you suggested above. However, I am facing a problem. Just to repeat, I want to constraint all dofs on the z surface to a single dof (named as *the_dof = 2249*) I apply the constraints the following way. ComponentMask cm(dim, false); cm.set(2,

Re: [deal.II] Bug? - FETools::extrapolate changes behavior in 8.5

2017-04-12 Thread Wolfgang Bangerth
Hi Winni, when adjusting my own code to deal.ii version 8.5. I noticed a change in the behavior of the FETools::extrapolate method. Since I could not find this in the list of changes, I would like to let you know as this might have been unintended. Description: Until 8.4.1

Re: [deal.II] Calculating gradients at arbitrary points

2017-04-12 Thread Wolfgang Bangerth
On 04/12/2017 09:32 AM, David F wrote: Hello, I'm trying to calculate the symmetric gradients of a solution vector at an arbitrary set of points. However, so far I have only found the way to do it at the quadrature points by means of /fe_values[fe_extractor].get_function_symmetric_gradient/. I

[deal.II] Calculating gradients at arbitrary points

2017-04-12 Thread David F
Hello, I'm trying to calculate the symmetric gradients of a solution vector at an arbitrary set of points. However, so far I have only found the way to do it at the quadrature points by means of *fe_values[fe_extractor].get_function_symmetric_gradient*. I would like to know what is the proper

[deal.II] Bug? - FETools::extrapolate changes behavior in 8.5

2017-04-12 Thread Winnifried Wollner
Dear all, when adjusting my own code to deal.ii version 8.5. I noticed a change in the behavior of the FETools::extrapolate method. Since I could not find this in the list of changes, I would like to let you know as this might have been unintended. Description: Until 8.4.1 FETools::extrapolate

[deal.II] Thank you!

2017-04-12 Thread Timo Heister
I would like to follow up on the announcement of this release by publicly saying again how much we appreciate the many contributions from those who have sent it code, bug reports, fixed grammar and typos, or have helped in any other way. Many thanks! The ChangeLog lists at least the following

[deal.II] deal.II Version 8.5.0 released

2017-04-12 Thread Timo Heister
Version 8.5.0 of deal.II, the object-oriented finite element library awarded the J. H. Wilkinson Prize for Numerical Software, has been released. It is available for free under an Open Source license from the deal.II homepage at

Re: [deal.II] Re: fe_enriched and step-47

2017-04-12 Thread Denis Davydov
> On 12 Apr 2017, at 14:54, Sean McGovern wrote: > > >> Would the general point be that if one can go to higher resolutions with >> distributed parallelism that the enhanced accuracy of the enrichment >> functions can be compensated for? > > Sorry, i did not get what you

Re: [deal.II] Re: fe_enriched and step-47

2017-04-12 Thread Sean McGovern
> Would the general point be that if one can go to higher resolutions with > distributed parallelism that the enhanced accuracy of the enrichment > functions can be compensated for? > > > Sorry, i did not get what you mean. > > Regards, > Denis. > > Hi Denis, Thanks for your prompt answers. I

Re: [deal.II] Re: fe_enriched and step-47

2017-04-12 Thread Sean McGovern
Hi, Thanks to you all for the feedback. The ANS paper and the plans for sophisticated quadrature are great. One further point of clarification: I am new to XFEM, but have been interested in level sets for a while, in particular in the distributed parallel context. Do I understand correctly