Re: [deal.II] Re: Assemble function, long time

Hermes, For large systems, Krylov solvers are faster and require less memory than direct solvers. Direct solvers scale poorly, in terms of memory and performance, with the number of unknowns. The only problem with Krylov solvers is that you need to use a good preconditioner. The choice of the prec

Re: [deal.II] University cluster installation issue: compilers

Stephanie, If you are still using the same compiler as before: icc version 19.0.4.243 (gcc version 4.8.5 compatibility) The problem is because of the standard library. The Intel compiler does not have its own standard library, instead it uses the library from gcc. In your case the compiler is i

Re: [deal.II] Re: Assemble function, long time

Hermes, I think Cuthill-McKee only works on symmetric matrices, is your matrix symmetric? Also, the goal of Cuthill-McKee is to help with the fill in of the matrix.There is no guarantee that it helps with the performance. If you don't know which preconditioner to use, you can use ILU (Incomplete L

Re: [deal.II] Re: Assemble function, long time

If your matrix is symmetric definite positive, you use CG . Otherwise, you use GMRES . Here is the page for ILU

Re: [deal.II] Re: Assemble function, long time

Yes, you should use your system_matrix. AdditionalData can be used to modify the parameters used by ILU. The interface of PreconditionILU should work very similarly to BlockJacobi see https://dealii.org/current/doxygen/deal.II/step_17.html#ElasticProblemsolve There are several tutorials that use p

Re: [deal.II] Re: Assemble function, long time

Hermes, Sorry, I don't use petsc. Maybe someone else can help you. Best, Bruno Le lun. 14 mars 2022 à 05:42, Hermes Sampedro a écrit : > Dear Bruno, > > I have been reading the examples and documents you pointed out. I tried to > use SolvereGREMS with PreconditionILU. However, I am getting a

[deal.II] Re: Error while installing dealii

Hello, Did you configure Trilinos with MPI support? Best, Bruno On Tuesday, March 15, 2022 at 7:09:55 AM UTC-4 punith...@gmail.com wrote: > I've compiled a deal.ii using CMake with the below options > > cmake -DCMAKE_INSTALL_PREFIX=/apps/codes/deal_II/9.3.3 > -DHDF5_DIR=/apps/libs/hdf5/gcc/1

[deal.II] Re: Solving separate FEM problems on subdomains

Corbin, Creating multiple Triangulations and mapping them is not easy to do. Instead, you want to use FE_Nothing to ignore some cells. Take a loot at step-46 https://dealii.org/current/doxygen/deal.II/step_46.html Best, Bruno On Monday, March 21, 2022 at 6:47:38 PM UTC-4 corbin@gmail.com

[deal.II] Re: range-based loop over cells

Niklas, If you want cell_iterators, you write for(auto cell : triangulation.cell_iterators ()) If you want active_cell_iterators, you write for(auto cell : triangulation.active_cell_iterators ()) The documentation is here

Re: [deal.II] Re: range-based loop over cells

Niklas, You can take a look here but basically if a class defines the begin() and end() iterators, you can do a range-based loop. That's just the C++11 standard. Best, Bruno Le ven. 25 mars 2022 à 08:21, Niklas Fehn a écrit : > so this is

[deal.II] Re: Testing: different results on local machine and docker container

Sebastian, You should get the same result with Docker and on your machine. I have had a few cases where the result was wrong in docker but correct on my machine. This was always related to a variable that was not initialized. For some reason, on my machine the variable would be set to zero but

Re: [deal.II] Supported matrix/vector to LinearOperator

Chen, Yes, LinearOperator does support Trilinos and PETSc matrices. LinearOperator even supports your own matrix type as long as you define vmult and Tvmult (see here ) Best, Bruno On Mond

Re: [deal.II] Supported matrix/vector to LinearOperator

> there any solution? > > best > Chen > > 在2022年4月18日星期一 UTC+8 22:07:36 写道： > >> Bruno Turcksin >> Thank you! it really helpful! >> >> best >> chen >> 在2022年4月18日星期一 UTC+8 21:07:45 写道： >> >>> Chen, >>> >>&

[deal.II] Re: Python bindings

Corbin, On Tuesday, April 19, 2022 at 3:04:39 PM UTC-4 corbin@gmail.com wrote: > I've built deal.ii version 9.3.0 with the > -DDEAL_II_COMPONENT_PYTHON_BINDING=ON configuration, yet when I run the > notebook, the Jupyter kernel is unable to locate the PyDealII module, that > is, the notebo

