[deal.II] Are there any tips for using Trilinos AMG preconditioner?

2025-01-28 Thread timo Hyvärinen
Hi, dealii devs and users I have been using the dealii Trilinos API of ML to generate an AMG preconditioner for MPI distributed sparse system matrix quite a while. It roubustely works for most computing of a damped Newton method on clusters. However, casually the linear solver of Newton update ve

Re: [deal.II] cmake configure does not find cgal

2024-01-05 Thread timo Hyvärinen
hi, praveen I think spack did cause something went wrong, building deal.ii with cgal without spack may be a good way. bset, Tim On Fri, Jan 5, 2024 at 4:25 PM Praveen C wrote: > So this error > > On 05-Jan-2024, at 7:28 PM, Praveen C wrote: > > [100%] Built target step-1 > $ ./step-1 > step-1

Re: [deal.II] cmake configure does not find cgal

2024-01-05 Thread timo Hyvärinen
ternal boost if I disable cgal. > > Thanks > praveen > > On 05-Jan-2024, at 3:11 PM, timo Hyvärinen > wrote: > > hi, praveen > > Did you check which Boost deal.ii was built from? The external one or > bundled one? If deal.ii didn't detect what "DBOOST_DI

Re: [deal.II] cmake configure does not find cgal

2024-01-05 Thread timo Hyvärinen
hi, praveen Did you check which Boost deal.ii was built from? The external one or bundled one? If deal.ii didn't detect what "DBOOST_DIR=`spack location -i boost` " is provided, then bundled Boost will be built. If this happens, according to " CGAL requires C++17 and an externally configured Boos

Re: [deal.II] cmake configure does not find cgal

2024-01-04 Thread timo Hyvärinen
hi, paraveen I think -DCGAL_DIR="/home/spack/opt/spack/linux-opensuse15-broadwell/gcc-12.3.0/cgal-5.5.2-6umjskepmayihkexoscjiovjmtgo5toy" will be better if spack location -i cgal return the library dictionary i.e., it contains include/, lib/. I'm not sure if shell evaluation (`` or $()) would be

[deal.II] cmake build linking error with source tree layout

2023-12-28 Thread timo Hyvärinen
Dear deal.ii developers and user community, Happy New Year! 🎊 I'm trying to split source code into different directories following to this documentation https://www.dealii.org/current/users/cmake_user.html Unfortunately, I keep getting the following error in linking step between libmylib.a and m

Re: [deal.II] Re: question about hybrid MPI and TBB thread parallelism

2023-10-02 Thread timo Hyvärinen
Bangerth wrote: > On 10/2/23 00:46, timo Hyvärinen wrote: > > (I) How to make the program run faster through Valgrind profiling? I > know I > > should use cachegrind and callgrind, but I don't know is what things I > should > > pay attention to from cachegrind/ca

Re: [deal.II] Re: question about hybrid MPI and TBB thread parallelism

2023-10-02 Thread timo Hyvärinen
sin wrote: > Tim, > > Valgrind is great but it is always slow. Instead you can use the > AddressSanitizer from clang or gcc. It's much faster than Valgrind but I > find the output is harder to read. > > Best, > > Bruno > > Le lun. 2 oct. 2023 à 02:47, timo Hyvär

Re: [deal.II] Re: what's the good way to use GPU acceleration in matrix-based code under dealii environment

2023-10-02 Thread timo Hyvärinen
Hi, Bruno, Thank you for your insightful and significant reply with the PR link. Assembly by GPU is still undergoing! This certainly is my most wished-for feature.😀 If currently the GPU is used for linear solver, then I think Kokkos leverage will be more natural for Trilinos ML library and Trili

[deal.II] what's the good way to use GPU acceleration in matrix-based code under dealii environment

2023-10-02 Thread timo Hyvärinen
Hi, deal community and developers, Thank you for all the help I ever got from you guys. I'm developing a MPI vector-valued non-linear PDEs solver. The way I chose is dumping newton iteration, which uses AMG preconditioner support from Trilinos 12.18 for linear solver. Currently, the program is CP

Re: [deal.II] Re: question about hybrid MPI and TBB thread parallelism

2023-10-01 Thread timo Hyvärinen
llo Wolfgang, > How are you doing today? > > Well your question is not quite clear. > So what exactly do you want to know? > > On Sun, Oct 1, 2023 at 3:40 PM Wolfgang Bangerth > wrote: > >> On 9/30/23 02:01, timo Hyvärinen wrote: >> > >> > So my q

Re: [deal.II] Re: question about hybrid MPI and TBB thread parallelism

2023-10-01 Thread timo Hyvärinen
: 7,738,278 allocs, 7,717,607 frees, 9,563,207,155 bytes allocated " I'm not sure if the mpi-valgrind suppress ever works, but there are about 10M leaks for one task. Tim, Sincerely On Mon, Oct 2, 2023 at 1:40 AM Wolfgang Bangerth wrote: > On 9/30/23 02:01, timo Hyvärinen wrote: >

Re: [deal.II] Re: question about hybrid MPI and TBB thread parallelism

2023-09-30 Thread timo Hyvärinen
want the debug > info. Without it, the information given by the profiler usually makes > little sense. So you compile in release mode but you use the following > option when configuring your deal.II "-DCMAKE_CXX_FLAGS=-g" > > Best, > > Bruno > > Le sam. 16 sept.

Re: [deal.II] Re: question about hybrid MPI and TBB thread parallelism

2023-09-16 Thread timo Hyvärinen
Hi Bruno, Thank you for your explanations. Seemingly, I should compile an optimized lib then do profiling. Sincerely, Timo On Fri, Sep 15, 2023 at 11:04 PM Bruno Turcksin wrote: > Timo, > > You will get vastly different results in debug and release modes for two > reasons. First, the compiler

Re: [deal.II] Re: question about hybrid MPI and TBB thread parallelism

2023-09-15 Thread timo Hyvärinen
hi, Marc, Thank you for the reply. I compiled the lib with debug mode, didn't try the optimized version. I didn't think this could be a significant issue, but I infer optimized lib could improve performance alot based on your question. Sincerely, Timo On Fri, Sep 15, 2023 at 8:21 PM Marc Fehlin

Re: [deal.II] question about hybrid MPI and TBB thread parallelism

2023-09-15 Thread timo Hyvärinen
e thread, but I believe it should be a significant experience. I will do profiling first and then talk in this thread later. Sincerely, Timo On Fri, Sep 15, 2023 at 3:51 PM Wolfgang Bangerth wrote: > On 9/15/23 03:17, timo Hyvärinen wrote: > > > > So here is my question, which tut

[deal.II] question about hybrid MPI and TBB thread parallelism

2023-09-15 Thread timo Hyvärinen
Dear dealii community and developers, I have used dealii framework (9.3.x) a while on HPC machine. My project involves solving vector-valued nonlinear PDE with nine components. Currently, I've implemented damping newton iteration with GMRES+AMG preconditioner with MPI on distributed memory archite