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

2016-12-13 Thread Wolfgang Bangerth
On 12/13/2016 10:13 AM, Justin Kauffman wrote: Upon further investigation you are correct. If I go beyond 4 processors I get the same error even with the default zero subdivisions. Justin -- I suspect that the class simply does not support parallel triangulations. Do you have any suspicion th

Re: [deal.II] Unsuccessful integration of p4est with dealii

2016-12-13 Thread Kartik Jujare
Thank you for your help. On Tuesday, December 13, 2016 at 8:39:40 AM UTC+1, Uwe Köcher wrote: > > Dear Kartik, > > well, candi has users around the world - if it reports a successful > installation, then everything should be fine. > But anyhow, you can go to each installation folder, e.g. deal.II

[deal.II] using FEValues extractor for a

2016-12-13 Thread Anup Basak
Hello all, 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 g

[deal.II] distributed body force

2016-12-13 Thread benhour . amirian66
Dear All, 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. Thanks, Benhour -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/foru

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

2016-12-13 Thread Hamed Babaei
Hi all, I have used* constraints.condense* function in a serial code to apply periodic constraints on system_matrix and system_rhs after their assembly in which I do not use *constraints.distribute_local_to_global* for cell_matrix and cell_rhs but I first assemble mass_matrix , laplace_matrix,

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

2016-12-13 Thread Justin Kauffman
Wolfgang, Upon further investigation you are correct. If I go beyond 4 processors I get the same error even with the default zero subdivisions. Thank you for your response. - Justin On Tuesday, December 13, 2016 at 12:00:18 PM UTC-5, Wolfgang Bangerth wrote: > > Justin > Is it really working

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

2016-12-13 Thread Wolfgang Bangerth
JustinIs it really working with zero subdivisions? Step 40 only uses the DataOut class, not the DataOutFaces class, so that may not be the best comparison. I don't recall whether the latter class works in parallel at all, regardless of subdivision. BestWolfgang Sent from Samsung tablet.

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

2016-12-13 Thread Wolfgang Bangerth
KyusikMaterial 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.  BestWolfgang Sent from Samsung tablet. Original message From: hanks0...@gmail.com Date: 12/13/2016 05:20 (GMT-07:00) T

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

2016-12-13 Thread Justin Kauffman
Hi all: I am writing an HDG code and trying to parallelize it in a fashion similar to step-40. I am having a problem with data_out.build_patches(degree). In step-40 the argument is the default value of zero, but I eventually want to run this code for higher order dg elements. I believe that the

[deal.II] Duplication access to the cell?

2016-12-13 Thread hanks0227
Dear All, I'm trying to calculate something differently depending on material id. So, I have to set material id before next iteration.(I'm going to use picard iteration) The way I have done is 1.the middle point of mesh would be selected. (it is indicated as max_cell in my code) 2.the neighb