Re: [deal.II] Extract some blocks from a sparse block matrix

2017-02-27 Thread lianzhang030592
Dear Timo Heister, Thank you very much. I have used the BlockMatrixArray object in my code and it works well. However, it seems that the BlockMatrixArray object is a kind of pointer so that we have to use the iteration method to solve the linear system, of course, we could develop precondition

Re: [deal.II] InternalError in VectorTools::compute_nonzero_normal_flux_constraints

2017-02-27 Thread Wolfgang Bangerth
Lars, I don't think anything particularly complicated is necessary -- you may just have to copy the support points from the base elements (or generalized points if the regular ones are not available) to the FESystem's ones. The interleaving will work just like for all of the oth

Re: [deal.II] write_vtk fail when using hp functionality

2017-02-27 Thread Wolfgang Bangerth
On 02/27/2017 12:24 PM, Weixiong Zheng wrote: Dear all, Might be a trivial questions to you guys. I get a scalar problem using hp class. When I tried to output the results using | std::vectorsolution_names; | solution_names.push_back ("phi"); std::vector data_component_interpretati

[deal.II] write_vtk fail when using hp functionality

2017-02-27 Thread Weixiong Zheng
Dear all, Might be a trivial questions to you guys. I get a scalar problem using hp class. When I tried to output the results using std::vector solution_names; solution_names.push_back ("phi"); std::vector data_component_interpretation(1, DataComponentInterpretation::component_i