c++のDebug Assertion Failedについて

Question

前提・実現したいこと

ファイルから座標を読み込んで、凸包を作りたいです。

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

Expression:stream!=nullptr

For information on how your program can cause an assertion failure,see the Visual C++ documentation on asserts

該当のソースコード

#define MAX_NUMBER 1000
#define FLAG_OFF 0
#define FLAG_ON 1
#define EPS 0.000001
#define V_NO_ERROR 0
#define V_ERROR 1

typedef struct {
	double x, y;
	int flag;
}POINT;

POINT p[MAX_NUMBER];
int number_of_points;

int find_point_with_minimum_x(void);

void make_2d_convex_hull(void);
void take_coordinate_of_point(int point_number, double v[]);
void determine_direction_of_base_line(double uv[], int point1, int point2);
void sort_using_distance(int ct[], int n, int p);
void vsub_(double v[], double v1[], double v2[], int *n);
void vunit_(double v[], double v0[], int *n, int *ierror);
double vinpro_(double v1[], double v2[], int *n);

int main(void)
{
	FILE *fp;
	int i = 0;
	int point_number;
	int n = 2;
	double v1[2], v2[2], v[2];
	double d;
	errno_t err;
	
	/*点データの入力*/
	if ((err = fopen_s(&fp,"point_table.dat.txt", "r")) == '0') {
		printf("cannot open point_table.dat\n");
		return 0;
	}
	while (fscanf_s(fp, "%lf%lf", &p[i].x, &p[i].y) != EOF) {
			p[i].flag = FLAG_OFF;
			i++;
		}
		fclose(fp);
	/*点の数の設定*/
	number_of_points = i;

	/*点データの出力*/
	for (i = 0; i < number_of_points; i++) {
		printf("%3d %7.3lf %7.3lf %d\n", i, p[i].x, p[i].y, p[i].flag);
	}
	/*２次元コンベックスハルの生成*/
	make_2d_convex_hull();

	return 0;
}

void make_2d_convex_hull(void)
{
	int start_point;
	int convex_point;
	int number_of_convex_points = 0;
	int cp[MAX_NUMBER];/*凸多角形を構成する点の並び*/
	int ctemp[10];
	int count;
	int n = 2;
	int ierror;
	int i, k;
	double v0[2], v1[2], v2[2], v[2], uv[2];
	double d, dmax;

	/*座標が最小の点を求める*/
	convex_point = find_point_with_minimum_x();
	printf("凸多角形を構成する１番目の点：%d\n", convex_point);
	start_point = convex_point;
	cp[number_of_convex_points] = convex_point;
	number_of_convex_points++;

	/*基準線の方向を設定する（ｙ軸の負の方向）*/
	uv[0] = 0.0; uv[1] = -1.0;
	/*凸多角形を構成する２番目以降の点を求める*/
	while (1) {
		take_coordinate_of_point(convex_point, v1);
		k = convex_point;
		dmax = -1.0;
		for (i = 0; i < number_of_points; i++) {
			if (i != convex_point || p[i].flag == FLAG_OFF) {
				take_coordinate_of_point(i, v2);
				vsub_(v0, v2, v1, &n);
				vunit_(v, v0, &n, &ierror);
				d = vinpro_(uv, v, &n);
				if (fabs(dmax - d) < EPS) {
					ctemp[count] = i;
					count++;
				}
				else if (dmax < d) {
					dmax = d;
					k = i;
					ctemp[0] = i;
					count = 1;
				}
			}
		}
		if (count >= 2) {
			/*距離を用いて点を並び替える*/
			sort_using_distance(ctemp, count, convex_point);
			for (i = 0; i < count; i++) {
				convex_point = ctemp[i];
				p[convex_point].flag = FLAG_ON;
				cp[number_of_convex_points] = convex_point;
				number_of_convex_points++;
				printf("凸多角形を構成する%d番目の点：%d\n", number_of_convex_points, convex_point);
			}
		}
		else {
			convex_point = k;
			p[convex_point].flag = FLAG_ON;
			cp[number_of_convex_points] = convex_point;
			number_of_convex_points++;
			printf("凸多角形を構成する%d番目の点：%d\n", number_of_convex_points, convex_point);
		}
		/*最初の点の戻れば終了*/
		if (convex_point == start_point)break;
		/*基準線の方向を設定する*/
		determine_direction_of_base_line(uv, cp[number_of_convex_points - 2], convex_point);
	}
}

void take_coordinate_of_point(int point_number, double v[])
{
	v[0] = p[point_number].x;
	v[1] = p[point_number].y;

}

int find_point_with_minimum_x(void)
{
	int point_number;
	int i;
	double min;

	point_number = 0;
	min = p[0].x;
	for (i = 1; i < number_of_points; i++) {
		if (min > p[i].x) {
			min = p[i].x;
			point_number = i;
		}
	}
	return point_number;
}

