Skip to content
Snippets Groups Projects
Commit 26c7b6e7 authored by Ruth Sabariego's avatar Ruth Sabariego
Browse files

normInfRightHandSide using abs ==> valid for double and complex

parent e07d69c0
No related branches found
No related tags found
No related merge requests found
...@@ -94,8 +94,8 @@ class linearSystemGmm : public linearSystem<scalar> { ...@@ -94,8 +94,8 @@ class linearSystemGmm : public linearSystem<scalar> {
double nor = 0.; double nor = 0.;
double temp; double temp;
for(unsigned int i=0;i<_b->size();i++){ for(unsigned int i=0;i<_b->size();i++){
temp = (*_b)[i]; temp = abs((*_b)[i]); // this is valid also for complex
if(temp<0) temp = -temp; //if(temp<0) temp = -temp;
if(nor<temp) nor=temp; if(nor<temp) nor=temp;
} }
return nor; return nor;
......
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