DataMatrixを作るためにECC200のエラーコード作成方法についてC->java

Question

JISX0512を参考にエラーコードの作成を使用としています。
C言語での作成方法は乗っていたのですが（後述）現在使用できるエディタがProcessingなのでそれで使えるように書き換えようとしています。

しかし、想定していた挙動がなされないため質問させていただきました。

JISX0512
https://www.kikakurui.com/x0/X0512-2015-01.html
(別紙Eに載っているコード。これをProcessingで再現したい)

/* "prod(x, y, log, alog, gf)" returns the product "x" times "y" */ 
int prod(int x, int y, int *log, int *alog, int gf) { 
　if (!x || !y) return 0; 
　else return alog[(log[x] + log[y]) % (gf-1)]; 
} 
/* "ReedSolomon(wd, nd, nc, gf.pp)" takes "nd" data codeword values in wd[] */ 
/* and adds on "nc" check codewords, all within GF(gf) where "gf" is a */ 
/* power of 2 and "pp" is the value of its prime modulus polynomial */ 
void ReedSolomon(int *wd, int nd, int nc, int gf, int pp) { 
　int i, j, k, *log,*alog,*c; 
/* allocate, then generate the log & antilog arrays: */ 
　log = malloc(sizeof(int) * gf); 
　alog = malloc(sizeof(int) * gf); 

　log[0] = 1-gf; alog[0] = 1; 
　for (i = 1; i < gf; i++) { 
　　alog[i] = alog[i-1] * 2; 
　　if (alog[i] >= gf) alog[i] ^= pp; 
　　log[alog[i]] = i; 
　} 
/* allocate, then generate the generator polynomial coefficients: */ 
　c = malloc(sizeof(int) * (nc+1)); 
　for (i=1; i<=nc; i++) c[i] = 0; c[0] = 1; 
　　for (i=1; i<=nc; i++) { 
　　　c[i] = c[i-1]; 
　　　for (j=i-1; j>=1; j--) { 
　　　　c[j] = c[j-1] ^ prod(c[j],alog[i],log,alog,gf);
　　　} 
　　　c[0] = prod(c[0],alog[i],log,alog,gf);
　　} 
　/* clear, then generate "nc" checkwords in the array wd[] : */ 
　for (i=nd; i<=(nd+nc); i++) wd[i] = 0; 
　for (i=0; i<nd; i++) { 
　　k = wd[nd] ^ wd[i] ; 
　　for (j=0; j<nc; j++) { 
　　　wd[nd+j] = wd[nd+j+l] ^ prod(k,c[nc-j-1],log, alog,gf); 
　　} 
　} 
　free(c); 
　free(alog); 
　free(log); 
} 

webでRuntimeErrorが出たもの
ideoneというサイトで回しました

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

int prod(int x,int y,int *log,int *alog,int gf){
	if(!x || !y)
	return 0;
	else
	return alog[(log[x] + log[y]) %(gf-1)];
}
 
void ReedSolomon(int *wd, int nd, int nc, int gf, int pp) { 
	int i,j,k,*log,*alog,*c;
	log = malloc(sizeof(int) * gf); 
	alog = malloc(sizeof(int) * gf); 
	log[0] = 1-gf; alog[0] = 1; 
	for (i = 1; i < gf; i++) { 
		alog[i] = alog[i-1] * 2; 
		if (alog[i] >= gf) alog[i] ^= pp; 
		log[alog[i]] = i; 
	} 
	c = malloc(sizeof(int) * (nc+1)); 
	
	for (i=1; i<=nc; i++) c[i] = 0; 
	c[0] = 1; 
	for (i=1; i<=nc; i++) { 
		c[i] = c[i-1]; 
		for (j=i-1; j>=1; j--) { 
			c[j] = c[j-1] ^ prod(c[j],alog[i],log,alog,gf); 
		} 
		c[0] = prod(c[0],alog[i],log,alog,gf); 
	} 
	for (i=nd; i<=(nd+nc); i++) {
		wd[i] = 0;
	} 
	for (i=0; i<nd; i++) { 
		k = wd[nd] ^ wd[i]; 
		for (j=0; j<nc; j++) { 
			wd[nd+j] = wd[nd+j+1] ^ prod(k,c[nc-j-1],log, alog,gf); 
		} 
	}
	free(c); 
	free(alog); 
	free(log); 
}
int main()
{
    printf("Hello, World!\n");
    int wd[3]  = {142,164,186};
	ReedSolomon(wd,3,5,28,301);
    return 0;
}

