>From tutorial 12, there is such a funtion called when assembling,
for (unsigned int i=0; i<dofs_per_cell; ++i)
for (unsigned int k=0; k<dofs_per_cell; ++k)
system_matrix.add(dofs[i], dofs_neighbor[k],
ue_vi_matrix(i,k));
"dofs[i]" and "dofs_neighbor[k]" both can be have a different elements of
each arrays.
So, I would change this funtion into
"ConstraintMatrix::distribute_local_to_global"
For main reason, I need to have a periodic constrain matrix for dg system.
A continous garekin system is fine at this momonent.
On Mon, August 3, 2009 4:02 pm, Wolfgang Bangerth wrote:
>
>> void ConstraintMatrix::distribute_local_to_global
>> ( const FullMatrix< double > & local_matrix,
>> const Vector< double > & local_vector,
>> const std::vector< unsigned int > & row_indices,
>> const std::vector< unsigned int > & col_indices,
>> MatrixType & global_matrix,
>> VectorType & global_vector
>> ) const [inline]
>>
>>
>> "const std::vector< unsigned int > & col_indices" is hopefully added.
>
> Commonly, the row and column indices are the same. How come this is not
> the
> case for you?
>
> Best
> W.
>
> --
> -------------------------------------------------------------------------
> Wolfgang Bangerth email: [email protected]
> www: http://www.math.tamu.edu/~bangerth/
>
_______________________________________________
dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
