Re: [deal.II] Matrix Free 1D Adaptive Mesh Problem

2024-11-06 Thread Sean Johnson
fers from what you *expect* it to output? It's often difficult for > those who didn't write the code to say what is expected, and how what > you get differs from the expectation. > > Best > WB > -- The deal.II project is located at http://www.dealii.org/ For maili

Re: [deal.II] Matrix Free 1D Adaptive Mesh Problem

2024-11-06 Thread Wolfgang Bangerth
e code you attach outputs, and how that differs from what you *expect* it to output? It's often difficult for those who didn't write the code to say what is expected, and how what you get differs from the expectation. Best WB -- The deal.II project is located at http://www.dealii.

Re: [deal.II] Printing PETScWrappers::MPI::SparseMatrix to a single file

2024-11-06 Thread Wolfgang Bangerth
s own file (say, matrix-xxx.txt, where xxx is std::to_string(Utilities::MPI::this_mpi_process(MPI_COMM_WORLD)) ). Then, after you run is finished, do cat matrix-*.txt > matrix.txt to collate all of these files into one. Best WB -- The deal.II project is located at http://www.dealii.org/ For

[deal.II] Matrix Free 1D Adaptive Mesh Problem

2024-11-05 Thread Sean Johnson
re. The code is super simplistic and will not give you enough information about coordinates if you step up to 3D. Thanks for any help you can give, Sean Johnson -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/for

Re: [deal.II] Printing PETScWrappers::MPI::SparseMatrix to a single file

2024-11-05 Thread Yann Jobic
    }       system_matrix.compress(VectorOperation::add);       /// Write matrix       matrixfile << "\n System matrix values:" <*Is there any way to print the system_matrix to a single file either by using 'print()' or ‘write_ascii()’ function? Or do I need

[deal.II] Printing PETScWrappers::MPI::SparseMatrix to a single file

2024-11-05 Thread Sparsh Sinha
ompress(VectorOperation::add); /// Write matrix matrixfile << "\n System matrix values:" <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 su

[deal.II] deal.II Newsletter #302

2024-11-05 Thread 'Rene Gassmoeller&#x27; via deal.II User Group
Hello everyone! This is deal.II newsletter #302. It automatically reports recently merged features and discussions about the deal.II finite element library. ## Below you find a list of recently proposed or merged features: #17842: Add a test. (proposed by bangerth) https://github.com/dealii

Re: [deal.II] mesh

2024-11-04 Thread Luca Heltai
I read the file using TopoDS_Shape > > <https://dealii.org/developer/doxygen/deal.II/classTensor.html> bow_surface > > <https://dealii.org/developer/doxygen/deal.II/classTensor.html> = > > OpenCASCADE::read_SLT > > <https://dealii.org/developer/dox

Re: [deal.II] mesh

2024-11-04 Thread Wolfgang Bangerth
not *guess* what your problem is if you don't show us what it is you are trying." In particular, there is literally nothing we can suggest when you say "an error occurs" without telling us what the error is, where it comes from, and what you tried to address it. There are curr

Re: [deal.II] mesh

2024-11-04 Thread 王艺康
mber 4, 2024 at 10:54:35 AM UTC+8 Wolfgang Bangerth wrote: > >> On 11/2/24 21:11, 王艺康 wrote: >> > Yes I did, I read the file using TopoDS_Shape >> > <https://dealii.org/developer/doxygen/deal.II/classTensor.html> >> bow_surface < >> https://d

Re: [deal.II] mesh

2024-11-03 Thread Wolfgang Bangerth
On 11/2/24 21:11, 王艺康 wrote: Yes I did, I read the file using TopoDS_Shape <https://dealii.org/developer/doxygen/deal.II/classTensor.html> bow_surface <https://dealii.org/developer/doxygen/deal.II/classTensor.html> = OpenCASCADE::read_SLT <https://dealii.org/developer/

Re: [deal.II] mesh

2024-11-02 Thread 王艺康
Yes I did, I read the file using TopoDS_Shape <https://dealii.org/developer/doxygen/deal.II/classTensor.html> bow_surface <https://dealii.org/developer/doxygen/deal.II/classTensor.html> = OpenCASCADE::read_SLT <https://dealii.org/developer/doxygen/deal.II/namespaceO

