-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

1

剩下的问题（1）是（三行）

mpf_set(Energy,0);

它的原型是

void mpf_set (mpf_t ROP, mpf_t OP)

anmpf_t是一个数组，

typedef __mpf_struct mpf_t[1];

所以你传递一个空指针来读取值。取消引用空指针是导致分段错误的另一个典型原因（除了访问超出数组边界）。

您打算使用mpf_set_d,mpf_set_ui或mpf_set_si那里（您使用哪一个并不重要）。

(1)我从修订版 2 中复制了代码并将两个大数组移动E到alpha文件范围以减轻堆栈的负担，然后mpf_set(variable,0)在第 94、147 和 271 行修复三个，它运行完成。

这并不一定意味着代码是正确的，在我阅读的输出末尾

Successes=0 Failures=1
Success Rate=0.00000

这暗示了代码逻辑中可能存在的问题，但我不知道期望什么输出。

于 2013-04-02T11:28:24.480 回答