проблема с gcc

[d9m0n]

у меня возникла проблема с программой. есть код (кусок):

PHP код:

for(int k=0;k<el;k++)
{
    if(k!=j)
    vects.push_back(newplane[k]-newplane[j]);
}
//cout<<"ab:"<<vects[0]<<endl;
//cout<<"ac:"<<vects[1]<<endl;
cout<<"ap:"<<vec<<endl;
double** matr = new double*[el-1];    
for (int ii=0; ii<el-1; ii++)
{
    matr[ii] = new double[el-1];
} 


здесь была необъяснимая проблема со строкой:
cout<<"ap:"<<vec<<endl;
Если эта строка есть, то выполнялось все нормально, если же её закомментировать, то ход выполнения программы изменялся.
Этот код располагается во вложенном цикле. Если убрать эту строку,то после break из этого цикла выполнение передавалось на начало внешнего цикла, но номер итерации внешнего цикла не увеличивался!!, хотя в заголовке цикла прописано увеличение счетчика после каждой итерации (for(int i=0;i<n;i++)).
Сначала я с этим мирился, но теперь, после внесения некоторых изменений в код, даже с этой строкой код просто зацикливается на одном и том же значении счетчика внешнего цикла.

Может быть, кто сталкивался с подобной проблемой?
Компилятор gcc 4.2.1, система - gentoo 2.6.27

[ISergeyN]

Весь код покажи..

[d9m0n]

вот код функции:

PHP код:

vector<point> qhull::extremums()
{
    std::vector <double> M(7);
    M[0]=12.011;//"c"
    M[1]=51.996;//"cr"
    M[2]=95.940;//"mo"
    M[3]=58.933;//"co"
    M[4]=26.982;//"al"
    M[5]=58.690;//"ni"
    M[6]=55.847;//"fe"
    std::vector <double> CONCENTRATION(7);
    CONCENTRATION[0]=2.e-3;   //"c"
    CONCENTRATION[1]=2.5e-2;  //"cr"
    CONCENTRATION[2]=1.4e-2;  //"mo"
    CONCENTRATION[3]=10e-2;   //"co"
    CONCENTRATION[4]=0.9e-2;  //"al"
    CONCENTRATION[5]=14e-2;   //"ni"
    CONCENTRATION[6]=1.-CONCENTRATION[0]-CONCENTRATION[1]-CONCENTRATION[2]-CONCENTRATION[3]-CONCENTRATION[4]-CONCENTRATION[5];//"fe"

    std::vector <double> CONCENTRATION_MOLE(7);
    CONCENTRATION_MOLE=CONCENTRATION_MASStoCONCENTRATION_MOLE(CONCENTRATION,M);//
    vector<double> coordpt;
    for(int i=0;i<el-1;i++) coordpt.push_back(CONCENTRATION_MOLE[i]);
    point pt(coordpt);
    vector<point> pts;
    int algoritm=0;
    if(el==2)
    {
    for(int i=0;i<qh.size();i++)
    {
        if(qh[i][0].function!=qh[i][1].function)
        {
        pts.clear();
        pts.push_back(qh[i][0]);
        pts.push_back(qh[i][1]);
        return pts;
        }
    }
    }
    else
    {
    hyperplane plane;
    int num;
    bool second=false;
    for(int i=0;i<qh.size();i++)
    {
        num=0;
        for(int j=0;j<el;j++)
        {
        if(qh[i][j].function!=qh[i][j+1].function)
        {
            algoritm=qh[i][j+1].function;
            num++;
            continue;
        }
        }
        if(num!=0)
        {
                 plane=qh[i];
        double angle=0;
        double temp;
        vector<double> newpoint;
                     std::vector <double> m(2);
        vector<point> pts1;
        for(int j=0;j<el;j++)
        {
            newpoint.clear();
               m[0]=1.0;
               m[1]=algoritm?3.0:1.0;
               newpoint.push_back(plane[j][0]*m[1]/(m[0]+m[1]-m[1]*(1.0-plane[j][0])));
               for (i = 1; i<el-1; i++)
                newpoint.push_back(plane[j][i]*m[0]/(m[0]+m[1]-m[1]*(1.0-plane[j][0])));
            pts1.push_back(point(newpoint,plane[j].function));
        }
        hyperplane newplane(pts1);
        for(int j=0;j<el;j++)
        {
            double M;
            vector<double> b_mas;
            vector<point> vects;
            vector<double> coef;
            point vec;
            vec=(pt-newplane[j]);
            for(int k=0;k<el;k++)
            {
            if(k!=j)
                vects.push_back(newplane[k]-newplane[j]);
            }
            //cout<<"ab:"<<vects[0]<<endl;
            //cout<<"ac:"<<vects[1]<<endl;
            cout<<"ap:"<<vec<<endl; // все дело в этой строке
                double** matr = new double*[el-1];
                for (int ii=0; ii<el-1; ii++)
                {
                matr[ii] = new double[el-1];
                }
            for(int j=0;j<el-1;j++)
            for(int k=0;k<el-1;k++)
            {
                double elem=0;
                for(int u=0;u<el-1;u++)
                elem+=vects[j][u]*vects[k][u];
                matr[j][k]=elem;
            }
            M=hyperplane::det(matr,el-1);
            for (int j=0; j<el-1; j++) 
                delete[] matr[j];
                delete[] matr;
            for (int j=0; j<el-1; j++)
            {
            double elem=0;
            for(int k=0;k<el-1;k++)
                elem+=vects[j][k]*vec[k];
            b_mas.push_back(elem);
            }
            for(int j=0;j<el-1;j++)
            coef.push_back(b_mas[j]/M);
            double coef_sum=0;
            bool lower=false;
            for(int j=0;j<el-1;j++)
            {
            if(coef[j]<=0.0)
            {
                lower=true;
                break;
            }
            coef_sum+=coef[j];
            }
            cout<<"test u and v"<<endl;
            if(lower||(coef_sum>1))break;
            cout<<"this is extremum place"<<endl;
            pts.clear();
            vector<point> line;
            for(int k=0;k<el;k++)
            {
            vector<double> tpt;
            if(newplane[k].function!=algoritm)
            {
                for(int u=0;u<el-1;u++)
                tpt.push_back(newplane[k][u]);
                pts.push_back(tpt);
            }
            else
            {
                line.clear();
                for(int u=0;u<el-1;u++)
                tpt.push_back(newplane[k][u]);
                line.push_back(tpt);
            }
            }
            line.push_back(pt);
            if(el==3) // dlya 3d
            {
            vector<double> coeffs1=hyperplane(pts).getCoeffs();
            vector<double> coeffs2=hyperplane(line).getCoeffs();
            double t=0;
            for(int k=0;k<el-1;k++)
                t+=coeffs1[k]*line[0][k];
            t+=coeffs1[el-1];
            point vect=line[0]-line[1];
            t=t/(coeffs1[0]*vect[0]+coeffs1[1]*vect[1]+coeffs1[2]);
            vector<double> newpt1;
            for(int k=0;k<el-1;k++)
                newpt1.push_back(line[0][k]+vect[k]*t);
            vector<double> newpt;
            pts.clear();
            pts.push_back(point(newpt1));
            pts.push_back(line[0]);
            return pts;
            }
        }
        }    
    }
    }
    return pts;
} 


