Dear All. Does it possible to make a funtion as
void ConstraintMatrix::distribute_local_to_global ( const FullMatrix< double > & local_matrix, const std::vector< unsigned int > & row_indices , const std::vector< unsigned int > & col_indices , MatrixType & global_matrix ) const [inline] The third argument should be added to what it was before, which point out to be column indices. I am doing with DG (based on tutorial 11) + periodic boundary via Constrant Matrix as given in deal.II library. I think that it surely needs to eliminate(or condense) the constraint lines after assembling global matrix. However, I could not find the relevant function for this one in deal.II docoment. The type of matrix is now "TrilinosWrappers::BlockSparseMatrix". COuld you help me? SUngho yoon Leeds UK _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