Re: [deal.II] Re: Python bindings

Corbin, That's great. Let us know if you encounter any other problems. What kind of bindings do you want to write? BTW, do you want to open a PR with the FindPython3 change? Best, Bruno Le mer. 20 avr. 2022 à 06:56, Corbin Foucart a écrit : > It looks like in the python-bindings CMakeLists.tx

Re: [deal.II] Re: Python bindings

Corbin, Le mer. 20 avr. 2022 à 19:10, Corbin Foucart a écrit : > > I'm hesitant to propose a change, as doing so may break the build process > for other people using different Python binding configurations. I'm also by > no means a CMake expert :-) However, if it's a safe assumption that anyone

[deal.II] Re: Help with modifying step 20; BlockMatrix structure

Ali, I am not sure I understand your question. Do you want to know how the local matrix, which is a FullMatrix, is correctly distributed in system_matrix, which is a BlockMatrix? This information is stored in the DoFHandler. You access this information by calling cell->get_dof_indices() Best,

[deal.II] Re: Setting Element Order Independently

Shane, Yes, it is possible. There are two different functions depending on what you want to do. If you want to set the order during setup you can use this function It is being used in step

[deal.II] Re: Periodic boundary conditions with AMR

Hi, There could be many reasons why your code is not working as expected. To find the problem, try to make the simplify the code as much as possible. For example, does the code fails if you don't use periodic boundaries? If it fails, try to simplify even more, for example, don't use periodic b

[deal.II] Re: SCALAPACK symbol check for pdsyevr_ and pssyevr_ failed

