c - 固有ベクトルを計算するために dstev を使用してゼロ値を取得しています

Question

私は gcc を使用して mac os x でコンパイルしています。Intel の mkl_lapack.h ライブラリがインストールされています。プログラムには NxN の三重対角行列があるため、行列の値を格納するために 2 つのベクトルを使用するだけです。「d」ベクトルは主対角で、副対角の値は「e」に格納されます。まず最初に値を初期化し、次に行列が 16x16 であるため (入力では値 16 を argv[1] として与えている)、ベクトルを 2 つのベクトルに分割します (dstev を一度だけ使用することもできますが、それは実験目的)、d[0] から d[N/2-1] まで d[N/2] から d[N-1] までの最初のベクトルがあります。したがって、「e」と「d」の値を初期化したら、dstev を 2 回呼び出します。しかし、dstev を 2 回呼び出した後、すべての「z」ベクトルで、8x8 のゼロ以外の値の 2 つの部分行列のみを持つ必要があります。しかし、「z」を入力しようとすると、一部の値が 0.0 になり、なぜこれが起こるのか説明できません。

#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#include <mpi.h>
#include "mkl_lapack.h"

int main(int argc, char **argv)
{
    int N,dim,info;
    double *d,*e,*z,*work;
    char jobz='V';
    switch(argc)
    {
        case 2:
            N=atoi(argv[1]);
            break;
        default:
            fprintf(stderr,"Errore nell' inserimento degli argomenti\n");
            exit(EXIT_FAILURE);
            break;
    }
    if(N%2!=0)
    {
        fprintf(stderr,"La dimensione della matrice deve essere pari\n");
        exit(EXIT_FAILURE);
    }
    dim=N/2;
    d=(double*)malloc(N*sizeof(double));
    e=(double*)malloc((N-1)*sizeof(double));
    z=(double*)malloc(N*N/2*sizeof(double));
    work=(double*)malloc((N-1)*2*sizeof(double));
    for(int i=0;i<N-1;i++)
    {
        d[i]=(double)(i+3);
        e[i]=1.0;
    }
    dstev(&jobz,&dim,d,e,z,&dim,work,&info);
    dim--;
    dstev(&jobz,&dim,&d[N/2],&e[N/2],&z[N*N/4],&dim,&work[N-1],&info);
    for(int i=0;i<(N*N/2);i++)
        printf("(%f) ",z[i]);
    return 0;
}

このことを明確に説明したことを願っています。不明な点がある場合はお知らせください。