Re: [deal.II] Can one FE be used over discontiguous domain

2017-01-05 Thread Weixiong Zheng
All right, I will try and come back if something's wrong. Thanks, sir. 在 2017年1月5日星期四 UTC-6下午9:08:56,Wolfgang Bangerth写道: > > On 01/05/2017 08:07 PM, Weixiong Zheng wrote: > > > > I am coding up a code such that DG is used in one subdomain, while CG is > used > > in all other subdomains. The i

Re: [deal.II] Can one FE be used over discontiguous domain

2017-01-05 Thread Wolfgang Bangerth
On 01/05/2017 08:07 PM, Weixiong Zheng wrote: I am coding up a code such that DG is used in one subdomain, while CG is used in all other subdomains. The implementation is with hp::FECollection class. As the CG-subdomains are discontiguous, i.e. they are spatially disconnected, can I still use on

[deal.II] Can one FE be used over discontiguous domain

2017-01-05 Thread Weixiong Zheng
Dear all, I am coding up a code such that DG is used in one subdomain, while CG is used in all other subdomains. The implementation is with hp::FECollection class. As the CG-subdomains are discontiguous, i.e. they are spatially disconnected, can I still use only one FE_Q for all disconnected C

[deal.II] Re: Periodic Boundary Condition

2017-01-05 Thread Hamed Babaei
Dear Daniel, You can obtain local DoF values from a global FE vector using > FEValuesBase::get_function_values [1], create the right-hand side locally > as you wish and use ConstraintMatrix::distribute_local_to_global to create > the global right-hand side. This is outlined in step-21 for exam

Re: [deal.II] Boundary ID in a parallel distributed triangulation & Neumann boundary conditions

2017-01-05 Thread Timo Heister
> The following line gives me a headache: > > std::cout << cell->face(face_number)->boundary_id() << std::endl; > > ... since I cannot see anything. It outputs strange symbols such as � or > some 0001 cube. just do std::cout << (int)cell->face(face_number)->boundary_id() << std::endl; -- Timo H

[deal.II] Boundary ID in a parallel distributed triangulation & Neumann boundary conditions

2017-01-05 Thread Seyed Ali Mohseni
Hi, Sorry to disturb again, I have a problem. :) The following line gives me a headache: std::cout << cell->face(face_number)->boundary_id() << std::endl; ... since I cannot see anything. It outputs strange symbols such as � or some 0001 cube. My aim is to somehow write my RHS for distribu

[deal.II] Re: Periodic Boundary Condition

2017-01-05 Thread Daniel Arndt
Hamed, [...] > However, to be able to use ConstraintMatrix::distribute_local_to_global, I > can assemble system_matrix in local level instead of global level. Having > that said, I can not do so for the system_rhs because solution and > old_solution vectors are global vectors. > You can obtai