-2

このコードにセグメンテーション違反 11 があります (編集: 表示されなくなりました)。

デバッガーを入手して配列を整理するように言われたので、配列を 0 から修正するようにコードを編集しました。また、valgrind をインストールして実行しましたが、OSX で作業する場合は制限されます。gbd を実行しましたが、問題は以下のとおりです。

これが恐ろしく不吉に見える場合、または何らかの前進を遂げるためにコードが必要な場合は、私に知らせてください。フォーラムが期待する形式に質問を編集できます。

コードで gdb を実行すると、次のようになります。

Program received signal EXC_BAD_ACCESS, Could not access memory.
Reason: KERN_PROTECTION_FAILURE at address: 0x00007fff5e4c27c8
0x0000000100001048 in main () at TammesA3.cc:23
23    mpf_set_default_prec(1024);

コードは次のとおりで、後でループが開始される前に問題が発生します。

#include <stdio.h> 
#include <string.h>
#include <math.h>
#include <iostream>
#include <iomanip>
#include <fstream>
#include <time.h>
#include <stdlib.h>
#include <sstream>
#include <gmpxx.h>
using namespace std;
#define PI 3.14159265358979323846

mpf_t E[100000];

int main()
{

  int a,b,c,d,f,i,j,k,m,n,s,Success,Fails;
  double p,theta,phi,Time,Averagetime,Distance,Length,DotProdForce,
         Forcemagnitude,ForceMagnitude[200],Force[200][3],Epsilon[3],
         x[200][3],new_x[200][3],y[200][3],A[200],alpha[200][200],degree,bestalpha[500];
  unsigned long int t;

  mpf_set_default_prec(1024);

  mpf_t Energy,energy,Power,D,En[500],Ep,e;

  mpf_init(Power);
  mpf_init(D);
  mpf_init(Ep);
  mpf_init(e);
  mpf_init(Energy);
  mpf_init(energy);

  for(i=0;i<100000;i++){
   mpf_init(E[i]);
  }

  for(i=0;i<500;i++){
    mpf_init(En[i]);
  }

  mpf_set_d(e,0.00000000000000000000000000000000000000000000000000000000001);

  clock_t t1,t2;
  t1=clock();

t=1000;

while(t<1001){

n=20;

while(n<21){

cout << "N=" << n << "\n";

  b=1;
  Time=0.0;

  while(b<2){

    clock_t t3,t4;
    t3=clock();

    if(n>200){
      cout << n << " is too many points for me :-( \n";
      exit(0);
    }

    srand((unsigned)time(0));  

    for (i=0;i<n;i++){
      x[i][0]=((rand()*1.0)/(1.0*RAND_MAX)-0.5)*2.0;
      x[i][1]=((rand()*1.0)/(1.0*RAND_MAX)-0.5)*2.0;
      x[i][2]=((rand()*1.0)/(1.0*RAND_MAX)-0.5)*2.0;

      Length=sqrt(pow(x[i][0],2)+pow(x[i][1],2)+pow(x[i][2],2));

      for (k=0;k<3;k++){
        x[i][k]=x[i][k]/Length;
      }
    }

    for(i=0;i<n;i++){
      for(j=0;j<3;j++){
        cout << "x[" << i << "][" << j << "]=" << x[i][j] << "\n";
      }
    }

    /* Points distributed randomly and normalised so they sit on unit sphere */

    mpf_set(Energy,0);

    for(i=0;i<n;i++){
      for(j=i+1;j<n;j++){
        Distance=sqrt(pow(x[i][0]-x[j][0],2)+pow(x[i][1]-x[j][1],2)
                     +pow(x[i][2]-x[j][2],2));

        mpf_set_d(D,Distance);

        mpf_pow_ui(Power,D,t);      
        mpf_ui_div(Power,1.0,Power);
        mpf_add(Energy,Energy,Power);

      }
    } etc.

この後、私のコードはループに入りますが、コメントで述べたコードは非常に長いため、要求されない限り再度挿入することはありません。

4

1 に答える 1