//追記
gfを28としていますが調べてみたところガロア体なのでGF（2^8=256)ではないかと考えています。

processing用に書き換えたもの

void draw(){
 int[] inst = {142,164,186};
 ReedSolomon(inst,3,5,28,301);
}

IntDict logDic = new IntDict();
IntDict alogDic = new IntDict();

int prod(int x, int y, IntDict log, IntDict alog, int gf) {
  if (x == 0 || y == 0) {
    return 0;
  } else {
    return alogDic.get( str(  ( logDic.get(str(x)) + logDic.get(str(y)) )%(gf-1)));
  }
}

void ReedSolomon(int[] wd, int nd, int nc, int gf, int pp) {
  int k;
  logDic.set(str(0), 1-gf);
  alogDic.set(str(0), 1);

  for (int i = 1; i < gf; i++) {
    alogDic.set(str(i), alogDic.get(str(i-1))*2);
    if (alogDic.get(str(i)) >=gf) 
      alogDic.set(str(i), alogDic.get(str(i))^pp);
      //logDic.set( str(i), (  alogDic.get( str(i) )  ) );
　　  logDic.set( str (  alogDic.get( str(i) )) ,(i) );
  }

  int[] c = new int[nc + 1];
  for (int i = 1; i <= nc; i++) {
    c[i] = 0;
  }
  c[0] = 1;
  println("c");
  for (int i = 1; i <= nc; i++) {
    c[i]=c[i-1];
    for (int j = i-1; j >= 1; j--) {
      c[j] = c[j-1]^prod(c[j], alogDic.get(str(i)), logDic, alogDic, gf);
　　　//このprod()で該当するx,yが存在しないとエラーが出る
    }
    c[0] = prod(c[0], alogDic.get(str(i)), logDic, alogDic, gf);
  }
  
  int[] wd2 = new int[nd + nc+1];
  for (int i = 0; i < nd; i++) {
    wd2[i] = wd[i];
  }
  for (int i=nd; i <= (nd+nc); i++) {
    wd2[i] = 0;
  }
  for (int i=0; i < nd; i++) { 
    k = wd2[nd] ^ wd2[i] ; 
    for (int j=0; j < nc; j++) { 
      wd2[nd+j] = wd2[nd+j+1] ^ prod(k, c[nc-j-1], logDic, alogDic, gf);
    }
  }
}
//wd2の中に142,164,186,114,25,5,88,102が入ってほしい

ーーー
現在起こっていること
Processingでは配列の動的確保の方法が異なるためDictionaryを利用。
配列の追加についても別の配列を新しく作成して代用
書き換えはおそらくできたが情報が足りないということでエラーが発生する。
No key named '823'
変数logに該当する値がないとのこと。

ちなみにそもそもCのほうがあっているかどうか試してみようと思いWebのほうではありますが試したところRuntimeErrorというのが出て検証できませんでした。

足りない情報等あれば指摘していただけると助かります。

Comments

[dodox86]

> ちなみにそもそもCのほうがあっているかどうか試してみよう

もともとのサンプルのC言語のソースファイルですが、間違っている（コンパイルエラー）ようです。下の方の以下のコード部分、
> for (j=0; j<nc; j++) {
> 　　　wd[nd+j] = wd[nd+j+l] ^ prod(k,c[nc-j-1],log, alog,gf);
> }

wd[nd + j + l] とありますが、l（小文字のL)は未定義の変数で、コンパイルエラーになります。数字の1（いち)の間違いなのか、それとも本当に未定義（それにしては他の部分で適切そうな値をセットしている部分が無い）なのか。この部分を直してみないとC言語のコードの実行結果と比べることはできなさそうです。