void determine_direction_of_base_line(double uv[], int point1, int point2)
{
	double v1[2], v2[2], v[2];
	int n = 2;
	int error;

	take_coordinate_of_point(point1, v1);
	take_coordinate_of_point(point2, v2);
	vsub_(v, v2, v1, &n);
	vunit_(uv, v, &n, &error);
}
void sort_using_distance(int ct[], int n ,int p)
{
	int i, j, k;
	double v1[2], v2[2];
	double dist[10];
	double d;
	double min, temp1;
	int temp2;

	/*２点間の距離の計算*/
	take_coordinate_of_point(p, v1);
	for (i = 0; i < n; i++) {
		take_coordinate_of_point(ct[i], v2);
		d = 0.0;
		for (j = 0; j < 2; j++) {
			d += (v2[j] - v1[j])*(v2[j] - v1[j]);
		}
		dist[i] = sqrt(d);
	}
	/*距離の昇順に点を並び替える*/
	for (i = 0; i < n - 1; i++) {
		min = dist[i];
		k = i;
		for (j = i + 1; j < n; j++) {
			if (dist[j] < min) {
				min = dist[j];
				k = j;
			}
		}
		temp1 = dist[i];
		dist[i] = dist[k];
		dist[k] = temp1;
		temp2 = ct[i];
		ct[i] = ct[k];
		ct[k] = temp2;
	}
}

/*単位ベクトル*/
void vunit_(double v[], double v0[], int *n, int *ierror)
{
	int i;
	double d = 0.0;
	*ierror = V_NO_ERROR;
	for (i = 0; i < *n; i++)d += v0[i] * v0[i];
	if (d < 0.00000001) {
		*ierror = V_ERROR;
		return;
	}
	d = sqrt(d);
	for (i = 0; i < *n; i++) v[i] = v0[i] / d;
	return;
}

/*ベクトルの差*/
void vsub_(double v[], double v1[], double v2[], int *n)
{
	int i;
	for (i = 0; i < *n; i++)v[i] = v1[i] - v2[i];
	return;
}

/*内積*/
double vinpro_(double v1[], double v2[], int *n)
{
	int i;
	double t = 0.0;
	for (i = 0; i < *n; i++) t+= v1[i] * v2[i];
	return t;
}

試したこと

ここに問題に対して試したことを記載してください。

補足情報（FW/ツールのバージョンなど）

このエラーメッセージについて教えていただきたいです。

Comments

[episteme]

何の説明もなされていない。やりなおし。

Answer 1

Expression:stream!=nullptr

コード中(ライブラリの中かも)で "stream が nullptr であってはいけない"ってassertに引っかかって
停止しています。心当たりはありませんか?

Comments

[namekuhito]

返信ありがとうございます。
ファイルから座標のデータを取得しているのですが、fscanf_sで使い方がいまいち理解できていないところがあり、その点が思い当たります。
エラーの内容はstreamが空であるということですか？

[episteme]

です。session は　assert に渡された変数名かなにかでしょう。

[namekuhito]

返信が遅くなってしまいすみません。
それは、fscanf_sがうまく動いていないからでしょうか？

[episteme]

おそらくそうでしょう。
fpがNULLなんじゃありませんか? であるなら、fopen_sでファイルを開けそこねていると思われます。

[namekuhito]

分かりました。
今、少し変えて以下のようにしました。
/*点データの入力*/
err = fopen_s(&fp, "a.txt", "r");
if (err == 0) {
printf("cannot open point_table.dat\n");
return 0;
}
else {
while (fscanf_s(fp, "%lf%lf", &p[i].x, &p[i].y) != EOF) {
p[i].flag = FLAG_OFF;
i++;
}
fclose(fp);
}

[namekuhito]

そして、他の方が貼っていただいたURLのように
/*点データの入力*/
err = fopen_s(&fp, "a.txt", "r");
if (err == 0) {
printf("cannot open point_table.dat\n");
return 0;
}
ここの、if文は'=='ではなく、'!='に変更するべきでしょうか？

[episteme]

試してみるのにいちいち許可もらうんですか?
# fopen_s は成功したら0を返します。

[namekuhito]

すみません
試してみたのですが、printf("cannot open point_table.dat\n");が表示されて、elseにはいってないようでした。

[episteme]

> if ((err = fopen_s(&fp,"point_table.dat.txt", "r")) == '0') {

オープンしたいファイルは　"point_table.dat.txt"　ですか? そんなファイルあるんですか?

[namekuhito]

すみません
/*点データの入力*/
err = fopen_s(&fp, "a.txt", "r");
if (err == 0) {
printf("cannot open a.txt\n");
return 0;
}
else {
while (fscanf_s(fp, "%lf %lf", &p[i].x, &p[i].y) != EOF) {
p[i].flag = FLAG_OFF;
i++;
}
fclose(fp);
}
に書き直しました。
テキストファイルは、プログラムと同じ場所にあります。

[episteme]

fopen_s は成功したら0を返します。

[namekuhito]

ありがとうございます。プログラムの関数の動作がよく分かってなかったみたいです。
そこで、テキストファイルを変更することで動作しました。
長い間、ありがとうございました。

Answer 2

if ((err = fopen_s(&fp,"point_table.dat.txt", "r")) == '0')

↓
if ((err = fopen_s(&fp,"point_table.dat.txt", "r")) != 0)
参考：fopen_sをfopenの代わりに使うときのやりかた

Comments

[cateye]

たぶん、fopen_s()が正常に終わっていない。
Windowsで"point_table.dat.txt"がファイル名として有効かどうか不明？

[namekuhito]

返信が遅くなってしまいすみません。
ファイル名の有効かどうかを確認する方法が分かりません。すみません。
point_table.datのファイル名を変更して、a.txtとかにしたほうが良いでしょうか？

[episteme]

openしたいファイルは point_table.dat.txt ですか? そんなファイルあるんですか? コードにはそう書かれてますよ?