Re: [deal.II] [Assigning boundary value to node]

2023-10-31 Thread Wolfgang Bangerth
Veerendra: I'm currently working on solving an elasticity problem based on the step-8 tutorial. In my simulation, I've successfully prescribed boundary values for an edge of the triangulation using the following code:  VectorTools::interpolate_boundary_values( mapping, dof_handler,      

[deal.II] Nan for stress values

2023-10-31 Thread Raghunandan Pratoori
Hello, I am running a code to study phase transformations in metals. The code runs without any problem in release mode, but gives nan when outputting stress values using the same code snippet as in Step-18. When I run the same code in debug mode, it gives the following error: An error occurred

Re: [deal.II] Using KINSOL with distributed vectors

2023-10-31 Thread Norihiro Watanabe
Hi Yu, The code I wrote are something like this: TrilinosWrappers::MPI::Vector vector_without_ghost; vector_without_ghost.reinit(jacobian.locally_owned_range_indices(), jacobian.get_mpi_communicator()); vector_without_ghost = vector_with_ghost; In the above, I used a jacobian matrix object to g

Re: [deal.II] "is protected within this context" error for BlockSparseMatrix iterators

2023-10-31 Thread Wolfgang Bangerth
On 10/31/23 13:34, Simon Wiesheier wrote: I have to stick to deal.ii 9.4.0 for a while which I installed via spack. Is there anything I can add to my user code or do I have to recompile the library with your changes? If you are able to edit the header files of your installation, you can ju

Re: [deal.II] "is protected within this context" error for BlockSparseMatrix iterators

2023-10-31 Thread Simon Wiesheier
Thank you for the fix! I have to stick to deal.ii 9.4.0 for a while which I installed via spack. Is there anything I can add to my user code or do I have to recompile the library with your changes? Best Simon Wolfgang Bangerth schrieb am Di., 31. Okt. 2023, 20:17: > > > dealii::BlockSparseMatr

Re: [deal.II] "is protected within this context" error for BlockSparseMatrix iterators

2023-10-31 Thread Wolfgang Bangerth
dealii::BlockSparseMatrix systemMatrix; auto systemMatrixIteratorEnd = systemMatrix.end(0); auto systemMatrixIteratorBegin = systemMatrix.begin(0); if (systemMatrixIteratorBegin != systemMatrixIteratorEnd) std::cout<<"Example"< Using deal.ii 9.4.0, this results in the compilation error inclu

[deal.II] "is protected within this context" error for BlockSparseMatrix iterators

2023-10-31 Thread Simon
Dear all, dealii::BlockSparseMatrix systemMatrix; auto systemMatrixIteratorEnd = systemMatrix.end(0); auto systemMatrixIteratorBegin = systemMatrix.begin(0); if (systemMatrixIteratorBegin != systemMatrixIteratorEnd) std::cout<< "Example"<::operator==(const dealii::BlockMatrixIterators::Accessor&)