[deal.II] Experimental data at output boundary

2018-03-06 Thread 'Joaquin' via deal.II User Group
Dear all, I am trying to insert experimental data from a file at the output boundary of a tube. Is it possible to do that in deal.II? I attach the file of experimental data, pressure vs time. Thanks for any help. Joaquín -- The deal.II project is located at http://www.dealii.org/ For mailing

Re: [deal.II] Are shell elements available in deal.II?

2018-03-06 Thread Wolfgang Bangerth
I am excited to learn that there are arguably better methods (mixed methods) to solve shell problems. I am not familiar with those; could you please kindly provide some pointers for me to start (anything, either example, code, or paper)? That all depends on what equation exactly you're lookin

Re: [deal.II] Are shell elements available in deal.II?

2018-03-06 Thread Yuxiang Wang
Thank you Professor Bangerth for your response! Yes that's exactly what I am interested in. I am excited to learn that there are arguably better methods (mixed methods) to solve shell problems. I am not familiar with those; could you please kindly provide some pointers for me to start (anything

Re: [deal.II] Re: Iterating over all the entries in a PETScWrapper::MPI::SparseMatrix in parallel

2018-03-06 Thread Feimi Yu
Hi Wolfgang, This time I used VectorOperation::insert and it didn't happen the memory error that I posted before. However, if I put the compress function after the loop, it only sets one entry on each rank then throws the exception "Object is in wrong state, not for unassembled matrix." Putti

Re: [deal.II] some problems concerning step-14

2018-03-06 Thread 曾元圆
Thank you so much Bangerth! Now I understand why we need to rewrite the error formula on a cell as residual times dual weight. But I'm still a little confused with the reason why we must introduce z_h. Just as you mentioned, if we introduce z_h, then z-z_h is a quantity that is only large wher

Re: [deal.II] Re: step-22 partial boundary conditions

2018-03-06 Thread Wolfgang Bangerth
Jane, I believe the formula is correct. the cubic term comes from p=z^3 being the pressure manufactured solution. so in (pI-2e) you get a z^3 term and indeed a linear term in the 2e portion. I see. I missed that there is of course also a pressure. This is how I have come to conclude that i

Re: [deal.II] Re: Iterating over all the entries in a PETScWrapper::MPI::SparseMatrix in parallel

2018-03-06 Thread Wolfgang Bangerth
On 03/05/2018 02:54 PM, Feimi Yu wrote: I changed my strategy to use set(r, c, v) function to set the values so that I can use the const iterators. also called compress after every add: for(autor =Abs_A_matrix->block(0, 0).local_range().first; r block(0, 0).local_range().second; ++r) { for(aut

Re: [deal.II] Are shell elements available in deal.II?

2018-03-06 Thread Wolfgang Bangerth
Shawn, I tried to search but did not find an implementation of shell elements in deal.II. Since this is a commonly used element, I'd like to make sure that it's not me missing it. Could you please help confirm? That depends on how exactly you define "shell elements". But if you mean "elemen

[deal.II] Re: Cell->measure returning Nan

2018-03-06 Thread Moraad Biagooi
I had a similar issue (returning NaN) when dealing with a 3D HyperBall mesh. I realized the problem was that the face I was trying to measure was not planar enough. There's a tolerance condition that has to be fulfilled at this line of the measure function. https://www.dealii.org/8.4.0/doxygen/

[deal.II] Re: Cell->measure returning Nan

2018-03-06 Thread Moraad Biagooi
I had a similar issue when dealing with a 3D HyperBall mesh. I realized the problem was that the face I was trying to measure was not planar enough. There's a tolerance condition that has to be fulfilled at this line of the measure function. https://www.dealii.org/8.4.0/doxygen/deal.II/tria__acc