[HogHog]

こちら１でした。
Webに書くとき私もエラーになりまして１に修正したのですが元のコードを載せるため参考ページより取り込みました。
今一度webで回せる状態のコードを書いていますが少々エラーが出ているのでお待ちください

[HogHog]

コードの方追記いたしました

Answer 1

直接の回答ではないです（算数どまりなので中身はさっぱりわかりません^^;

ProcessingはJavaコードも動きます（Java5相当）ので、Javaで調べたほうがいいと思います。

現在リードソロモン符号をjavaで書こうとしています。 - や... - Yahoo!知恵袋

こちらの回答のようにZXingから関連コードを拝借してテストしてみました。
５つのファイルを追加時に、.javaと拡張子を付けるのを忘れないようにしてください。

テストコードだけ質問に合わせました。

import java.util.Arrays;

final int ECC_COUNT = 5;
final int[] MESSAGE = { 142, 164, 186 };

int[] data = new int[MESSAGE.length + ECC_COUNT];
System.arraycopy(MESSAGE, 0, data, 0, MESSAGE.length);
ReedSolomonEncoder encoder = new ReedSolomonEncoder(GenericGF.DATA_MATRIX_FIELD_256);
encoder.encode(data, ECC_COUNT);
println("encoded:\n" + Arrays.toString(data));

data[0] = 0;
data[1] = 0;
//data[5] = 0;
println("changed:\n" + Arrays.toString(data));

ReedSolomonDecoder decoder = new ReedSolomonDecoder(GenericGF.DATA_MATRIX_FIELD_256);
try {
  decoder.decode(data, ECC_COUNT);
  println("decoded:\n" + Arrays.toString(data));
} 
catch (ReedSolomonException e) {
  println("デコードに失敗しました");
}

encoded:
[142, 164, 186, 114, 25, 5, 88, 102]
changed:
[0, 0, 186, 114, 25, 5, 88, 102]
decoded:
[142, 164, 186, 114, 25, 5, 88, 102]

Answer 2

できなかった原因はガロア値を28にしていたことによるものでした。
下に修正済みコードを入れます。

void setup(){
  int[] instaa = {142, 164, 186};
  ReedSolomon(instaa,3,5,256,301);
}
void ReedSolomon(int[] wd, int nd, int nc, int gf, int pp) {
  int k;
  IntDict logDic = new IntDict();
  IntDict alogDic = new IntDict();
  logDic.set(str(0), 1-gf);
  alogDic.set(str(0), 1);
  for (int i = 1; i < gf; i++) {
    alogDic.set(str(i), alogDic.get(str(i-1)) * 2);
    if (alogDic.get(str(i)) >= gf) {
      alogDic.set(str(i), alogDic.get(str(i)) ^ pp);
    }
    logDic.set( str  (  alogDic.get( str(i) )), (i));
  }

  int[] c = new int[nc + 1];
  for (int i = 1; i <= nc; i++) {
    c[i] = 0;
  }
  c[0] = 1;

  for (int i = 1; i <= nc; i++) {
    c[i]=c[i-1];
    for (int j = i-1; j >= 1; j--) {
      c[j] = c[j-1]^prod(c[j], alogDic.get(str(i)), logDic, alogDic, gf);
    }
    c[0] = prod(c[0], alogDic.get(str(i)), logDic, alogDic, gf);
  }

  int[] wd2 = new int[nd + nc + 1];
  for (int i = 0; i < nd; i++) {
    wd2[i] = wd[i];
  }
  for (int i=nd; i <= (nd + nc); i++) {
    wd2[i] = 0;
  }
  for (int i=0; i < nd; i++) { 
    k = wd2[nd] ^ wd2[i];
    for (int j=0; j < nc; j++) { 
      wd2[nd+j] = wd2[nd+j+1] ^ prod(k, c[nc-j-1], logDic, alogDic, gf);
    }
  }
  println(wd2);
}

wd2の配列が一つ多くなっているのは減らすと一番下のforで突っかかってしまうためです。
まだ改善すべき点はありますが一応解決したので終わりとさせていただきます。
ありがとうございます。