Re: [deal.II] How to compute derivatives of the composite function with the reference cell mapping

2023-08-05 Thread Anton Evgrafov
Thank you Wolfgang, this does look relevant! /Anton On Sat, Aug 5, 2023 at 12:37 PM Wolfgang Bangerth wrote: > > On 8/5/23 04:30, Anton Evgrafov wrote: > > I have a problem involving a > > field living in R^{dim^2} and a field in R^{dim}. > > This may or may not be relev

Re: [deal.II] How to compute derivatives of the composite function with the reference cell mapping

2023-08-05 Thread Anton Evgrafov
lized for dimensions 1..3. Right now I am doing a "quick and dirty" test implementation, where I am basically manually working on a tensor product mesh with tensor product FEM etc, which is highly inelegant. If you are doing/interested in research, I would be happy to discuss the prob

Re: [deal.II] How to compute derivatives of the composite function with the reference cell mapping

2023-08-04 Thread Anton Evgrafov
I was contemplating is to apply the NonMatching::DiscreteQuadratureGenerator on a fake triangulation with just one element. Then I simply would need to interpolate (evaluate) my real-space function at this one element and let the class' implementation do the magic for me. This is probably cheaper and

[deal.II] How to compute derivatives of the composite function with the reference cell mapping

2023-08-03 Thread Anton
Hello, I would like to use the NonMatching::QuadratureGenerator class directly, as I need to generate multiple (as in many!, unfortunately) nonmatching quadratures for the same cell. Therefore I would like to avoid the route of interpolating the level set function onto the whole mesh etc, as

Re: [deal.II] Re: Eigenproblem and creating a preconditioner out of a linear operator, Eigensolver Selection

2021-09-17 Thread Anton Ermakov
Dear Daniel, I am working on a similar problem (planetary acoustic oscillations). I am interested in looking at your tutorial for an electromagnetic cavity, but it seems that the pull request was deleted. I wonder if it can be revived. Thank you, Anton. On Monday, July 8, 2019 at 3:51:08 AM

Re: [deal.II] Vector eigenvalue problem with SLEPc

2021-08-25 Thread Anton Ermakov
state* *[0]PETSC ERROR: Matrix is missing diagonal entry 3* I presume this happens due to a zero block in one (or both) of the original matrices. Hmmm, should I then just explicitly see all zero diagonal values to 0.0. Anton. On Wednesday, August 25, 2021 at 12:20:38 PM UTC-6 Wolfgang Bangerth

[deal.II] Vector eigenvalue problem with SLEPc

2021-08-25 Thread Anton Ermakov
nomial filtering transformation. Thanks a lot, Anton. -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal

[deal.II] Re: shift-invert transformation problem with slepc

2021-08-24 Thread Anton Ermakov
Ok, I see - it is that EPS_TARGET_MAGNITUDE flag that was missing. It works now! Thanks a lot, Simon. Anton. On Tuesday, August 24, 2021 at 12:14:20 AM UTC-6 simon...@gmail.com wrote: > Hi, > > I think you need to set both which eigen pair and the target value. For > example,

[deal.II] shift-invert transformation problem with slepc

2021-08-23 Thread Anton Ermakov
]PETSC ERROR: Shift-and-invert requires a target 'which' (see EPSSetWhichEigenpairs), for instance -st_type sinvert -eps_target 0 -eps_target_magnitude Any help to resolve this error is greatly appreciated. I also attached the full code of modified step-36. Thank you, Anton. -- The deal.II project

Re: [deal.II] treatment of handing nodes in the heat conduction problem

2021-03-21 Thread Anton Ermakov
Yes, that part now works. Thank you, Wolfgang. Anton. On Wednesday, March 17, 2021 at 3:13:17 PM UTC-7 Wolfgang Bangerth wrote: > On 3/17/21 10:41 PM, Anton Ermakov wrote: > > Thanks for the reply. I think I got it. Basically, I had to manually > place the > > local matric

Re: [deal.II] treatment of handing nodes in the heat conduction problem

2021-03-17 Thread Anton Ermakov
ndices, mass_matrix);* Anton. On Tuesday, March 16, 2021 at 3:22:54 AM UTC-7 Wolfgang Bangerth wrote: > On 3/16/21 6:12 AM, Anton Ermakov wrote: > > > > I attach the minimally changed Step-26 code to illustrate the problem. > It is > > different from the original step

Re: [deal.II] project_boundary_values for components

2019-04-09 Thread Anton
Wolfgang, Thank you for the help. I think I have it working in the non-distributed case, but project_boundary_values segfaults on me without any explanations (compiled in debug mode) when run with a distributed triangulation. Is this "normal"? /Anton On Monday, April 8, 2019 a

Re: [deal.II] project_boundary_values for components

2019-04-06 Thread Anton
row number in the ConstraintMatrix - I cannot immediately figure this out from DoFHandler interface? /Anton On Saturday, April 6, 2019 at 5:11:09 PM UTC+2, Wolfgang Bangerth wrote: > > On 4/5/19 11:25 AM, Anton wrote: > > I need to apply a boundary condition (discontinuous function) to

[deal.II] project_boundary_values for components

2019-04-05 Thread Anton
s because the interpolation point falls exactly between the elements. Is there a way of applying project_boundary_values instead to fix this issue? Basically I cannot seem to find support for component selection in connection with project_boundary_values. Sincerely, /Anton -- The deal.II proje

[deal.II] DataPostprocessor for two DoFHandler s

2017-05-05 Thread Anton
one have to combine the two vectors into one by creating an FESystem with an associated DoFHandler first? Sincerely, /Anton -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received