Re: [deal.II] Augmented Lagrangian Method for Asymetric Stokes Problem
On 6/28/19 6:36 AM, Félix Bunel wrote: > > I tried a direct solver and it sucessfully solved my problem. That way I can > start working on the physical problem and I will see later for optimization. Glad to hear that works now. I think using a direct solver for first experiments is generally a good idea. Get things to work, and worry about making it fast at a later time. > As a thank you, here is a nice velocity field outputed by my code ! :-) Looks like a spiral galaxy :-) > I hope I will be able to finish this project and publish something with the > results. I will be sure to cite dealii if that happens again. (last project > with dealii is already published DOI:10.1103/PhysRevE.96.012705 ) Thanks in advance! (We already list your last project, I see.) Best W. -- Wolfgang Bangerth email: bange...@colostate.edu www: http://www.math.colostate.edu/~bangerth/ -- 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.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/1c0452a9-b5a1-3f38-cad1-48a29d666dcc%40colostate.edu. For more options, visit https://groups.google.com/d/optout.
[deal.II] How I see open source software
All, if you want to read a bit about how I got into open source software, and how I see my role in deal.II, here's something that just came out: https://bssw.io/blog_posts/leading-a-scientific-software-project-it-s-all-personal Cheers W. -- Wolfgang Bangerth email: bange...@colostate.edu www: http://www.math.colostate.edu/~bangerth/ -- 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.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/ec84578d-ac32-7f9e-9620-a25219c087e6%40colostate.edu. For more options, visit https://groups.google.com/d/optout.
[deal.II] Re: Storing additional information at the DOF unrelated to the system you are solving
Bruno, Hello everyone, I hope you are well. > I have a quick question which I cannot seem to wrap my hear around. I > think it will sound confused, but anyway. > Let's say I am solving for a scalar equation. I would have a > triangulation, then a DoFHandler and a FESystem with a multiplicity of one. > I can loop over the cells and assemble my equations. > If you only have a scalar equation, there is no need to use FESystem but you can do that of course. > Now, let's say I wanted to have an additional property that would be a > vector defined at the nodes of my triangulation using the positions of the > nodes. > However, this vector would not be used to be interpolated at the gauss > points, but would be use to cut the cells and do other geometric > manipulations. > I could define another DoFHandler , and FESystem, but then when I loop > over the cells, there would be nothing to make sure that my DOF on both > equations are numbered in the same way. > However, I cannot defined another vector using the FESystem of my scalar > equation, because the multiplicity is not the same. > What would be the best way to proceed? > I am not quite sure to understand what you are asking for. The support points only coincide for FE_Q(1) elements with the nodes of the triangulation. Normally, you would associate the entries of a Vector with the finite element vector defined by the degrees of freedom (and the ansatz space), but you can also store different information. You are talking about two equations. What is the second equation and how is it related to the first one? Can you give an example how you want to use the Vectors? Best, Daniel -- 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.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/c0f86511-a8d5-41df-8236-b7fc95ac2d00%40googlegroups.com. For more options, visit https://groups.google.com/d/optout.
[deal.II] Storing additional information at the DOF unrelated to the system you are solving
Hello everyone, I hope you are well. I have a quick question which I cannot seem to wrap my hear around. I think it will sound confused, but anyway. Let's say I am solving for a scalar equation. I would have a triangulation, then a DoFHandler and a FESystem with a multiplicity of one. I can loop over the cells and assemble my equations. Now, let's say I wanted to have an additional property that would be a vector defined at the nodes of my triangulation using the positions of the nodes. However, this vector would not be used to be interpolated at the gauss points, but would be use to cut the cells and do other geometric manipulations. I could define another DoFHandler , and FESystem, but then when I loop over the cells, there would be nothing to make sure that my DOF on both equations are numbered in the same way. However, I cannot defined another vector using the FESystem of my scalar equation, because the multiplicity is not the same. What would be the best way to proceed? TLDR: I would like to store additional information at the position of the DOFs, without including them in my vector of unknown, but I would still want to be able to know their position. Is this doable? Thanks! Bruno -- 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.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to dealii+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/8b4d0f2f-8889-4d97-bc3d-9d04c93c0167%40googlegroups.com. For more options, visit https://groups.google.com/d/optout.