配列の初期化がうまくいかない

Question

前提・実現したいこと

配列、pheromoneを初期化したい。

発生している問題・エラーメッセージ

initial_pheromoneでの
セグメントエラー

該当のソースコード

ファイル as.c
#include <stdio.h>
#include <math.h>
#include <limits.h>    
#include <assert.h>
#include <string.h>
#include <stdlib.h>
#include <time.h>

#include "sub.h"

double pheromone_0;
 double alpha;
 double beta;
 double rho;
 double q_0; 
 int n;
 int nant;

void erace(void){        

 alpha=1.0;
 beta=1.0;                       //ヒューリスティック情報の優先度
 rho=0.5;                         //フェロモン蒸発率
 q_0=0;

 pheromone_0=1.0;
printf("pheromone_0;%f\n",pheromone_0);
 initial_pheromone( pheromone_0);
printf("%f\n",pheromone[0][0]);
}

int main(int argc, char *argv[]){

 int i;

 n=51;

 nant=10;                         //アリの数


 start_program();

  erace();





 
 for(i=0;i<n;i++){
   free(pheromone[i]);
   free(total[i]);
 }
 free(total);
 free(pheromone);
   return(0);
}

sub.c

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <time.h>
#include <limits.h>

#include "tesb.h"

ant_struct *ant;      
ant_struct *bestant;

double **pheromone;
double **total;

void initial_pheromone(double initial_trail) //フェロモンの初期化

{
    long int i, j;

    for ( i = 0 ; i < n ; i++ ) {
    for ( j =0 ; j <=i ; j++ ) {
      pheromone[i][j] = initial_trail;
        pheromone[j][i]=pheromone[i][j];

    }
    }
}


void start_program(void)
{
  int i,j;
   if((ant=malloc(sizeof(ant_struct)*nant+sizeof(ant_struct *)*nant))==NULL){
   printf("out of memory,exit");
  exit(1);
 }

  for(i=0;i<nant;i++){
    ant[i].tour=calloc(n+1,sizeof(int));
    ant[i].visited=calloc(n,sizeof(char));
  }
  if((bestant=malloc(sizeof(ant_struct)))==NULL){
   printf("out of memory,exit");
  exit(1);
 }
  bestant->tour=calloc(n+1,sizeof(int));
  bestant->visited=calloc(n,sizeof(char));

 pheromone=malloc(sizeof(double *)*n);
  for(i=0;i<n;i++){
    pheromone[i]=malloc(sizeof(double*)*n);
  }

  total=malloc(sizeof(double *)*n);
  for(i=0;i<n;i++){
    total[i]=malloc(sizeof(double*)*n);
  }
}

ファイル　sub.h
typedef struct {
  int  *tour;
  char      *visited;
  int  tour_length;
} ant_struct;


extern ant_struct *ant;      
extern ant_struct *bestant;   

extern double   **pheromone;

extern int nant;


extern double rho;      //フェロモン蒸発率
extern double alpha;    //フェロモンの重要性
extern double beta;     //ヒューリスティックの重要性
extern double q_0;      //最良の道を選ぶ確率

int *best_ant_score;

void compute_total(void);



void initial_pheromone(double initial_trail);

void start_program(void);


extern double pheromone_0;

extern int n;

extern double **pheromone;
extern double **total;

試したこと

配列が問題なのかと思って配列を表示してみたりしたのですが配列のメモリもしっかり用意できていました。
修正してみたのですがまだエラーがでてしまいます。明日またmallocを使わない方法でやってみようと思います

Answer 1

double **pheromone;

pheromone=malloc(sizeof(int *)*n);

なんでdoubleなのにintポインタのサイズで確保している？
そこらへんがそもそもダメダメですね

Comments

[ro_tose]

確かにそこからおかしいですね。ご指摘ありがとうございます

[kazuma-s]

コードの意味は変ですが、double * も int * も char * も void * もポインタのサイズは同じなので悪影響はありません。

[y_waiwai]

おなじ「はず」だから、そんでいいとかやってたら、
> for(i=0;i<n;i++){
> pheromone[i]=malloc(sizeof(double*)*n);
> }
これもサイズが「おなじ」だからOKOKとするんでしょうか
まあ、そこらへんはうごいてんだからいいぢゃん、という考え方もありますがｗ

Answer 2

あなたの提示されたソースをコンパイルしたのですが、以下のエラーがでます。
(CentOs7.5 gcc バージョン 4.8.5)
gcc -o as as.c sub.c
as.c:29:1: エラー: 不明な型名 ‘iint’ です
iint main(int argc, char *argv[]){
^
as.c: 関数 ‘main’ 内:
as.c:39:15: エラー: ‘optimal’ が宣言されていません (この関数内での最初の使用)
scanf("%\n",&optimal);
^
as.c:39:15: 備考: 未宣言の識別子は出現した各関数内で一回のみ報告されます
as.c:56:9: エラー: ‘total’ が宣言されていません (この関数内での最初の使用)
free(total[i]);

エラーの発生しないソースを提示していただけませんでしょうか。

Comments

[ro_tose]

すいません、端折る場所が悪かったです。
今は外出先ですので他の方の指摘していただいたところを修正してまた提示させていただきます。

[tatsu99]

再度、コンパイルしましたが、エラーになりました。
sub.hはtesb.hにリネームしてコンパイルしました。
as.c:36: `nant' に対する定義されていない参照です
sub.c:20: `n' に対する定義されていない参照です
nant及びnをas.cの頭で定義するようにしました。
int nant;
int n;
このようにしてコンパイル＆実行すると、以下のようになりました。
as
1.000000

[ro_tose]

すいません、そこも間違っていました。その1.000のあとのinitial_pheromoneでセグメントエラーが出るのですがそちらでは起きずに実行できましたか？

[tatsu99]

はい。最後まで正常に終了しました。

[ro_tose]

ありがとうございます。こちらが変なことをしていてそれが原因でエラーが起きていそうなので確認させていただきます。

[ro_tose]

解決できました、ありがとうございます！

Answer 3

start_program()内のpheromone=malloc(sizeof(int *)*n);の時点での，nの値は想定した値になっていますか？

---

ところで，malloc()でメモリ確保するのはmustなのでしょうか？
とりあえず最初は各要素の数を固定した状態で作ってみる段階を経てみてはどうでしょうか．
例えば，

#define N 30
#define N_ANT 10

double pheromone[N][N];
ant_struct ant[N_ANT];

みたいな形．

Comments

[ro_tose]

malloc()についてまだよくわかっていないのでそちらで試して見るのもいいですかね。
他の方のアドバイスも試してダメでしたらそちらの方法も試してみます。

[ro_tose]

固定した値でやってみたのですがそれでもエラーが出たのでもしかしたら配列自体が用意されていないとかではなくて別の部分でエラーがでているんでしょうか・・・