Re: [deal.II] mesh

2024-11-02 Thread Wolfgang Bangerth
: 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

Re: [deal.II] mesh

2024-11-01 Thread 王艺康
- > Wolfgang Bangerth email: bang...@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.googl

Re: [deal.II] Introducing a Chatbot Assistant for deal.II Tutorials

2024-10-31 Thread Mehdi Memar
Dear deal.II Community, I am excited to introduce the second version of the deal.II Assistant, now covering 90 tutorials as well as the test suite available in the library’s source code. This update is possible thanks to Professor Heltai’s suggestion. As you may know, LLMs are available at

Re: [deal.II] Getting zero determinant for the stokes problem

2024-10-30 Thread Abbas Ballout
erion that can be tested numerically because floating point > numbers cannot distinguish whether the determinant is (0.1)^1000 or in > fact zero. > > Best > W. > -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.c

Re: [deal.II] Getting zero determinant for the stokes problem

2024-10-30 Thread Wolfgang Bangerth
rminant is (0.1)^1000 or in fact zero. Best W. -- 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&quo

Re: [deal.II] Getting zero determinant for the stokes problem

2024-10-30 Thread Abbas Ballout
of course also the issue that computing the determinant is an > exceptionally expensive operation (at least n^3). > > If you're interested in the smallest eigenvalue, compute that rather > than the determinant. > > Best > W. > -- The deal.II project is located at

Re: [deal.II] Missing follower force stiffness term in Step 44

2024-10-29 Thread Yiliang Wang
t; large-deformation modeling and not having written the program. But >> perhaps >> someone else knows an answer to your question. >> >> Best >> W. >> >> -- >> >> Wolfgang B

Re: [deal.II] Missing follower force stiffness term in Step 44

2024-10-29 Thread Jean-Paul Pelteret
t having written the program. But perhaps > someone else knows an answer to your question. > > Best > W. > > -- > > Wolfgang Bangerth email: bange...@colostate.edu > www: http://www.math.colost

Re: [deal.II] mesh

2024-10-28 Thread Wolfgang Bangerth
ostate.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

Re: [deal.II] Introducing a Chatbot Assistant for deal.II Tutorials

2024-10-28 Thread Mehdi Memar
Dear Bruno, Thank you for the kind words! I'm glad to hear the positive response. Looking forward to continuing our work together toward our shared goals. After all, we all love deal.II and its brothers! :) Best regards, Mehdi On Mon, Oct 28, 2024 at 7:19 PM blais...@gmail.com wrote: &

Re: [deal.II] Getting zero determinant for the stokes problem

2024-10-28 Thread Wolfgang Bangerth
the smallest eigenvalue, compute that rather than the determinant. Best W. -- 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 &qu

Re: [deal.II] Re: access coordinates of quadrature points

2024-10-28 Thread Wolfgang Bangerth
issue. My initial thought reading your question was that you're looking at the wrong statement and that the statement you pointed out probably did not cause the error. It seems like that was the right guess. Best W. -- The deal.II project is located at http://www.dealii.org/ For mai

[deal.II] Re: access coordinates of quadrature points

2024-10-28 Thread Yann Jobic
dim=2), i've got the error : An error occurred in line <2778> of file dealii/9.5.2_with_petsc/include/deal.II/base/table.h> in function     typename dealii::AlignedVector::reference dealii::Table<2, ElementType>::operator()(dealii::Table<2, ElementType>::size_type, dealii

[deal.II] access coordinates of quadrature points

2024-10-28 Thread Yann Jobic
#3 ./step-22: main I don't understand why !! How can i access to those coordinates ? Many thanks, Yann -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d

Re: [deal.II] Introducing a Chatbot Assistant for deal.II Tutorials

2024-10-28 Thread blais...@gmail.com
s may be a corpus that is too big for RAG, but it may be possible to >> do fine tuning with it, and it would be tremendously useful. >> >> Would you like to give this a shot? If this works out well, I think we >> could point to the chatbot from the main webpage of the library.

Re: [deal.II] Introducing a Chatbot Assistant for deal.II Tutorials

