Re: [deal.II] Re: ifem

Dear Jean-Paul, Thank you for your reply. I tried with the second alternative you give me. According to the example, contract<0,2> (t1, t2) , I modify the code as follows: PeFT = contract<2,2> (Pe[qs], F[qs]); I interprete that PeFT is the result, Pe[qs] would be the tensor t1, and

[deal.II] Re: ifem

Dear Joaquin, The warning appears because you are calling a function that is marked as deprecated, so although it remains in the library for now it will be removed in the near future. There are two a

Re: [deal.II] Re: Tensor output via ConditionalOStream pcout

On 12/14/2016 03:42 PM, seyedali88 via deal.II User Group wrote: Thanks for the hint, but I am still confused. Did no one ever try to output a sparse matrix in parallel? Isn't there a single example of how to output the system_matrix? How should someone debug or check the code properly... There

[deal.II] ifem

Dear all, I am trying to run the IFEM code on FSI problems, developed by Luca Heltai. I am using deal.II 8.4.1. I edited some parts of the code due to warnings because it is considered deprecated for the version I am using. I stopped when warning appears for line 1038: 1037 if((!p

[deal.II] Re: Tensor output via ConditionalOStream pcout

Thanks for the hint, but I am still confused. Did no one ever try to output a sparse matrix in parallel? Isn't there a single example of how to output the system_matrix? How should someone debug or check the code properly... BR, Seyed Ali Mohseni -- The deal.II project is located at http://www

[deal.II] Re: Best application to plot

Hi, Don't know about Visit, but in Paraview you can also plot-over-line and then export the data and re-plot it using Gnuplot with right styles, lines, fonts, etc. But otherwise, I never cam across any instance where Vist and/or Paraview would not be enough. Cheers, Denis. On Wednesday, Dece

Re: [deal.II] Best application to plot

On 12/14/2016 12:42 PM, Sudarshan Kumar wrote: Apart from VisIt and Paraview, could any one tell, which is the best to plot .vtu files of 2D data. Why would you not want to use these two? I sometimes use gnuplot to get a quick overview, but visit/paraview are so much better for real pictur

[deal.II] Best application to plot

Apart from VisIt and Paraview, could any one tell, which is the best to plot .vtu files of 2D data. Thanks a lot, Best -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this me

Re: [deal.II] build_patches(degree) for higher order dg elements in parallel

Wolfgang, I've attached a small test case that simply initializes a parallel vector and then attempts to write it to file via DataOutFaces. The finite element used only lives on the faces of the cells so DataOutFaces in appropriate. This test case provided me with some insight into my problem.

Re: [deal.II] Re: BlockSparseMatrix and BlockVector with a nested hierarchy of blocks

On 12/14/2016 03:01 AM, Jean-Paul Pelteret wrote: I think that the BlockMatrixArray class might offer the functionality that you're looking for? Yes, that's exactly the way to go: chop the matrix into its smallest compo

Re: [deal.II] using FEValues extractor for a

Anup, I am solving for a 3 dimensional displacement vector (*dim+1* dimensional) in a *dim=2* dimensional space (for a planer problem in elasticity). Hence I consider an FEvaluesExtractor: const FEValuesExtractors::Vector u_fe; and initialize it as u_fe(0) Now I need to compute the gradie

[deal.II] Re: Tensor output via ConditionalOStream pcout

Dear S.A. Mohseni, You need to look carefully at the documentation for the specific SparseMatrix class. Its unclear to me whether LA::MPI::SparseMatrix represents a TrilinosWrappers::SparseMatrix

Re: [deal.II] distributed body force

On 12/13/2016 03:28 PM, benhour.amiria...@gmail.com wrote: It would be very kind of you if you let me know how I can imply distributed load on the inner radius of a ring instead of using nodal load. Benhour -- as I and others have pointed out on previous occasions, this question is impossible

Re: [deal.II] Alternative for ConstraintMatrix::condense in parallel codes?

On 12/13/2016 02:54 PM, Hamed Babaei wrote: The problem is that*constraints.condense* doesn't work with *TrilinosWrappers::SparseMatrix *so the question is when we can not use *constraints.distribute_local_to_global *how we can apply constraints in the assembly for parallel codes. In the https:/

[deal.II] Re: Tensor output via ConditionalOStream pcout

Hi again, It's not quite clear to me yet. Is it possible to pcout a MPI SparseMatrix or Vector? If I do e.g. pcout << system_matrix << std::endl; Then I receive the address as output since the system matrix is initialized like this: LA::MPI::SparseMatrix system_matrix; How is it possible to

[deal.II] Re: Literature sources of Krylov subspace solvers

Hello Bastian, thank you for your advice. Like stated in an answer before, I was already recommended to this book. But it doesn't hurt to ask, so I did it. Kind regards, Ingo Am Mittwoch, 14. Dezember 2016 14:40:33 UTC+1 schrieb Bastian: > > Hi Ingo, > > this one is in German: > > http://link.

[deal.II] Re: Literature sources of Krylov subspace solvers

Hello Bruno, thank you for your fast answer. I already have been recommended to a book (A. Meister, Numerik linearer Gleichungssysteme) by a colleague, but I just wanted to be sure - if you were wondering why I posed the question in the first place. Yours sincerely, Ingo Am Mittwoch, 14. Dez

Re: [deal.II] Duplication access to the cell?

Thank you very much!!! The problem is solved, due to your help!! Kyusik. 2016년 12월 14일 수요일 오전 1시 55분 53초 UTC+9, Wolfgang Bangerth 님의 말: > > Kyusik > Material ids are only stored as unsigned char objects, which can only > represent values up to 255. You may want to use the user_index instead. >

[deal.II] Re: Literature sources of Krylov subspace solvers

Hi Ingo, this one is in German: http://link.springer.com/book/10.1007%2F978-3-658-07200-1 Best, Bastian -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you ar

[deal.II] Re: Literature sources of Krylov subspace solvers

Ingo, Krylov solvers have been around for a while and the implementation in deal.II is pretty standard. Here is a reference that a few people use https://web.stanford.edu/class/cme324/saad.pdf but you can just use google//google scholar to get more references. Best, Bruno On Wednesday, Decem

[deal.II] Literature sources of Krylov subspace solvers

Hello, I'm going to write my master thesis about efficient preconditioning and solving of a saddle point-like system arising in thermoelastic modelling of drilling, implemented in deal.ii. Therefore I need references to shortly review the appropriate Krylov subspace methods (BiCGStab, GMRES and

[deal.II] Re: BlockSparseMatrix and BlockVector with a nested hierarchy of blocks

Dear Oded, I think that the BlockMatrixArray class might offer the functionality that you're looking for? If so then I'm not quite sure what the equivalent class for the BlockVector would be. If there is none, then you h