[deal.II] Re: Neumann boundary condition for the curved domain

Dear Wolfgang, It should not be zero. In addition, the local quadrature point is for the cell. Thanks, Benhour On Thursday, December 8, 2016 at 8:16:03 PM UTC-6, benhour@gmail.com wrote: > > Dear All, > I still have problem with implying boundary condition on the curved > domain. consider

Re: [deal.II] Doubt with set_boundary_id()

Hello Professor, The problem was because of the mesh. I regenerated the mesh and this time there is no such issue that arises. I think the options used to export and make the mesh earlier were wrong. Thanks for your help again. On Thursday, December 8, 2016 at 11:26:49 PM UTC-5, RAJAT ARORA

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

If someone is interested, I also wrote some code for rank 4 tensors since the pcout version only outputs everything to a single line: template void SolidMechanics::print_tensor (const SymmetricTensor<4, dim>& tensor, char* name) { pcout << std::endl << name << std::endl; for (unsigned int i =

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

Interesting. The ConitionalOstream pcout can directly output such things. I thought it is just created for replacement of cout to reduce parallelized output. Very nice feature. Thank you, I didn't know. BR, Seyed Ali Mohseni -- The deal.II project is located at http://www.dealii.org/ For

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

Interesting. The ConitionalOstream pcout can directly output such things. I thought it is just created for replacement of cout. Very nice feature. Thank you, I didn't know. BR, Seyed Ali Mohseni -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see

[deal.II] Re: Periodic boundary conditions seem to be not applied

Dear Daniel, > I would suggest you to look at the constraints in the ConstraintMatrix > using ConstraintMatrix::print[1] for a mesh consisting of has few cells as > possible. > Using DoFTools::map_dofs_to_support_points[2] you can then verify that > these are correct. > Issuing "

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

Dear S. A. Mohseni, I found these definitions within symmetric_tensor.h. How can I access them? > They seem to allow output for my case. > [...] > You would just use the << operator, i.e. pcout << tensor Best, Daniel -- The deal.II project is located at http://www.dealii.org/ For mailing

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

I found these definitions within symmetric_tensor.h. How can I access them? They seem to allow output for my case. /** * Output operator for symmetric tensors of rank 2. Print the elements * consecutively, with a space in between, two spaces between rank 1 * subtensors, three between rank 2

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

Yes, well TrilinosWrappers::SparseMatrix only has a print() function. I took from what you wrote about that you were using dealii:: SparseMatrix. As could be expected, you

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

I just wrote a code for rank 2 tensors, maybe it's useful for someone: template void SolidMechanics::print_tensor (const SymmetricTensor<2, dim>& tensor, char* name) { pcout << std::endl << name << std::endl; for (unsigned int i = 0; i < dim; ++i) { for (unsigned int j = 0; j < dim; ++j) pcout

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

It says the following after I run your suggestion: class "dealii::TrilinosWrappers::SparseMatrix" has no member "print_formatted" system_matrix.print_formatted(pcout.get_stream(),2,true,2,"0",1); Am Freitag, 9. Dezember 2016 16:46:26 UTC+1 schrieb mohseng...@gmail.com: > > Hi, > > I am

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

Dear S. A. Mohseni, >From the class documentation , I figure that you should be able reach the equivalent output with (something like) the following steps: if (pcout.is_active()) { system_rhs.print(pcout.get_stream());

[deal.II] Tensor output via ConditionalOStream pcout

Hi, I am wondering, if there is a possibility to output matrices and vectors for parallel code. For instance I already read in previous topics here that you can use: system_matrix.print_formatted(cout,2,true,2,"0",1) or system_rhs.print(cout). Unfortunately, it is not possible to use pcout