Changeset 368


Ignore:
Timestamp:
Jun 30, 2005 9:44:06 PM (14 years ago)
Author:
claird
Message:

Fixes in AddVectorQuotient? - the "a" value was not included in many branches of the code.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/LinAlgImpl/Serial/IpDenseVector.cpp

    r361 r368  
    993993        // then s is not homogeneous
    994994        for (Index i=0; i<Dim(); i++) {
    995           values_[i] = dense_z->scalar_ / values_s[i];
     995          values_[i] = a * dense_z->scalar_ / values_s[i];
    996996        }
    997997      }
     
    999999        // then z is not homogeneous
    10001000        for (Index i=0; i<Dim(); i++) {
    1001           values_[i] = values_z[i] / dense_s->scalar_;
     1001          values_[i] = values_z[i] * a / dense_s->scalar_;
    10021002        }
    10031003      }
    10041004      else {
    10051005        for (Index i=0; i<Dim(); i++) {
    1006           values_[i] = values_z[i] / values_s[i];
     1006          values_[i] = a * values_z[i] / values_s[i];
    10071007        }
    10081008      }
     
    10131013        // then s is not homogeneous
    10141014        for (Index i=0; i<Dim(); i++) {
    1015           values_[i] = val + dense_z->scalar_ / values_s[i];
     1015          values_[i] = val + a * dense_z->scalar_ / values_s[i];
    10161016        }
    10171017      }
     
    10191019        // then z is not homogeneous
    10201020        for (Index i=0; i<Dim(); i++) {
    1021           values_[i] = val + values_z[i] / dense_s->scalar_;
     1021          values_[i] = val + values_z[i] * a / dense_s->scalar_;
    10221022        }
    10231023      }
    10241024      else {
    10251025        for (Index i=0; i<Dim(); i++) {
    1026           values_[i] = val + values_z[i] / values_s[i];
     1026          values_[i] = val + a * values_z[i] / values_s[i];
    10271027        }
    10281028      }
     
    10331033        if (homogeneous_s) {
    10341034          for (Index i=0; i<Dim(); i++) {
    1035             values_[i] = c*values_[i] + dense_z->scalar_/dense_s->scalar_;
    1036           }
    1037         }
    1038         else {
    1039           for (Index i=0; i<Dim(); i++) {
    1040             values_[i] = c*values_[i] + dense_z->scalar_/values_s[i];
     1035            values_[i] = c*values_[i] + a * dense_z->scalar_/dense_s->scalar_;
     1036          }
     1037        }
     1038        else {
     1039          for (Index i=0; i<Dim(); i++) {
     1040            values_[i] = c*values_[i] + a * dense_z->scalar_/values_s[i];
    10411041          }
    10421042        }
     
    10451045        if (homogeneous_s) {
    10461046          for (Index i=0; i<Dim(); i++) {
    1047             values_[i] = c*values_[i] + values_z[i]/dense_s->scalar_;
    1048           }
    1049         }
    1050         else {
    1051           for (Index i=0; i<Dim(); i++) {
    1052             values_[i] = c*values_[i] + values_z[i]/values_s[i];
     1047            values_[i] = c*values_[i] + values_z[i] * a /dense_s->scalar_;
     1048          }
     1049        }
     1050        else {
     1051          for (Index i=0; i<Dim(); i++) {
     1052            values_[i] = c*values_[i] + a * values_z[i]/values_s[i];
    10531053          }
    10541054        }
Note: See TracChangeset for help on using the changeset viewer.