поиск нуля в матрице

[ciaonataha]

Код какой был, такой и остался. Я совершенно ничего не меняла. Только почему то в out.тхт выводит для векторов индексов 0000
+
А я хотела бы в главной фукции распечативать эти вектора. Но , когда я их распечтиваю они аффишируют индексы и элементы матрицы, которые не равны 0 и дозопалняутся нулями из-за их первоночальных размеров.

  Array <double,1> L1(IMAX*JMAX*KMAX*FMAX);                     
  Array <double,1> iL1(IMAX*JMAX*KMAX*FMAX);                    
  Array <double,1> jL1(IMAX*JMAX*KMAX*FMAX);                    
  Array <double,1> kL1(IMAX*JMAX*KMAX*FMAX);                    
  Array <double,1> lL1(IMAX*JMAX*KMAX*FMAX);

Вот мне и надо поменять их размеры в функции поиска не нулевых элементов

[Sazary]

Код какой был, такой и остался. Я совершенно ничего не меняла. Только почему то в out.тхт выводит для векторов индексов 0000

Ну если бы вы совсем ничего не меняли, то все бы работало.

А я хотела бы в главной фукции распечативать эти вектора.

Показываю на своем коде (из прошлой темы):

for(i=0;i<N1;i++)
  for(j=0;j<N2;j++)
   for(k=0;k<N3;k++)
    for(l=0;l<N4;l++)
     if(MAT[i][j][k][l]!=0)
      {
       V1[vk] = i;
       V2[vk] = j;
       V3[vk] = k;
       V4[vk] = l;
       Vz[vk++] = MAT[i][j][k][l];
      }

Это основной кусок, который и заполняет вектора. А vk, в конечном итоге, будет содержать реальный размер каждого вектора (количество ячеек, заполненных нужными числами).

В вашем случае, как я понимаю, это vel.

[ciaonataha]

Все в точности тоге самое и у меня написану. Тем не менее проблема остается.

 int i1 = Afinal.extent(0);
 int j1 = Afinal.extent(1);
 int f1 = Afinal.extent(2);
 int l1 = Afinal.extent(3);
 int i, j , f, l;
 int vel = 0;
    
for(i = 0; i < i1; i++)
 for(l = 0; l < l1; l++)
  for (j = 0; j < j1; j++)
    for (f = 0; f < f1; f++)
    if (Afinal(i,j,f,l)!=0)
     {
      iA(vel)=i;
      jA(vel)=j;       
      kA(vel)=f;
      lA(vel)=l;   
      A(vel++)=Afinal(i,j,f,l); 
     }

[Sazary]

А вот это:

int i1 = Afinal.extent(0);
 int j1 = Afinal.extent(1);
 int f1 = Afinal.extent(2);
 int l1 = Afinal.extent(3);

что значит?
Может, проблема при чтении массива?