2024-10-26 Thread Mehdi Memar
f the library. > > L. > > > On 14 Oct 2024, at 13:34, Mehdi Memar wrote: > > > > Dear deal.II Community, > > > > I hope this message finds you well. I am excited to announce that I have > developed a chatbot assistant designed to assist users with questions >

[deal.II] deal.II Newsletter #301

2024-10-25 Thread 'Rene Gassmoeller&#x27; via deal.II User Group
Hello everyone! This is deal.II newsletter #301. It automatically reports recently merged features and discussions about the deal.II finite element library. ## Below you find a list of recently proposed or merged features: #17806: Chunking for WorkStream with Taskflow (proposed by RyanMoulday

Re: [deal.II] Introducing a Chatbot Assistant for deal.II Tutorials

2024-10-25 Thread Luca Heltai
useful. Would you like to give this a shot? If this works out well, I think we could point to the chatbot from the main webpage of the library. L. > On 14 Oct 2024, at 13:34, Mehdi Memar wrote: > > Dear deal.II Community, > > I hope this message finds you well. I am excited t

Re: [deal.II] mesh

2024-10-24 Thread Wolfgang Bangerth
and then use the cylindrical coordinate system to move the ideal mesh I would like to enquire about the feasibility of this As I already said, there is a function in deal.II that *does* the transformation, but you have to tell it *what the transformation is*. It is not enough to just provide

Re: [deal.II] mesh

2024-10-24 Thread 王艺康
model. If you >> have >> such a function, you can use >> GridTools::transform() >> to transform one mesh into the other. >> >> Best >> W. >> >> -- >> >> Wol

Re: [deal.II] mesh

2024-10-22 Thread 王艺康
ponding point in the heart model. If you have > such a function, you can use > GridTools::transform() > to transform one mesh into the other. > > Best > W. > > -- > > Wolfgang Bangerth email:

Re: [deal.II] Missing follower force stiffness term in Step 44

2024-10-22 Thread Wolfgang Bangerth
ostate.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 Group

Re: [deal.II] mesh

2024-10-22 Thread Wolfgang Bangerth
o the other. Best W. -- Wolfgang Bangerth email: bange...@colostate.edu www: http://www.math.colostate.edu/~bangerth/ -- The deal.II project is located at http://www.deali

[deal.II] Adding Transparency in VisIt Without Displaying Element Faces

2024-10-22 Thread Laryssa Abdala
Hi everyone, I have a question about postprocessing simulation data. Specifically, I’m running simulations with deal.ii and visualizing the VTU outputs in VisIt (LLNL). I’m trying to add transparency without displaying the element faces. In ParaView, this can be achieved using the "Clean to

[deal.II] Missing follower force stiffness term in Step 44

2024-10-22 Thread Yiliang Wang
total stiffness matrix in element level? If that term is ignored, I guess the convergence rate will be affected. I am not if the accuracy of the results will also be affected if that term is ignored. Best, Yiliang -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum op

Re: [deal.II] MPI with METIS causes memory leak

2024-10-22 Thread 'Alexander Greiner&#x27; via deal.II User Group
der Greiner' via deal.II User Group wrote: > > > > For step-17 I get the same error message indicating that the memory > > allocation failed, for step-18 and my personal code it just fills up the > > system RAM until one aborts it. (Everything works using only 1 core >

[deal.II] mesh

2024-10-22 Thread 王艺康
teful for any suggestions! -- 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 unsubs

Re: [deal.II] MPI with METIS causes memory leak

