Re: [deal.II] Time dependent Step 46, implementing non-homogeneous kinetic coupling

2024-12-18 Thread Wolfgang Bangerth
On 12/18/24 16:00, Sclah wrote: for (unsigned int i = 0; i < local_face_dof_indices_fluid.size(); ++i) { if (fe_fluid->face_system_to_component_index(i).first < dim) // select just fluid velocity (the solution is stored as (fluid_velocity + fluid_pressure + solid_displacement) with (dim + 1 +

Re: [deal.II] Time dependent Step 46, implementing non-homogeneous kinetic coupling

2024-12-18 Thread Sclah
I need them just on the faces that separate the materials (fluid and solid). The piece of code is the following, on an interface cell face I have: cell->face(face_number)->get_dof_indices(local_face_dof_indices_fluid, fluid_fe_index); // select the fluid dofs sets from hp-dofs cell->face(face_

Re: [deal.II] Time dependent Step 46, implementing non-homogeneous kinetic coupling

2024-12-18 Thread Sclah
I need them just on the faces that separate the materials (fluid and solid), the get_dof_indices. The piece of code is the following, on an interface cell I have: cell->face(face_number)->get_dof_indices(local_face_dof_indices_fluid, fluid_fe_index); // select the fluid dofs sets from hp-dofs c

Re: [deal.II] Handling NaN Values in TrilinosWrappers::SparseMatrix with Amesos_SuperluDist

2024-12-18 Thread Wolfgang Bangerth
Simon: Unfortunately, the program results in a segmentation fault. Specifically, no exception (TrilinosWrappers::SolverDirect::ExcTrilinosError) is thrown when the matrix containing NaN values is passed to the solver's solve function. If you run the program in a debugger, can you find out w

[deal.II] Handling NaN Values in TrilinosWrappers::SparseMatrix with Amesos_SuperluDist

2024-12-18 Thread Simon
Dear all, In the attached example, I deliberately assign NaN values to a TrilinosWrappers::SparseMatrix to test how the Amesos_SuperluDist direct solver handles such inputs. The code is run using deal.II version 9.4.0. Unfortunately, the program results in a segmentation fault. Specifically, n

[deal.II] Support for Adaptive Mesh Refinement on Simplex Meshes (triangles) in deal.II

2024-12-18 Thread YI-YUNG YANG
Hi everyone, I was wondering if deal.II currently supports adaptive mesh refinement for simplex meshes, specifically for triangular elements? If so, are there any examples or documentation that demonstrate its implementation? Thanks for your help! Best, Yi-Yung Yang -- The deal.II projec