ファイルオフセットの計算
指定されたデータ値のオフセット(ファイル内の位置)を計算するには、指定された変数型nc_typeに適切なデータ値の一つの外部サイズ(バイト単位)をexternal_sizeofとします。
NC_BYTE 1
NC_CHAR 1
NC_SHORT 2
NC_INT 4
NC_FLOAT 4
NC_DOUBLE 8
NF_OPEN (もしくはNF_ENDDEF)呼び出しは前もってvar_arrayと示された変数配列内をスキャンし、recsizeを計算するために"記録”変数のvsizeの和を計算します。
変数の次元サイズの積を右から左にとっていき、記録変数の最も左の(記録)次元は飛ばし、各変数についての結果をproduct配列に格納します。例えば:
Non-record variable:
dimension lengths: [ 5 3 2 7]
product: [210 42 14 7]
Record variable:
dimension lengths: [0 2 9 4]
product: [0 72 36 4]
この時点では、最も左にある積を次の4の倍数に丸めたものが変数サイズ、すなわち、上の文法においてはvsizeになります。例えば、上記の非記録変数では、vsizeフィールド値は212(210を次の4の倍数に丸めた値)です。記録変数に対しては、vsizeの値はちょうど72です。なぜならば、72は既に4の倍数であるからです。
求めるデータ値の座標の配列をcoordとし、求める結果をoffsetとします。この時、offsetは単に、求める変数の最初のデータ値のファイルオフセット(そのbeginフィールド)にcoordとproductベクトルの内積を変数の各データのサイズ(バイト単位)を掛けたものを加えた値となります。最後に、もしその変数が記録変数であれば、記録数‘coord[0]’と記録サイズrecsizeとの積が加算され、最終的なオフセット値が導かれます。
擬似Cコードにおけるoffsetの計算は次のようになります。
for (innerProduct = i = 0; i < var.rank; i++)
innerProduct += product[i] * coord[i]
offset = var.begin;
offset += external_sizeof * innerProduct
if(IS_RECVAR(var))
offset += coord[0] * recsize;
したがって、(外部表現法の)データ値を取得するには、次のようになります。
lseek(fd, offset, SEEK_SET);
read(fd, buf, external_sizeof);
特殊例:記録変数が一つしかない場合には、各記録が4倍との境界に合っていなければならないという制限は外れるので、この場合には記録の詰め込みが行なわれません。
Quadralay Corporation http://www.webworks.com Voice: (512) 719-3399 Fax: (512) 719-3606 sales@webworks.com |