2024-10-21 Thread Wolfgang Bangerth
On 10/21/24 03:58, 'Alexander Greiner' via deal.II User Group wrote: For step-17 I get the same error message indicating that the memory allocation failed, for step-18 and my personal code it just fills up the system RAM until one aborts it. (Everything works using only 1 co

[deal.II] Postdoc stipend in computational mathematics

2024-10-21 Thread Siyang Wang
project, we aim to develop a novel and unconventional finite element method based on the Galerkin difference framework and resolve geometrical features by using the fictitious domain approach. We will implement the method in the state-of-the-art C++ software library deal.II. For details and

[deal.II] MPI with METIS causes memory leak

2024-10-21 Thread 'Alexander Greiner&#x27; via deal.II User Group
eally expect there is a solution for that, other than trying to recompile everything or switching to parallel::distributed::Triangulation. Still, I'm very grateful for any suggestions! All the best, Alex -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum op

Re: [deal.II] Re: cmake error when building with arbox

2024-10-20 Thread Marco Feder
> Primitives &, const SpaceFillingCurve &) [with > MemorySpace=Kokkos::HostSpace, BoundingVolume=ArborX::Box, Enable=void, > ExecutionSpace=Kokkos::DefaultHostExecutionSpace, > Primitives=std::vector, > std::allocator>>, > SpaceFillingCurve=ArborX::Experimental:

Re: [deal.II] Membrane problem and locking free method with deal.ii

2024-10-18 Thread Wolfgang Bangerth
On 10/18/24 03:09, Léonhard YU wrote:  I am working on a bending problem of membrane or thin plate, which is involved to the locking problem. I found a discussion in deal.ii gmail group in 2018 that Prof. Bangerth mentioned the mixed method which was available in deal.ii. I am wondering if

[deal.II] Membrane problem and locking free method with deal.ii

2024-10-18 Thread Léonhard YU
Hello everyone, I am working on a bending problem of membrane or thin plate, which is involved to the locking problem. I found a discussion in deal.ii gmail group in 2018 that Prof. Bangerth mentioned the mixed method which was available in deal.ii. I am wondering if there is any example for

Re: [deal.II] FENothing and `VectorTools::integrate_difference`

2024-10-17 Thread Davit Gyulamiryan
of > > the built-in functions, but I was interested if I have a bug somewhere > > in the code. > > > > Has anyone ever had this or similar problem, or know how to use the > > built-in funcionality of deal.ii for my case? > > It's of course hard to tell wi

Re: [deal.II] Re: cmake error when building with arbox

2024-10-17 Thread Kaushik Das
hy(const ExecutionSpace &, const Primitives &, const SpaceFillingCurve &) [with MemorySpace=Kokkos::HostSpace, BoundingVolume=ArborX::Box, Enable=void, ExecutionSpace=Kokkos::DefaultHostExecutionSpace, Primitives=std::vector, std::allocator>>, SpaceFillingCurve=ArborX::Experimental::Mor

Re: [deal.II] Re: cmake error when building with arbox

2024-10-17 Thread Kaushik Das
0:45:34 AM UTC-4 k.d...@gmail.com > wrote: > >>>>> > >>>>> Hi, > >>>>> If I turn on Arborx then cmake is also looking for ADOLC. > >>>>> Is Arborx support dependent on ADOLC? > >>>>> I have attached a log fr

Re: [deal.II] residual for a static condensation problem

2024-10-16 Thread Wolfgang Bangerth
d in data structures, but that they are conceptually different. You simply cannot multiply [Kpu] by a collection of velocity vectors; [Kpu] needs to be multiplied by a vector of nodal values. Best W. -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see

Re: [deal.II] residual for a static condensation problem

2024-10-16 Thread Yann Jobic
in step-57, but i do it on a simpler problem first (stokes). Thus if i define {R(u,p)} = {-2*div u + \nabla p -f}, I want to compute ||R(u,p)||. Yann, so the residual for you is a function of x, and what you're looking to compute is an integral (for the norm). Integrals in deal.II are al

Re: [deal.II] Re: cmake error when building with arbox

2024-10-16 Thread Daniel Arndt
;>> >>>> On Tuesday, October 15, 2024 at 10:45:34 AM UTC-4 k.d...@gmail.com wrote: >>>>> >>>>> Hi, >>>>> If I turn on Arborx then cmake is also looking for ADOLC. >>>>> Is Arborx support dependent on ADOLC? >>>&

Re: [deal.II] FENothing and `VectorTools::integrate_difference`

2024-10-16 Thread Wolfgang Bangerth
working so can just use it instead of the built-in functions, but I was interested if I have a bug somewhere in the code. Has anyone ever had this or similar problem, or know how to use the built-in funcionality of deal.ii for my case? It's of course hard to tell without seeing the com

Re: [deal.II] Re: cmake error when building with arbox

2024-10-16 Thread Kaushik Das
>>> >>>> Hi, >>>> If I turn on Arborx then cmake is also looking for ADOLC. >>>> Is Arborx support dependent on ADOLC? >>>> I have attached a log from a docker build. Arborx was built first then >>>> I used -DDEAL_II_WITH_ARBORX=ON and -

[deal.II] FENothing and `VectorTools::integrate_difference`

2024-10-16 Thread Davit Gyulamiryan
test with adaptive refinement of the boundary and the interior of the physical domain. Hanging nodes constraints are applied where needed. Error calculation using deal.ii `VectorTools::integrate_difference`: const ComponentSelectFunction velocity_mask( std::make_pair(0

Re: [deal.II] Re: cmake error when building with arbox

2024-10-16 Thread Bruno Turcksin
en cmake is also looking for ADOLC. >>> Is Arborx support dependent on ADOLC? >>> I have attached a log from a docker build. Arborx was built first then I >>> used -DDEAL_II_WITH_ARBORX=ON and -DARBORX_DIR="/usr/local/Arborx" but >>> cmake thorws an er

Re: [deal.II] Re: cmake error when building with arbox

2024-10-16 Thread Kaushik Das
om a docker build. Arborx was built first then I >> used -DDEAL_II_WITH_ARBORX=ON and -DARBORX_DIR="/usr/local/Arborx" but >> cmake thorws an error "ADOLC_INCLUDE_DIR not found!" >> >> Thanks, >> Kaushik >> > -- > The deal.II project is l

[deal.II] Re: cmake error when building with arbox

2024-10-16 Thread Bruno Turcksin
d -DARBORX_DIR="/usr/local/Arborx" but > cmake thorws an error "ADOLC_INCLUDE_DIR not found!" > > Thanks, > Kaushik > -- 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 --

[deal.II] deal.II Newsletter #300

2024-10-16 Thread 'Rene Gassmoeller&#x27; via deal.II User Group
Hello everyone! This is deal.II newsletter #300. It automatically reports recently merged features and discussions about the deal.II finite element library. ## Below you find a list of recently proposed or merged features: #17785: Fix cubic simplex constraints. (proposed by drwells) https

Re: [deal.II] Import a mesh with generated grain structure

2024-10-15 Thread Kirill Romanov
Thanks for the answers! I've been thinking about solving my question along roughly the same lines. вт, 15 окт. 2024 г. в 18:11, Praveen C : > I think you can get this with > > https://dealii.org/current/doxygen/deal.II/structCellData.html > > > https://dealii.org/

Re: [deal.II] residual for a static condensation problem

2024-10-15 Thread Wolfgang Bangerth
want to compute ||R(u,p)||. Yann, so the residual for you is a function of x, and what you're looking to compute is an integral (for the norm). Integrals in deal.II are always computed using FEValues, and what you want here is going to look something like the following: double integra

Re: [deal.II] Error GridIn::read_msh(), mesh with multiple transfinite regions

2024-10-15 Thread Wolfgang Bangerth
-- 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 --

Re: [deal.II] Error GridIn::read_msh(), mesh with multiple transfinite regions

2024-10-15 Thread Matteo Malvestiti
sh. > > best > praveen > > -- 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" gro

Re: [deal.II] Error GridIn::read_msh(), mesh with multiple transfinite regions

2024-10-15 Thread Praveen C
The original mesh did have ids on interior edges. You need to set physical tags to all boundary edges and to cells also. Use attached filed and make sure to actually regenerate the mesh. best praveen -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options

Re: [deal.II] Error GridIn::read_msh(), mesh with multiple transfinite regions

2024-10-15 Thread Wolfgang Bangerth
. -- 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

Re: [deal.II] Error GridIn::read_msh(), mesh with multiple transfinite regions

2024-10-15 Thread Matteo Malvestiti
dd this and try > > Physical Curve("others", 3) = {2,4,5,6}; > > best > praveen > > On 15 Oct 2024, at 2:43 PM, Matteo Malvestiti wrote: > > [Errata corrige] > Hello everyone. > I'm getting an error when importing msh files in deal.II with >

[deal.II] cmake error when building with arbox

2024-10-15 Thread Kaushik Das
INCLUDE_DIR not found!" Thanks, Kaushik -- 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" g

Re: [deal.II] Error GridIn::read_msh(), mesh with multiple transfinite regions

2024-10-15 Thread Praveen C
It is because you have not set ids for all boundary faces. Add this and try Physical Curve("others", 3) = {2,4,5,6}; best praveen > On 15 Oct 2024, at 2:43 PM, Matteo Malvestiti wrote: > > [Errata corrige] > Hello everyone. > I'm getting an error when impor

Re: [deal.II] Import a mesh with generated grain structure

2024-10-15 Thread Praveen C
I think you can get this with https://dealii.org/current/doxygen/deal.II/structCellData.html https://dealii.org/current/doxygen/deal.II/classCellAccessor.html#ae4769702cd7ab67a61b25778ea3021b2 cell->material_id() This should give the physical tag you have set to each cell in your gmsh g

[deal.II] Error GridIn::read_msh(), mesh with multiple transfinite regions

2024-10-15 Thread Matteo Malvestiti
[Errata corrige] Hello everyone. I'm getting an error when importing msh files in deal.II with GridIn::read_msh(). For the sake of simplicity I brought a very basic grid: - Rectangle 2x1 - Divided in two squares 1x1, both uniformly meshed, transfinite, with different granularity When I i

[deal.II] Error GridIn::read_msh(), simple mesh with two transfinite regions

2024-10-15 Thread Matteo Malvestiti
When I try to import it in deal.II with: const std::string filename = "mesh.msh"; std::ifstream input_file(filename); Triangulation<2> triangulation; GridIn<2> grid_in; grid_in.attach_triangulation(triangulation); grid_in.read_msh(input_file); I get this error: termina

Re: [deal.II] residual for a static condensation problem

2024-10-15 Thread Yann Jobic
be ordered as [Kup], [Kpu], ... Yann -- the term "residual" is used in different meanings in numerical analysis. Can you define precisely what it is you want to compute? Best  W. -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https:/

Re: [deal.II] residual for a static condensation problem

2024-10-14 Thread Wolfgang Bangerth
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 G

Re: [deal.II] Import a mesh with generated grain structure

2024-10-14 Thread Wolfgang Bangerth
On 10/14/24 07:55, Kirill Romanov wrote: ** We are developing crystal plasticity material models in a team and we want to embed our models in deal.ii. So I have a question – are there any simple tools to import a mesh (we use a simplex mesh) with generated grain structure into deal.ii (e.g

[deal.II] residual for a static condensation problem

2024-10-14 Thread Yann Jobic
ambda the penalisation parameter. Thus the weak form looks like (v shape function of u, and phi shape function of p), 2*a(grad v, grad u)-a(div v,p)-a(phi,div u)+a(phi,\lambda p) = l(v*f) then {P}=-[Kpp]^{-1}[Kpu]*{U} and finally ([Kuu]-[Kup][Kpp]^-1[Kpu]){U}={f} -- The deal.II project is located at

[deal.II] Import a mesh with generated grain structure

2024-10-14 Thread Kirill Romanov
We are developing crystal plasticity material models in a team and we want to embed our models in deal.ii. So I have a question – are there any simple tools to import a mesh (we use a simplex mesh) with generated grain structure into deal.ii (e.g. from Neper)? I have not been able to find a

Re: [deal.II] Parallization Step-40

2024-10-14 Thread James Shi
:47写道: > Dear Prof. Timo Heister and Wolfgang Bangerth > > I am trying to solve my multiphysics problem by deal.ii and it's > paralization. > > I started from STEP-40 and encountered a problem regarding paralization > and MPI. > > The program didn't pro

[deal.II] Introducing a Chatbot Assistant for deal.II Tutorials

2024-10-14 Thread Mehdi Memar
Dear deal.II Community, I hope this message finds you well. I am excited to announce that I have developed a chatbot assistant designed to assist users with questions specifically related to the tutorials available in the deal.II documentation. This assistant sources information from all 90

Re: [deal.II] Parallization Step-40

2024-10-10 Thread Wolfgang Bangerth
ther processes that are already done assembling on their cells. You end up with a deadlock in these situations, and that's what you observe. There is a little bit more information about "Collective operations" in the glossary: https://dealii.org/developer/doxygen/deal.II/DEALG

[deal.II] Parallization Step-40

2024-10-10 Thread Junxiang Wang
Dear Prof. Timo Heister and Wolfgang Bangerth I am trying to solve my multiphysics problem by deal.ii and it's paralization. I started from STEP-40 and encountered a problem regarding paralization and MPI. The program didn't proceed after I made a subrout

Re: [deal.II] Re: FSI with separate grids.

2024-10-08 Thread Luca Heltai
-- > 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 ma

Re: [deal.II] Re: FSI with separate grids.

2024-10-08 Thread Wolfgang Bangerth
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/ F

[deal.II] Re: FSI with separate grids.

2024-10-08 Thread Tim Braun
es or perhaps using a > connector to map one solution to the other? I suppose I could resort to > preCICE but imagine there's a way to run this within a deal ii "program". > > I hope my questions are somewhat clear enough, because I'm learning as I > go. Any fi

Re: [deal.II] Should cubic FE_SimplexP work in an FESystem on master?

2024-10-08 Thread Kyle Schwiebert
ail: 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 > --- &

[deal.II] FSI with separate grids.

2024-10-08 Thread Tim Braun
ps using a connector to map one solution to the other? I suppose I could resort to preCICE but imagine there's a way to run this within a deal ii "program". I hope my questions are somewhat clear enough, because I'm learning as I go. Any first thoughts are appreciated -- The

Re: [deal.II] Should cubic FE_SimplexP work in an FESystem on master?

2024-10-08 Thread Wolfgang Bangerth
/~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 g

[deal.II] Should cubic FE_SimplexP work in an FESystem on master?

2024-10-07 Thread Kyle Schwiebert
rstand why this is happening. Unfortunately, I wouldn't have any comment on which specific commit might have caused the issue (if it is an issue) as I just switched to my own build of deal.II and was previously using 9.5.1. I'm not sure if this new instalation could be the issue, so

Re: [deal.II] Suggestions for preconditioners of implicit time discretized wave equations

2024-10-04 Thread 'yy.wayne&#x27; via deal.II User Group
Thank you Prof. Bangerth for your time and very insightful suggestions. I'll definitely do more computational experiments and check literatures . Best, Wayne 在2024年10月5日星期六 UTC+8 11:49:54 写道: > On 10/4/24 21:44, 'yy.wayne' via deal.II User Group wrote: > > ** > >

Re: [deal.II] Suggestions for preconditioners of implicit time discretized wave equations

2024-10-04 Thread Wolfgang Bangerth
On 10/4/24 21:44, 'yy.wayne' via deal.II User Group wrote: ** I decrease dt by half for each mesh refinement. Now the iterations remain constant (20, 21, and 23), same as you predicted. I guess in practical I should keep dt and h constant (based on frequency and wave speed) for a

Re: [deal.II] Suggestions for preconditioners of implicit time discretized wave equations

2024-10-04 Thread 'yy.wayne&#x27; via deal.II User Group
test the number of iterations? 在2024年10月5日星期六 UTC+8 11:33:36 写道: > I didn't know that dt should be propotional to h/c for accuracy > results. Thank you for > pointing it out. > > 在2024年10月5日星期六 UTC+8 11:15:07 写道: > >> On 10/4/24 20:39, 'yy.wayne' via de

Re: [deal.II] Suggestions for preconditioners of implicit time discretized wave equations

2024-10-04 Thread 'yy.wayne&#x27; via deal.II User Group
I didn't know that dt should be propotional to h/c for accuracy results. Thank you for pointing it out. 在2024年10月5日星期六 UTC+8 11:15:07 写道: > On 10/4/24 20:39, 'yy.wayne' via deal.II User Group wrote: > > > > I'm solving a time-dependent wave equ

Re: [deal.II] Suggestions for preconditioners of implicit time discretized wave equations

2024-10-04 Thread Wolfgang Bangerth
On 10/4/24 20:39, 'yy.wayne' via deal.II User Group wrote: I'm solving a time-dependent wave equation, and dt is the time step. For backward Euler time integration, I get M+dt^2*K = (v,u) + (\delta T)^2*(\nabla v, \nabla u). The time step dt does not scale with mesh size h, s

Re: [deal.II] Suggestions for preconditioners of implicit time discretized wave equations

2024-10-04 Thread 'yy.wayne&#x27; via deal.II User Group
M + dt K > which will (if you choose dt=h) have a condition number that grows like > h^{-1}. > > As a consequence, I'm a *bit* surprised that your iterations grow by so > much, but it's been a long time since I've solved the wave equation and > so I don't reca

Re: [deal.II] Suggestions for preconditioners of implicit time discretized wave equations

2024-10-04 Thread Wolfgang Bangerth
7;m a *bit* surprised that your iterations grow by so much, but it's been a long time since I've solved the wave equation and so I don't recall whether what you see is expected or not. Best W. -- The deal.II project is located at http://www.dealii.org/ For mailing list/

Re: [deal.II] Suggestions for preconditioners of implicit time discretized wave equations

2024-10-04 Thread 'yy.wayne&#x27; via deal.II User Group
yy.wayne' via deal.II User Group wrote: > > > > I've considered multigrid preconditioners as the iterations are > relatively > > constant, but isn't it suits diffusion problems more, > > as mentioned here (A former discussion > > < > http

Re: [deal.II] Suggestions for preconditioners of implicit time discretized wave equations

2024-10-04 Thread Wolfgang Bangerth
On 10/4/24 08:10, 'yy.wayne' via deal.II User Group wrote: I've considered multigrid preconditioners as the iterations are relatively constant, but isn't it suits diffusion problems more, as mentioned here (A former discussion <https://groups.google.com/g/dealii/c/det

Re: [deal.II] Suggestions for preconditioners of implicit time discretized wave equations

2024-10-04 Thread 'yy.wayne&#x27; via deal.II User Group
s expensive for such 'easy' problem(when dt is small). The only drawback for jacobi preconditioner is the iteration number scales with problem size ... 在2024年10月4日星期五 UTC+8 22:02:25 写道: > On 10/4/24 01:38, 'yy.wayne' via deal.II User Group wrote: > > ** > > &g

Re: [deal.II] Suggestions for preconditioners of implicit time discretized wave equations

2024-10-04 Thread Wolfgang Bangerth
On 10/4/24 01:38, 'yy.wayne' via deal.II User Group wrote: ** I tried by SSOR gives same behavior as jacobi. For mesh with 0, 1, and 2 times of refinement, it needs (jacobi: 20, 31, 67) and (SSOR: 10, 15, 31) iterations. Besides, is SSOR not applicable in matrix-free framework?

Re: [deal.II] Problems during the installation of an older version of deal.II (version 8.4)

2024-10-04 Thread Mahdi Gharehbaygloo
know it is not exact, but it conveys the message. 😅 From: dealii@googlegroups.com on behalf of Nicolò Lontani Sent: Thursday, October 3, 2024 9:57:04 AM To: deal.II User Group Subject: [deal.II] Problems during the installation of an older version of deal.II

Re: [deal.II] Suggestions for preconditioners of implicit time discretized wave equations

2024-10-04 Thread 'yy.wayne&#x27; via deal.II User Group
I tried by SSOR gives same behavior as jacobi. For mesh with 0, 1, and 2 times of refinement, it needs (jacobi: 20, 31, 67) and (SSOR: 10, 15, 31) iterations. Besides, is SSOR not applicable in matrix-free framework? 在2024年10月4日星期五 UTC+8 11:54:41 写道: > On 10/3/24 01:01, 'yy.wayne'

Re: [deal.II] run the step-1 tutorial in Visual Studio 2019

2024-10-03 Thread Wolfgang Bangerth
late already defined step1 C:\dealii-9.6.0\deal.II\bundled\kokkos-3.7.00\tpls\desul\include\desul\atomics\Generic.hpp 326 Error C2995: “T desul::Impl::atomic_fetch_oper(const Oper &,T *const ,dont_deduce_this_parameter::type,MemoryOrder,MemoryScope)”: function template already defined step1 C

  1   2   3   4   5   6   7   8   9   10   >