Re: [deal.II] FullMatrix to SparseMatrix

2020-10-08 Thread Nikki Holtzer
When I do the above, with the modifications you have provided, I still obtain the following: error: no match for ‘operator+’ (operand types are ‘const dealii::SparseMatrix’ and ‘dealii::SparseMatrix’) or error: no matching function for call to ‘dealii::SparseMatrix::add(double,

Re: [deal.II] Making a pair of extracted data and cast into a set

2020-10-08 Thread Daniel Arndt
Behrooz, Can you elaborate some more on what you are trying to achieve? The (pseudo-)code you posted looks OK apart from vertices being declared in the innermost loop. You probably want it to be outside the loop over the vertices of a face. The API for Point can be found at

Re: [deal.II] FullMatrix to SparseMatrix

2020-10-08 Thread Nikki Holtzer
Hello everyone, I was just curious if anyone had any suggestions for the above 'const' qualifier problem I have stated above? Thanks! Nikki On Thursday, October 8, 2020 at 9:16:24 AM UTC-4 Nikki Holtzer wrote: > When I do the above, with the modifications you have provided, I still > obtain

Re: [deal.II] FullMatrix to SparseMatrix

2020-10-08 Thread Wolfgang Bangerth
On 10/8/20 7:16 AM, Nikki Holtzer wrote: error: no match for ‘operator+’ (operand types are ‘const dealii::SparseMatrix’ and ‘dealii::SparseMatrix’) SparseMatrix has no operator+, so this will not work. error: no matching function for call to ‘dealii::SparseMatrix::add(double,

Re: [deal.II] Making a pair of extracted data and cast into a set

2020-10-08 Thread Behrooz Karami
Dear Daniel, Thanks very much for your reply. Actually what I want to do is very simple. Consider that the above code gives me 3 nodal indices as: 1, 4 and 5. First, I need to be able to cast them into a set or vector. We can call it set1. Then assume that I have another set, e.g. set2, which

[deal.II] "locked" vtu outputs

2020-10-08 Thread David Montiel Taboada
Hello everyone, We are using PRISMS-PF, a phase field framework finite element built upon the deal.II library (with the MPI and p4est) dependencies. The code uses the deal.II DataOut class (see code snippet at the end of the email) to output simulation results into vtu files at certain time

Re: [deal.II] Use of QGauss object in tutorials

2020-10-08 Thread Xuefeng Li
On Wed, Oct 7, 2020 at 7:19 PM Wolfgang Bangerth wrote: > > Xuefeng, > > > Still on the course of learning dealii by reading through the tutorials. > Came > > up with the following observation. > > > > 1. In Step-4, a single QGauss object was created in assemble_system(), > and > > was used

[deal.II] Making a pair of extracted data and cast into a set

2020-10-08 Thread Behrooz Karami
Hi everyone, I am trying to cast extracted vertex_indices into a set (or vector). Before that it is needed to make a pair of those indices (just doubling) and to cast the pairs into a set for further manipulations. I get the vertex information mainly through following lines of code: for