座標の値の重複を削除したいです C言語

Question

前提・実現したいこと

C言語で2つのdatファイルを読み込み、結合して新しいファイルに出力するプログラムを作ってます。2つのdatファイルには、節点のx座標、y座標、z座標があり、それぞれのファイルから値を二次元配列に保存するところまではできたのですが、新しいファイルに出力する際に、重複をなくして出力しないといけないのですが、重複をなくしてすっきりと表示できません。
二次元配列には　接点1のx座標は[0][0]、y座標は[0][1]、z座標は[0][2]、接点2のx座標は[1][0]、y座標は[1][1]、z座標は[1][2]、のように保存してあります。

今回が最初の質問なので伝わりにくかったらごめんなさい。お願いします。
コードもきれいとは言えないので改善点もあればお願いします。

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

for文とif文を使ったのですが、重複を取り除いた後の座標の数がとても多くなってしまいました。
きれいに取り除けず、取り除いた値が重複してしまいます。

該当のソースコード

#include <stdio.h>

int main(){
int point1[1024][16];
int point2[1024][16];
double epoint1[1024][2];
double epoint2[1024][2];
double ecomb2[1024][2];

/* ファイル読み込みなどは割愛 */


//読み取って配列に保存
for(int i = 1; i < point1[0][0]+1; i++){
  //point1[0][1]、point2[0][1]はヘッダー部分に記された接点の数
  //すでに配列に保存済み。
    for(int j = 0; j < 3; j++){
      fscanf(readf_in1,"%lf",  //lfは倍精度実数で読み取る
        &(epoint1[i][j]),
        &(epoint1[i][j+1]),
        &(epoint1[i][j+2])
      );
    }
  }
  for(int i = 1; i < point2[0][0]+1; i++){//配列に保存
    for(int j = 0; j < 3; j++){
      fscanf(readf_in2,"%lf",  //lfは倍精度実数で読み取る
        &(epoint2[i][j]),
        &(epoint2[i][j+1]),
        &(epoint2[i][j+2])
      );
    }
  }

/* -- ここからわかりません -- */
  int k = 0;
  for(int i = 0; i < point2[0][0]; i++){
    for(int j =0; j < point1[0][0]; j++){
      if(epoint1[i][0] == epoint2[j][0] && epoint1[i][1] == epoint2[j][1] && 
epoint1[i][3] == epoint2[j][3] ){
          printf("match\n");//確認用
      }
      else{
        k = k + 1;
        ecomb2[k][0] = epoint2[j][0];
        ecomb2[k][1] = epoint2[j][1];
        ecomb2[k][2] = epoint2[j][2];
      }
    }
  }
/* -- -- */
  //確認用
  for(int i = 0; i < k; i++){
   fprintf(writef_in,"%24.13E%24.13E%24.13E\n",
   ecomb2[i][0], ecomb2[i][1], ecomb2[i][2]
  );
}

Comments

[kazuma-s]

2つのファイルそれぞれには重複はないのですか？

[退会済みユーザー]

重複はないです。

Answer 1

出力したものだけを憶えておいて、
入力の値がそこにあるかどうかをチェックすればよいでしょう。

#include <stdio.h>
#include <stdlib.h>  // exit

int find_add(double e[1024][3], int n, const double *p)
{
	if (n >= 1024) { puts("too many data"); exit(1); }
	for (int i = 0; i < 3; i++) e[n][i] = p[i];
	for (int i = 0; i < n; i++)
		if (e[i][0] == p[0] && e[i][1] == p[1] && e[i][2] == p[2]) return 0;
	return 1;
}

int main(void)
{
	FILE *f1 = fopen("f1.txt", "r");
	if (!f1) return 1;
	FILE *f2 = fopen("f2.txt", "r");
	if (!f2) return 2;
	FILE *f3 = fopen("f3.txt", "w");
	if (!f3) return 3;
	int point1[1][1] = { 4 }, point2[1][1] = { 3 };
	double epoint[1024][3], p[3];
	int n = 0;
	while (fscanf(f1, "%lf%lf%lf", p, p+1, p+2) == 3)
		if (find_add(epoint, n, p)) {
			fprintf(f3, "%24.13E%24.13E%24.13E\n", p[0], p[1], p[2]);
			n++;
		}
	while (fscanf(f2, "%lf%lf%lf", p, p+1, p+2) == 3)
		if (find_add(epoint, n, p)) {
			fprintf(f3, "%24.13E%24.13E%24.13E\n", p[0], p[1], p[2]);
			n++;
		}
	fclose(f1), fclose(f2), fclose(f3);
}

f1.txt

44.444 44.444 44.444
33.333 33.333 33.333
22.222 22.222 22.222
11.111 11.111 11.111

f2.txt

55.555 55.555 55.555
33.333 33.333 33.333
66.666 66.666 66.666

実行後の f3.txt

     4.4444000000000E+01     4.4444000000000E+01     4.4444000000000E+01
     3.3333000000000E+01     3.3333000000000E+01     3.3333000000000E+01
     2.2222000000000E+01     2.2222000000000E+01     2.2222000000000E+01
     1.1111000000000E+01     1.1111000000000E+01     1.1111000000000E+01
     5.5555000000000E+01     5.5555000000000E+01     5.5555000000000E+01
     6.6666000000000E+01     6.6666000000000E+01     6.6666000000000E+01

追記
それぞれのファイルに重複がないのなら、
1つめのファイルの内容を全部覚えておいて、
2つめのファイルの値がそこにあるかどうかをチェックすればよいでしょう。

#include <stdio.h>

int main(void)
{
	FILE *f1 = fopen("f1.txt", "r");
	if (!f1) return 1;
	FILE *f2 = fopen("f2.txt", "r");
	if (!f2) return 2;
	FILE *f3 = fopen("f3.txt", "w");
	if (!f3) return 3;
	int point1[1][1] = { 4 }, point2[1][1] = { 3 };
	double ep[1024][3], p[3];
	int n = 0, i;
	for (; n < 1024 && fscanf(f1, "%lf%lf%lf", ep[n], ep[n]+1, ep[n]+2) == 3; n++)
		fprintf(f3, "%24.13E%24.13E%24.13E\n", ep[n][0], ep[n][1], ep[n][2]);
	if (n == 1024) { puts("too many data"); return 4; }
		
	while (fscanf(f2, "%lf%lf%lf", p, p+1, p+2) == 3) {
		for (i = 0; i < n; i++)
			if (ep[i][0] == p[0] && ep[i][1] == p[1] && ep[i][2] == p[2]) break;
		if (i == n)
			fprintf(f3, "%24.13E%24.13E%24.13E\n", p[0], p[1], p[2]);
	}
	fclose(f1), fclose(f2), fclose(f3);
}

Comments

[退会済みユーザー]

ありがとうございます。
とてもすっきりしました！
ぜひ機会があればまたお願いします。

[退会済みユーザー]

よければ回答してください。
二つ目のコードのpは一つ目と同じくdouble 型で宣言するだけで良いのでしょうか。

Answer 2

質問は、

/* -- ここからわかりません -- */

ここの二重ループのことでしょうか?
という以前にその前のコードに不明点があります。
まず、point1[0][0], point1[0][0]が何かが分からない。きっと、/* ファイル読み込みなどは割愛 */にあるかも思いますが、質問文だけでは不明。

for(int i = 1; i < point1[0][0]+1; i++){//配列に保存

iが 1から始まってますが、Cの場合、普通は、0からです。その一方、for(int i = 0; i < point2[0][0]; i++){ とこちらと範囲が違ってます。もう一つの epoint2[][]に入力している方も同様。
と思ったら、for(int i = 0; i < point2[0][0]; i++){のループは、iが、point2[0][0]にも関わらず、epoint1[i]で参照、組合せが違ってませんか?

また、データの読み込みですが、

  fscanf(readf_in1,"%lf",  //lfは倍精度実数で読み取る

で、一度に、1つのデータしか読み込まない設定なのに、引数が 3つとなっています。後半の 2つには値が入らない、、、と思ったら、 for()ループの都合で、結果 OKか? 配列のサイズも x, y, z の 3つ分で良いのに、epoint1[1024][16] で 16個分は余分ですね。

と、全体的に、見直すべきところは多いように思えます。

Comments

[退会済みユーザー]

回答ありがとうございます。質問も含めて見直します。
もう少し、コードのほうを載せたいと思います。