Luca, The problem is that your version of scalapack is missing functions that we need (see: https://github.com/dealii/dealii/issues/7774) We don't know why some versions of scalapack don't have these symbols. You will probably need to compile scalapack yourself. Best, Bruno On Thursday, June

[deal.II] Re: step-77 in parallel - SUNDIALS::KINSOL

Jose, What's the error that you get? Is there a message? Is it a segfault? Best, Bruno On Wednesday, June 8, 2022 at 1:03:06 PM UTC-4 jose.a...@gmail.com wrote: > Dear dealii community, > > I am solving a nonlinear problem that has a regularized sign function. The > classic Newton-Raphson met

Re: [deal.II] Can I use 3D prism elements (extruded triangles) in deal.ii?

Amit, I think that GridGenerator::quarter_hyper_shell is what you want. Best, Bruno On Friday, June 24, 2022 at 9:02:11 AM UTC-4 Paras Kumar wrote: > HI Amit, > > Focusing on the questi

[deal.II] Re: WorkStream::run error for hp iterator

Elyn, hp support has been added a couple of month ago. This means that you need to update deal.ii to 9.4.0 which has just been released. Best, Bruno On Monday, June 27, 2022 at 1:11:09 AM UTC-4 Yilin Luo wrote: > Hi, > > I'm trying to implement multithreading for a fluid-solid interaction >

Re: [deal.II] python binding

Chen, Yes, that should be it. It doesn't work? Best, Bruno On Tuesday, June 28, 2022 at 3:53:11 AM UTC-4 hkch...@gmail.com wrote: > Dear all, > > I compiled deal.ii with python binding. and I want to run the example of > step-49.ipynb. I have installed jupyter by pip3 install jupyter. Now, my

Re: [deal.II] Re: step-77 in parallel - SUNDIALS::KINSOL

Jose, First, I have never used KINSOL so I don't know what it requires. With that being said, this x.reinit(locally_relevant_dofs, MPI_COMM_WORLD); is wrong. Either you only need the locally_owned_dofs and you write x.reinit (locally_owned_dofs, MPI_COMM_WORLD); or you need the locally_relevant_d

Re: [deal.II] Re: step-77 in parallel - SUNDIALS::KINSOL

Jose, Sorry, you can give reinit a single argument to create ghosted vector when using a Trilinos vector. I usually use deal.II's own vector and they behave slightly differently. The problem when you give a single argument is that there is no way to know what is locally owned and thus you canno

[deal.II] Re: error while installing dealii using candi

Hello, This error often happens when you don't have enough memory to compile the program. Try using only one or two cores. Best, Bruno On Wednesday, July 6, 2022 at 4:01:45 AM UTC-4 me20...@iittp.ac.in wrote: > Dear all, > please find the error attached. > > thank you in advance.[image: PHOTO

[deal.II] Re: spack load dealii@9.3.3 show error: list index is out of range

Hello, This is a spack problem. You should open an issue on their github page. Best, Bruno On Friday, July 8, 2022 at 11:08:57 AM UTC-4 sj20...@gmail.com wrote: > *Steps to reproduce* > > I installed dea...@9.3.3 with command spack install -v > dealii~examples~threads. After installation, I l

[deal.II] Re: error installing dealii

Hello, It looks like you installed deal.II in /home/ajbar/Bureau/dealii How did you configure it? To get the code to compile, you need to reconfigure, recompile, and reinstall deal.II This time when you configure deal.II, make sure that the configuration line contains -DDEAL_II_WITH_MPI=ON and

Re: [deal.II] Re: error installing dealii

Mahmoud, You should add -DDEAL_II_WITH_MPI=ON and -DDEAL_II_WITH_TRILINOS=ON to the cmake command. Then, you can compile and install deal.II You should not modify deal.IIConfig.cmake (or any other file for that matter). Th cmake like should look like: cmake -DDEAL_II_WITH_MPI_MPI=ON -DDEAL_II_W

Re: [deal.II] Re: error installing dealii

Mahmoud, Make sure you write the cmake command from an empty directory. Otherwise you can sometimes get weird results. If that doesn't help, please send the CMakeOutput.log and the CMakeError.log files. Best, Bruno Le jeu. 28 juil. 2022 à 04:32, Mahmoud Ajbar a écrit : > Hello Bruno, > > Tha

Re: [deal.II] Re: error installing dealii

It looks like there is a problem with MUMPS dependencies. The compiler complains that it cannot find /usr/lib/x86_64-linux-gnu/libptscotch.so Can you check that this file exists. If it doesn't, try sudo apt install libmumps-ptscotch-dev or sudo apt install libptscotch-dev. Bruno Le jeu. 28 juil

[deal.II] Re: interface between dealii and matlab for data transfer (using mex files?)

Simon, You could just write/read files to transfer the data and you can launch your deal.II code from Matlab using https://www.mathworks.com/help/matlab/ref/system.html This way Matlab does not need to know anything about deal.II. You just write and read files, and you execute a bash command.

Re: [deal.II] Re: interface between dealii and matlab for data transfer (using mex files?)

Simon Le lun. 15 août 2022 à 08:52, Simon Wiesheier a écrit : > > I also thought about transferring the data by reading from/writing to a > file. > But I was looking for a more general approach without having to create > files. > > Writing files is no less general than coupling directly Matlab a

[deal.II] Re: Iterating over mesh cells in a custom order

Corbin, It's possible to do it using WorkStream::run (see here ) However, you need to create the ordering manually by "coloring" the cells. All the cells in the same color can be worked on

[deal.II] Re: Solving the linear system of equations using PETSc BlockSparseMatrix

Hi, If you search for "block solver" here https://dealii.org/developer/doxygen/deal.II/Tutorial.html, you will see all the tutorials that use block solvers. I think that only deal.II's own solvers support BlockSparseMatrix directly. Best, Bruno On Monday, August 22, 2022 at 9:02:28 AM UTC-4

Re: [deal.II] Re: Solving the linear system of equations using PETSc BlockSparseMatrix

They do. Take for example step-70, the matrix is a LA::MPI::BlockSparseMatrix which is a PETSc or a Trilinos BlockSpaceMatrix depending on how deal.II was configured. Bruno Le mar. 2

Re: [deal.II] Re: Solving the linear system of equations using PETSc BlockSparseMatrix

I am not sure what you mean by solving for diagonal blocks. The solver is general; it is the preconditioner that is block diagonal. If you want to use MPI, I assume that you have a large system. In this case, direct solvers are not very effective and you need to use an iterative solver. deal.II's i

[deal.II] Re: Add/remove some cells "one by one"

Olivier, I don't think what you want to do is possible. deal.II delegates the load balancing and the re-distribution to p4est. My understanding is that p4est cannot do what you ask and therefore deal.II can't do it either. Best, Bruno On Friday, August 26, 2022 at 8:36:13 AM UTC-4 olivier...@

[deal.II] Re: Python Binding

Felix, First, you need to configure deal.II using -DDEAL_II_COMPONENT_PYTHON_BINDINGS=ON. Then you need to re-compile and re-install the library. Now the wrappers are built but you still need to tell python where to find them. You can do that by setting the PYTHONPATH to the directory where py

[deal.II] Re: Trying to increase precision of floating point values in step3 of tutorial

Hello, What do you mean exactly by "deal.ii outputs the solution vector values upto 7 decimal"? Which function are you using? If you are using a function like std::cout to inspect the values in the vector, the default precision is 6 (I think) but it can be changed using std::setprecision (see

Re: [deal.II] Re: Python Binding

Felix, Glad to hear it is working and thank you for letting us know. Best, Bruno Le mer. 5 oct. 2022 à 11:17, Felix Felix a écrit : > Hi > > Thank you Bruno for your help I have resolved the issue, it turns out I > had many problems throughout the installations with my directories. > I reinst

Re: [deal.II] Re: Trying to increase precision of floating point values in step3 of tutorial

You want to do something like this: DataOut<2> data_out; data_out.attach_dof_handler(dof_handler); data_out.add_data_vector(solution, "solution"); data_out.build_patches(); std::ofstream output("solution.vtk"); output.set_precisition(12); data_out.write_vtk(output); Best, Bruno L

Re: [deal.II] Re: Trying to increase precision of floating point values in step3 of tutorial

Sorry, it should be output.precision(12) Bruno Le mar. 11 oct. 2022 à 09:03, Wasim Niyaz Munshi ce21d400 < ce21d...@smail.iitm.ac.in> a écrit : > Did you mean > output.set_precision(12); > > It gives the following error > error: ‘std::ofstream’ {aka ‘class std::basic_ofstream’} has no > member

[deal.II] Re: MUMPS solver with Trilinos

Rahul, On Sunday, October 16, 2022 at 11:42:50 AM UTC-4 grrah...@gmail.com wrote: > > 1. Does mumps has to be any specific version to compile with amesos. > Trilinos documentation available on net suggest 4.7.3. However, this is > fairly old version. > I don't know the answer to this question

[deal.II] Re: measuring cpu and wall time for assembly routine

Simon, The best way to profile a code is to use a profiler. It can give a lot more information than what simple timers can do. You say that your code is not parallelized but by default deal.II is multithreaded . Did you set DEAL_II_NUM_THREADS=1? That could explain why CPU and Wall time are di

Re: [deal.II] Re: measuring cpu and wall time for assembly routine

Simon, Le mer. 19 oct. 2022 à 09:33, Simon Wiesheier a écrit : > Thank you for your answer! > > " Did you set DEAL_II_NUM_THREADS=1?" > > How can I double-check that? > ccmake . > only shows my the variables CMAKE_BUILD_TYPE and deal.II_DIR . > But I do do knot if this is the right place to loo

[deal.II] Re: MacOS Installation - Monterey - M1 - spack

Francesco, You can do spack spec dealii%clang to see what spack will install. You can then check which package requires OpenGL (it is not deal.II) and then remove support for this package if you don't need it. For instance, let's say that it's trilinos that's pulling OpenGL and you don't nee

[deal.II] Re: Solving elasticity and laplace(damage) equation, one after the other.

Wasim, It's hard to say without seeing any code. Are you working with one or two DoFHandler? If you are working with two DoFHandler, are you sure that you are using the correct one? Best, Bruno On Monday, November 21, 2022 at 5:54:22 AM UTC-5 ce21...@smail.iitm.ac.in wrote: > Hello everyone

Re: [deal.II] Re: Solving elasticity and laplace(damage) equation, one after the other.

Wasim, Before the line that errors out can you do std::cout<get_fe().get_name()< a écrit : > I am working with 2 DoFHandlers. I checked my code. I am using the correct > DoFHandlers for the 2 equations. > > On Mon, Nov 21, 2022 at 6:58 PM Bruno Turcksin > wrote: > >

Re: [deal.II] Re: Solving elasticity and laplace(damage) equation, one after the other.

g to me, comes in this line in H function : > fe_values_elastic.reinit(cell);* > I need to *reinit *the fe_values_elasticity for each cell, but the > problem is that the *cell *I am passing as an argument lives on the > traingulation_damage while fe_values_elasticity lives on > triangulation_elasticity

Re: [deal.II] Re: Solving elasticity and laplace(damage) equation, one after the other.

all boundaries might have same bcs but for elasticity different bcs > exist at different faces of the square domain). > > On Tue, Nov 22, 2022 at 2:25 AM Bruno Turcksin > wrote: > >> Wasim, >> >> Why do you have two different triangulations? Can you have single >

[deal.II] Re: How to solve 2d elasticity problem for two different materials.

Deepika, You can take a look at step-28. We use material_id to represent different materials. Best, Bruno On Monday, November 28, 2022 at 1:56:22 AM UTC-5 deepika...@iitgoa.ac.in wrote: > Hello Everyone, > > I have solved the 2d elasticity problem for DBCs and NBCs for a single > material.

[deal.II] Re: Importing .msh file from CUBIT

Himal, It looks like there is a problem opening the file. How did you export the mesh? Are you using read_ucd or read_abaqus? Best, Bruno On Tuesday, December 6, 2022 at 10:36:30 AM UTC-5 075bme0...@pcampus.edu.np wrote: > Warm Greetings!, > I am trying to import .msh file that I have create

Re: [deal.II] Re: Importing .msh file from CUBIT

.inp format file and used read_ucd here. I was actually > trying to import a .msh file. However, the same error occurs for the > read_msh function as well. > > Himal. > > > On Tue, Dec 6, 2022 at 9:41 PM Bruno Turcksin > wrote: > >> Himal, >> >> It lo

[deal.II] Re: Helmholtz solve on GPUs with multicomponent vector

Gourab, Multicomponent vector is currently not supported. Part of the infrastructure is already present as you saw but you cannot use MatrixFree on GPU with multicomponent vectors at the moment. Best, Bruno On Wednesday, December 14, 2022 at 1:19:19 AM UTC-5 pgo...@iisc.ac.in wrote: > I am t

[deal.II] Re: Error in Reading Gmsh File in Dealii

Deepika, Can you try using a linear mesh instead of curvilinear one? Does that work? Best, Bruno On Thursday, December 22, 2022 at 1:47:32 AM UTC-5 deepika...@iitgoa.ac.in wrote: > Hello Everyone, > > I have successfully read a gmsh file for one material in my code. Now I am > trying to read

[deal.II] Re: Regarding deal ii installation with PETSC

Anupama, The error message tells you that PETSc was compiled with MPI support but deal.II was configured without MPI support. Either recompile PETSc without MPI or reconfigure deal.II with MPI support using -DDEAL_II_WITH_MPI=ON Best, Bruno On Tuesday, December 27, 2022 at 3:30:49 AM UTC-5 an

[deal.II] Re: Running non-cuda code when dealii and petsc are compiled with cuda support

Stefano, Currently in deal.II we always assume that the PETSc vector is on the CPU but I don't know what PETSc does internally when it is compiled with GPU support. I don't think that the deal.II developers who use PETSc, use GPU so it's not tested. If you can, compile deal.II with GPU support

[deal.II] Re: Extracting all the nodes on a particular boundary

Wasim, You need to loop over the cells. That's the only way to get the information that you want. Best, Bruno On Thursday, January 19, 2023 at 4:43:34 AM UTC-5 ce21...@smail.iitm.ac.in wrote: > Hello everyone. > I am trying to extract all the nodes present on a boundary with a > particular

Re: [deal.II] Extracting all the nodes on a particular boundary

Nice, I looked for such a function in GridTools but I forgot about DoFTools. Bruno Le jeu. 19 janv. 2023 à 11:44, Wolfgang Bangerth a écrit : > On 1/19/23 02:43, Wasim Niyaz Munshi ce21d400 wrote: > > I am trying to extract all the nodes present on a boundary with a > > particular id. > > Would

[deal.II] Re: Monitor memory usage of the program

Hello, Personally, I prefer to use profilers when monitoring the memory usage. This is because what really matters is how much memory the operating system is allocating for you code. For multiple reasons this can be different than what your code require. Best, Bruno On Tuesday, February 7, 2

[deal.II] Re: facing problem in deal.ii Installation

Hello, It looks like you are trying to install deal.II in /usr/local and you need to be sudo to do it. If you do "ccmake ." in the directory where you configured deal.II you should be able to see what's the path in CMAKE_INSTALL_PREFIX. Make sure that it's what you expect. Best, Bruno On Fri

[deal.II] Re: Creating a plot in dealii

Hello, DataOut is made to add data to a mesh. It cannot be used to do what you want. The easiest way to do what you want is to write the vectors to a file using print() and then do the plot using matlab or python. Best, Bruno On Saturday, February 11, 2023 at 9:27:06 AM UTC-5 ce21...@smail.ii

[deal.II] Re: Spack error with cgal and arbox

Hi, Don't set any variable when using spack. Please post the exact spack command you are using to install deal.II with spack. In general, you want to let spack picks the dependencies for you. I would also advise that you restart from a clean state and remove the configuration files created by

[deal.II] Re: trying to install dealii 9.1.1 using candi.

Hello, It looks like the error is from boost. Which version are you using? What's your compiler? Best, Bruno On Tuesday, March 7, 2023 at 12:58:10 AM UTC-5 me20...@iittp.ac.in wrote: > Dear Sir/Ma'am, > trying to install dealii-9.1.1 version using candi command. i just change > the version i

Re: [deal.II] Re: trying to install dealii 9.1.1 using candi.

Do you need to use Boost in your own project? If you don't, try adding -DDEAL_II_FORCE_BUNDLED_BOOST=ON to your configuration. Best, Bruno Le mer. 8 mars 2023 à 00:02, ME20D503 NEWTON a écrit : > Dear Bruno, > i am using ubuntu 22.04.1 and g++ (Ubuntu 11.3.0-1ubuntu1~22.04) 11.3.0. > > Thank y

[deal.II] Re: Status of hybrid shared/distributed parallelism

Lucas, Why are you interested in hybrid parallelism? Are you hoping to improve the performance of your code or is it simply something you want to try? If the solver is the bottleneck in your code, you should focus on finding a better preconditioner. With that being said, matrix-free methods ten

Re: [deal.II] writing data to a file changes program output although this file is never used

Hello, Usually when I have this kind of bug, there are two possibilities: 1. I am using an un-initialized value 2. I am writing out of bound What I do is using valgrind with my code in Debug mode and without TBB enabled otherwise you get difficult to understand backtrace like here. However,

[deal.II] Re: Linking error when compiling examples: error: cannot find -lhdf5-shared

Amit, Did you load the hdf5 package? You can check using `spack find --loaded`. Best, Bruno On Friday, March 31, 2023 at 4:51:02 AM UTC-4 amitshar...@gmail.com wrote: > Hello everyone, > > I installed dealii using spack manager. I started learning dealii using > examples but I am getting a li

Re: [deal.II] writing data to a file changes program output although this file is never used

Simon, Le ven. 31 mars 2023 à 08:58, Simon Wiesheier a écrit : > > > Given that it works fine in debug mode, there are likely no out-of-bound > violations. I just use std::vector and Vector as data structures > which would be detected in debug mode. > > That's optimistic :) It's not because you

[deal.II] Re: How can I get the jump in my solution if my solution vector is distributed?

Hello, You are using a non-ghosted vector. You need to created a ghosted vector and they assign the non-ghosted vector to the ghosted one. Best, Bruno On Thursday, April 20, 2023 at 7:07:09 AM UTC-4 Abbas wrote: > I am trying to querry for the jump in my solution with > "get_jump_in_functio

[deal.II] Re: How can I visit a face exactly once when running with MPI?

Hello, On Thursday, April 27, 2023 at 7:56:50 AM UTC-4 abbas.b...@gmail.com wrote: This works with a single MPI but not for more. Is it because the cell->face( face)->clear_user_flag() doesn't clear the user flag in the other distributed meshes? That's correct. The flag is only cleared locall

Re: [deal.II] Re: How can I visit a face exactly once when running with MPI?

Right, it won't clear it everywhere. You would need to use MPI to clear the flag on the other processor. If you would do that you would create a huge amount of messages and your code would be extremely slow. Best, Bruno Le jeu. 27 avr. 2023 à 10:09, Abbas Ballout a écrit : > This works!! Thank

[deal.II] Re: Problem compiling deal.II with CUDA

Hello, It looks like the documentation is not up to date. We do not support GPU your GPU, deal.II requires compute capability of at least 6.0. Sorry about that. Best, Bruno On Saturday, April 29, 2023 at 1:46:44 PM UTC-4 sljoh...@gmail.com wrote: > Hey everyone, > > Sorry I have been workin

Re: [deal.II] Re: Problem compiling deal.II with CUDA

No, there is nothing to do. Best, Bruno Le mer. 3 mai 2023 à 14:37, Sean Johnson a écrit : > Thanks! No problem with the documentation not being changed on the compute > capability. With all the changes being made this is really just the > smallest thing. Do I need to mark this solved or anyth

[deal.II] Re: Poor parallel efficiency of FEValues::reinit() with Workstream

Jonas, What you are seeing is expected. There is an overhead associated with using multiple threads. Since the overhead is constant it looks like it disappears when you do more work. You probably want to give more work do to your threads to amortize the constant cost. Best, Bruno On Thursday

[deal.II] Re: CUDA facilities

Giuseppe, I am not sure if google sent my previous message, so I will send it again. I apologize if you get it twice. We are currently replacing the CUDA code with Kokkos. These will allow us to support Nvidia, AMD, and Intel GPUs. It will also make it easier to add new functionalities. Unfor

Re: [deal.II] Re: CUDA facilities

Kokkos keywords’, but I cannot figure out from a rapid glance if it > can be already executed with a version of deal.II installed with > ‘DEAL_II_WITH_KOKKOS’ flag on or one has to wait the new release. > > > > Best > > > > Giuseppe > > > > *Da: *Bruno

[deal.II] Re: evaluation FE function outside domain

Math, The basis functions are not defined outside of an element, so you cannot extrapolate them outside of the mesh. What are you trying to achieve? Best, Bruno On Monday, June 5, 2023 at 4:59:28 PM UTC-4 Mathieu wrote: > Hello everyone, > > I was wondering if deal.II has built-in functions

Re: [deal.II] Re: evaluation FE function outside domain

I do not know if that's reasonable. > > Best, > > Math > > Bruno Turcksin schrieb am Di., 6. Juni 2023, > 21:55: > >> Math, >> >> The basis functions are not defined outside of an element, so you cannot >> extrapolate them outside of the mesh. What

[deal.II] Re: dealii v 9.2.0 installation difficulty (also posted on spack group)

Andrew, What's the error message that you get when using `spack install dealii@9.2.0 %gcc@12.2.0` ? spack writes the exact error in an output file. I don't think changing the C++ standard is the right fix, instead I would change the version of boost. Best, Bruno On Friday, June 23, 2023 at 1

Re: [deal.II] Re: dealii v 9.2.0 installation difficulty (also posted on spack group)

The error is from here In file included from ./boost/math/tools/mp.hpp:341, from ./boost/math/policies/policy.hpp:11, from ./boost/math/special_functions/math_fwd.hpp:31, from ./boost/math/special_functions/airy.hpp:11, from ./boo

Re: [deal.II] Re: dealii v 9.2.0 installation difficulty (also posted on spack group)

It looks like the error you get when using cxx=11 is from ginkgo so if you don't need ginkgo, you could try something like spack install dea...@9.2.0 ~ginkgo ^bo...@1.76.0 cxxstd=11 %g...@12.2.0 Bruno Le ven. 23 juin 2023 à 14:47, Andrew Li a écrit : > Unfortunately dealii@9.2.0 seems to be in

[deal.II] Re: Searching for libs in BLAS_LIBRARY_DIRS=''

Phil, did you forget to change MKL=OFF to MKL=ON, in candi.cfg (https://github.com/dealii/candi/blob/master/candi.cfg#L68)? This should change the name from blas to whatever is the name for mkl. Best, Bruno On Thursday, July 6, 2017 at 9:28:22 AM UTC-4, Phil H wrote: > > Hello I'm using candi

Re: [deal.II] Re: Searching for libs in BLAS_LIBRARY_DIRS=''

gt; No, I changed the MKL to ON, which is why this is strange. > > Cheers, > > Phil > > > On Thursday, July 6, 2017 at 10:14:53 AM UTC-4, Bruno Turcksin wrote: >> >> Phil, >> >> did you forget to change MKL=OFF to MKL=ON, in candi.cfg >> (https://git

[deal.II] Re: CMake extra libraries?

Victor, On Thursday, July 6, 2017 at 3:22:04 PM UTC-4, Victor Eijkhout wrote: > > add the flag that you need in SET(Trilinos_CXX_COMPILER_FLAGS) > > build@build-BLDCHROOT:12.10.1> find . -name \*.cmake -exec grep > COMPILER_FLAGS {} \; -print >

[deal.II] Re: Example for three coupled equations in dealII?

Hey, On Friday, July 7, 2017 at 5:03:05 AM UTC-4, Maxi Miller wrote: > > I have the following three PDEs (reduced from a more complex version): > > >

Re: [deal.II] Re: CMake extra libraries?

This is very strange. Do you set the AVX flags yourself? @Martin have you ever seen something like that? Bruno 2017-07-07 3:50 GMT-04:00 Victor Eijkhout : > Hm. That second compile line should be: > > > > > > > > > > [ 63%] Building CXX object > source/numerics/CMakeFiles/obj_numerics_debug.dir/

Re: [deal.II] Re: CMake extra libraries?

2017-07-07 15:42 GMT-04:00 Victor Eijkhout : > It's weird. I just noticed taht I didn't have any -O2 flags in trilinos, so > I added. And now the deal install has compile lines with O2 and the AVX > flags, but also a bunch of compile lines with O0. And only sometimes those > blow up: Can you config

[deal.II] Re: Combination of Trilinos and deal.II fails at compilation

Hi, do you have the same problem with openmpi or mpich? If I remember correctly, we had problems in the past with Intel MPI Best, Bruno On Saturday, July 8, 2017 at 1:55:46 PM UTC-4, Maxi Miller wrote: > > When compiling deal.II without Trilinos, the compilation works. When > including Trili

Re: [deal.II] Re: Combination of Trilinos and deal.II fails at compilation

2017-07-08 16:40 GMT-04:00 'Maxi Miller' via deal.II User Group : > Should I test that (i.e. recompile all libraries using OpenMPI), or are > there other possible fixes? Yeah, I would try using OpenMPI. Honestly, I have never seen this error before, so I don't know what else to do. Best, Bruno -

Re: [deal.II] Re: CMake extra libraries?

Victor, 2017-07-10 10:28 GMT-04:00 Victor Eijkhout : > Et cetera. We have boost 1.64. Any idea? Make sure that deal.II is compiled with only one version of boost (the same than the one you used for Trilinos if you compiled Trilinos with boost). It already happened to me that I had two versions of

Re: [deal.II] Re: CMake extra libraries?

Victor, I have never seen anything like that before :-/ The only thing I found is this (very) old thread without a reply https://cmake.org/pipermail/cmake/2006-April/008893.html Someone got a similar problem than yours when updating cmake. So maybe try another version of cmake. Best, Bruno --

[deal.II] Re: Modification of right hand side of example 52 for multiple equations leads to linker error

Maxi, like Jean-Paul already said we don't instantiate the function for Table so you need to add #include This way you can use any data type you want. Best Bruno On Wednesday, July 19, 2017 at 10:43:23 AM UTC-4, Maxi Miller wrote: > > > I tried to modify the right side of example 52 from one

Re: [deal.II] Re: Modification of right hand side of example 52 for multiple equations leads to linker error

2017-07-19 14:40 GMT-04:00 'Maxi Miller' via deal.II User Group : > But is that even a good idea to put the three solution-vectors into one > Table, and return them, or should I rather put all three vectors in one > large vector, and then slice the respective results out of the large vector? Yeah u

[deal.II] Re: Defining "complex" geometry for hand drawn mesh

Philip, On Monday, July 31, 2017 at 1:54:26 PM UTC-4, phillip mobley wrote: > > 1) Based on this post, do you think that I should attempt to use dealli > internal mesh generation or would it be better to utilize gmsh source code > and intergrate it into my code? (As a side note, I am not sure ye

[deal.II] Re: How to clear data of triangulation completely

Felix, it should work. I have done something similar in the past. Can you send a small example that shows the problem. Best, Bruno On Tuesday, August 1, 2017 at 1:26:00 PM UTC-4, Felix Lorenz wrote: > > Hello everybody, > I have the following problem: > > Class C{ > ... > > Triangulation <

[deal.II] Re: about deal.II with CUDA C programming acceleration

Chih-Che, The CUDA support in deal.II is very new. It is only in the development version of deal (https://github.com/dealii/dealii). You can see our current development plan here https://github.com/dealii/dealii/projects/2 Right now, we have support for vector and partial support for matrix-fre

[deal.II] Re: errors when installing dealii-8.5.1 after successful cmake configuation

Hello, On Wednesday, August 16, 2017 at 7:45:25 AM UTC-4, 曾元圆 wrote: > ### > -- Configuring done > -- Generating done > -- Build files have been written to: /mnt/hgfs/My code/dealii/build > On linux it is a bad idea to have a blank space in the name of a directory. Instead of /mnt/hgfs/My code/

<    1   2   3   4   5   6   7   8   >