Skip to content
Snippets Groups Projects
Commit 37f8307e authored by Markus Blatt's avatar Markus Blatt
Browse files

Use norms correctly for unsymmetric dependency to make systems work.

Fixes flyspray #490

[[Imported from SVN: r987]]
parent a4993487
No related branches found
No related tags found
No related merge requests found
......@@ -1256,21 +1256,21 @@ namespace Dune
{
maxValue_ = std::min(- std::numeric_limits<typename Matrix::field_type>::max(), std::numeric_limits<typename Matrix::field_type>::min());
row_ = index;
diagonal_ = matrix_->operator[](row_)[row_];
diagonal_ = norm_(matrix_->operator[](row_)[row_]);
}
template<class M, class N>
inline void Dependency<M,N>::examine(const ColIter& col)
{
maxValue_ = std::max(maxValue_,
-*col);
-norm_(*col));
}
template<class M, class N>
template<class G>
inline void Dependency<M,N>::examine(G& graph, const typename G::EdgeIterator& edge, const ColIter& col)
{
if(-*col >= maxValue_ * alpha()) {
if(-norm_(*col) >= maxValue_ * alpha()) {
edge.properties().setDepends();
edge.properties().setInfluences();
}
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment