[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 (auto

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 https://www.dealii.org

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 is

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

2020-10-09 Thread Behrooz Karami
By the way, to be more specific, a simple command like the following : std::pair , Point> mypair (vertices[v], vertices[v]); (in conjunction with vertices[v] = cell->face(f)->vertex(v);) makes a pair of coordinates of vertices. However, when I modify it towards vertex_index(v), it doesn't work: v

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

2020-10-09 Thread Jean-Paul Pelteret
Hi Behrooz, I must admit that I’ve not been able to follow why you want to store Points when you’re interested in global vertex indices, but I’ll try to give a little help where I can. In case its of any help to you, I believe that you should be able to get the coordinate of the vertex with a k

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

2020-10-13 Thread Behrooz Karami
Dear Jean-Paul, Thanks very much for your explanations. I managed to achieve my objective. You were right, one of the crucial points I had to take into account was the compatibility of data structures. Defining a new operator was another to the point advice, as some of the errors were laid on t