Nach dem Aufruf einer Bibliotheksfunktion kann der gewünschte Zeigerwert nicht aufgerufen/geändert werden.

Question

Ich verwende eine externe Bibliothek https://people.sc.fsu.edu/~jburkardt/c_src/jacobi_eigenvalue/jacobi_eigenvalue.c, um den Eigenwert einer Matrix zu berechnen. Nach dem Aufruf der Bibliotheksfunktion habe ich jedoch einen Segmentierungsfehler erhalten, als ich versuchte, auf die gewünschten Zeigerwerte zuzugreifen.

void asphericity(particle *p, int N, float aveLTail[MAX2][MAX2], double *Nominator, double *Denominator, double *Asp ,double *cmx, double *cmy, double *cmz){ 
int i, j, k,m,n,M; 
double x[10000]; 
double y[10000]; 
double z[10000]; 
double *Q; 
double gyration[3][3]; 
double v[M*M]; 
double d[M]; 
float ut, up; 
float bs; 
double hdistance; 
double tdistance; 
double rog; 
int hcount=0; 
int tcount=0; 
int it_max; 
int it_num; 
int rot_num; 
double hx=0; 
double hy=0; 
double hz=0; 
double tx=0; 
double ty=0; 
double tz=0; 
m=0; 


it_max=100; 
for(i = 1; i<=N; i++) { 
    if (p[i].mol == p[i-1].mol){ 
     if (p[i-1].type==1 ||p[i-1].type==2 || p[i-1].type==3 || p[i-1].type==6){ 
      hcount++; 
      hx += p[i-1].x; 
      hy += p[i-1].y; 
      hz += p[i-1].z; 
     } 
     else if (p[i-1].type==4 && p[i-2].type==9){ 
      tcount++; 
      tx += p[i-1].x; 
      ty += p[i-1].y; 
      tz += p[i-1].z; 
      hcount += 2; 
      hx += p[i-2].x+p[i-3].x; 
      hy += p[i-2].y+p[i-3].y; 
      hz += p[i-2].z+p[i-3].z; 
     } 
     else if (p[i-1].type==4 && p[i-2].type!=9){ 
      tcount++; 
      tx += p[i-1].x; 
      ty += p[i-1].y; 
      tz += p[i-1].z; 
     } 

    } 
    else if (p[i].mol != p[i-1].mol && p[i-1].mol!=0) { 
      tcount++; 
      tx += p[i-1].x; 
      ty += p[i-1].y; 
      tz += p[i-1].z; 
      (hx) /= hcount; /*geometric center of head beads*/ 
      (hy) /= hcount; 
      (hz) /= hcount; 
      (tx) /= tcount; 
      (ty) /= tcount; 
      (tz) /= tcount; 
      hcount=0; 
      tcount=0; 
      hdistance=sqrt(pow((hx-*cmx),2)+pow((hy-*cmy),2)+pow((hz-*cmz),2)); 
      tdistance=sqrt(pow((tx-*cmx),2)+pow((ty-*cmy),2)+pow((tz-*cmz),2)); 
      if(hdistance>tdistance){ 
       x[m]=hx; 
       y[m]=hy; 
       z[m]=hz; 
       m++; 
      } 
      hx=0; 
      hy=0; 
      hz=0; 
      tx=0; 
      ty=0; 
      tz=0; 
     }    
    } 


for(j = 0; j<m; j++) { 
    gyration[0][0]+= pow((x[j]-*cmx),2); 
    gyration[0][1]+=(x[j]-*cmx)*(y[j]-*cmy); 
    gyration[0][2]+=(x[j]-*cmx)*(z[j]-*cmz); 
    gyration[1][0]+=(y[j]-*cmy)*(x[j]-*cmx); 
    gyration[1][1]+=pow((y[j]-*cmy),2); 
    gyration[1][2]+=(y[j]-*cmy)*(z[j]-*cmz); 
    gyration[2][0]+=(z[j]-*cmz)*(x[j]-*cmx); 
    gyration[2][1]+=(z[j]-*cmz)*(y[j]-*cmy); 
    gyration[2][2]+=pow((z[j]-*cmz),2); 
} 
for(i =0; i<3; i++){ 
    for(j=0; j<3; j++){ 
     gyration[i][j]/=m; 
     } 
} 

#define M 3 
double Gyration[M*M] = { 
      gyration[0][0],gyration[0][1],gyration[0][2], 
      gyration[1][0],gyration[1][1],gyration[1][2], 
      gyration[2][0],gyration[2][1],gyration[2][2] }; 
printf("%.6f\t%.6f\t%.6f\n", *Denominator, *Nominator, *Asp); 
jacobi_eigenvalue(M, Gyration, it_max, v, d, &it_num, &rot_num); 
printf("%.6f\t%.6f\t%.6f\n", *Denominator, *Nominator, *Asp); 
rog = sqrt(d[0]+d[1]+d[2]); 
*Denominator= 2*pow(rog,4); 
*Nominator= pow((d[2]-d[1]),2)+pow((d[2]-d[0]),2)+pow((d[1]-d[0]),2); 
* Asp = (*Nominator)/(*Denominator); 
# undef M 
} 

Denominator, Nominator und Asp wurden als Doppel Variablen in Hauptfunktion definiert und als Zeiger übergeben. Allerdings habe ich einen Segmentierungsfehler erhalten, als ich versuchte, nach dem Aufruf der Bibliotheksfunktion auf sie zuzugreifen. Ich konnte auf die Ausgabe v [] und d [] zugreifen.

Der erste Druck zeigt * Nenner, * Nominator und * Asp Null, während der Segmentierungsfehler beim zweiten Druck erscheint. Also ist etwas in der Bibliotheksfunktion nicht korrekt. Aber ich bin nicht sicher, wo genau das Problem ist und wie es zu beheben ist.

Vielen Dank im Voraus.

Answer 1

Eine Zusammenfassung Ihres Codes zeigt mehrere potenzielle "dividiere durch null" -Fehler. Die Konzentration auf * Nenner ist das folgende eine mögliche Gefahr:

* Asp = (*Nominator)/(*Denominator); 

Wenn * Nenner Null ist, dann werden Sie (Division durch Null) zum Absturz bringen. Außerdem haben Sie dies in einer for-Schleife:

int hcount=0; 

// some code 

if (p[i].mol == p[i-1].mol){ 

    if (p[i-1].type==1 ||p[i-1].type==2 || p[i-1].type==3 || p[i-1].type==6){ 

    hcount++; 

    } 
    else { 
    // hcount incremented 
    } 
    // some code 
} 
else if (p[i].mol != p[i-1].mol && p[i-1].mol!=0) { 
     tcount++; 
     tx += p[i-1].x; 
     ty += p[i-1].y; 
     tz += p[i-1].z; 
     (hx) /= hcount; /*geometric center of head beads*/ 
     (hy) /= hcount; 
     (hz) /= hcount; 
     (tx) /= tcount; 
     (ty) /= tcount; 
     (tz) /= tcount; 
    } 

Beachten Sie, dass in der anderen, teilen Sie durch HZÄHLUNG; hcount wird bei der Deklaration auf Null initialisiert und nur in der ersten if-Anweisung innerhalb der for-Schleife inkrementiert; Wenn diese erste Iteration diese erste if-Bedingung nicht erfüllt und ihren Weg in das else findet, dann teilen Sie durch null.

Auch nicht nötig

p[i].mol != p[i-1].mol 

in

else if (p[i].mol != p[i-1].mol && p[i-1].mol!=0) 

weil es folgt notwendigerweise, wenn der erste, wenn die Bedingung fehlschlägt.