[ISergeyN]

cout<<"ap:"<<vec<<endl; // все дело в этой строке

Значит проблема в этом операторе

std::ostream &operator << (std::ostream &lhs, const point &rhs){
что здесь?
}

[d9m0n]

я думаю, что с оператором все нормально, я его проверил, оттестировал. Проблем с использованием не возникало. вот его код:

PHP код:

ostream& operator<<(std::ostream &os, const point pt)
{
    os<<"(";
    for (int i=0;i<pt.dimensions();i++)// пока меньше размерности пространства
    {
        os<<(double)pt[i]; // вывод координаты
        if (i!=pt.dimensions()-1)// вывод ";" между координатами
            os<<";";
    }
    os<<"("<<pt.function<<")";
    os<<")";
    return os;
} 


operator[] в классе point переопределен. Возвращает значение из массива внутри класса по индексу.

[ISergeyN]

Здесь точно нехватает &

ostream& operator<<(std::ostream &os, const point &pt)

Если в point нету оператора копирования (и конструктора копирования) то проблема может быть из за этого..

[d9m0n]

проверил, не помогает.
Как в этом может быть проблема? оператор вывода ничего не меняет в объекте. И не может поменять, т.е. объект передается по значению.
В классе нет указателей и виртуальных функций, значит конструктор копирования по умолчанию вполне подходит.
Раньше у меня была проблема только при отсутствии строки вывода значения, сейчас же - и с ней те же самые проблемы. Отладчик как то не очень корректно переходит при выходе из цикла по break. Сначала показывает выполнение следующего за break оператора, а затем переходит к началу внешнего цикла. Я думаю, что проблема с компилятором. Попробую обновить его.