Clazz.declarePackage("J.jvxl.readers");
Clazz.load(["J.jvxl.readers.SurfaceReader"], "J.jvxl.readers.VolumeDataReader", ["JU.AU", "$.SB", "J.jvxl.data.JvxlCoder", "JU.Logger"], function(){
var c$ = Clazz.decorateAsClass(function(){
this.dataType = 0;
this.precalculateVoxelData = false;
this.allowMapData = false;
this.point = null;
this.ptsPerAngstrom = 0;
this.maxGrid = 0;
this.useOriginStepsPoints = false;
Clazz.instantialize(this, arguments);}, J.jvxl.readers, "VolumeDataReader", J.jvxl.readers.SurfaceReader);
Clazz.makeConstructor(c$, 
function(){
Clazz.superConstructor (this, J.jvxl.readers.VolumeDataReader, []);
});
Clazz.overrideMethod(c$, "init", 
function(sg){
this.initVDR(sg);
}, "J.jvxl.readers.SurfaceGenerator");
Clazz.defineMethod(c$, "initVDR", 
function(sg){
this.initSR(sg);
this.useOriginStepsPoints = (this.params.origin != null && this.params.points != null && this.params.steps != null);
this.dataType = this.params.dataType;
this.precalculateVoxelData = true;
this.allowMapData = true;
}, "J.jvxl.readers.SurfaceGenerator");
Clazz.defineMethod(c$, "setup", 
function(isMapData){
this.jvxlFileHeaderBuffer =  new JU.SB().append("volume data read from file\n\n");
J.jvxl.data.JvxlCoder.jvxlCreateHeaderWithoutTitleOrAtoms(this.volumeData, this.jvxlFileHeaderBuffer);
}, "~B");
Clazz.overrideMethod(c$, "readVolumeParameters", 
function(isMapData){
this.setup(isMapData);
this.initializeVolumetricData();
return true;
}, "~B");
Clazz.overrideMethod(c$, "readVolumeData", 
function(isMapData){
try {
this.readSurfaceData(isMapData);
} catch (e) {
if (Clazz.exceptionOf(e, Exception)){
System.out.println(e.toString());
{
}return false;
} else {
throw e;
}
}
return true;
}, "~B");
Clazz.defineMethod(c$, "readVoxelDataIndividually", 
function(isMapData){
if (isMapData && !this.allowMapData) return;
if (!isMapData || this.volumeData.sr != null) {
this.volumeData.setVoxelDataAsArray(this.voxelData = null);
return;
}this.newVoxelDataCube();
for (var x = 0; x < this.nPointsX; ++x) {
var plane = JU.AU.newFloat2(this.nPointsY);
this.voxelData[x] = plane;
var ptyz = 0;
for (var y = 0; y < this.nPointsY; ++y) {
var strip = plane[y] =  Clazz.newFloatArray (this.nPointsZ, 0);
for (var z = 0; z < this.nPointsZ; ++z, ++ptyz) {
strip[z] = this.getValue(x, y, z, ptyz);
}
}
}
}, "~B");
Clazz.defineMethod(c$, "setVolumeData", 
function(){
});
Clazz.defineMethod(c$, "setVolumeDataParams", 
function(){
if (this.params.volumeData != null) {
this.setVolumeDataV(this.params.volumeData);
return true;
}if (!this.useOriginStepsPoints) {
return false;
}this.volumetricOrigin.setT(this.params.origin);
this.volumetricVectors[0].set(this.params.steps.x, 0, 0);
this.volumetricVectors[1].set(0, this.params.steps.y, 0);
this.volumetricVectors[2].set(0, 0, this.params.steps.z);
this.voxelCounts[0] = Clazz.floatToInt(this.params.points.x);
this.voxelCounts[1] = Clazz.floatToInt(this.params.points.y);
this.voxelCounts[2] = Clazz.floatToInt(this.params.points.z);
if (this.voxelCounts[0] < 1 || this.voxelCounts[1] < 1 || this.voxelCounts[2] < 1) return false;
this.showGridInfo();
return true;
});
Clazz.defineMethod(c$, "showGridInfo", 
function(){
JU.Logger.info("grid origin  = " + this.params.origin);
JU.Logger.info("grid steps   = " + this.params.steps);
JU.Logger.info("grid points  = " + this.params.points);
this.ptTemp.x = this.params.steps.x * this.params.points.x;
this.ptTemp.y = this.params.steps.y * this.params.points.y;
this.ptTemp.z = this.params.steps.z * this.params.points.z;
JU.Logger.info("grid lengths = " + this.ptTemp);
this.ptTemp.add(this.params.origin);
JU.Logger.info("grid max xyz = " + this.ptTemp);
});
Clazz.defineMethod(c$, "setVoxelRange", 
function(index, min, max, ptsPerAngstrom, gridMax, minPointsPerAngstrom){
var nGrid;
var d;
if (min - max >= -1.0E-4) {
min = -10;
max = 10;
}var range = max - min;
var resolution = this.params.resolution;
if (resolution != 3.4028235E38) {
ptsPerAngstrom = resolution;
minPointsPerAngstrom = 0;
}nGrid = Clazz.doubleToInt(Math.floor(range * ptsPerAngstrom)) + 1;
if (nGrid > gridMax) {
if ((this.dataType & 256) > 0) {
if (resolution == 3.4028235E38) {
if (!this.isQuiet) JU.Logger.info("Maximum number of voxels for index=" + index + " exceeded (" + nGrid + ") -- set to " + gridMax);
nGrid = gridMax;
} else {
if (!this.isQuiet) JU.Logger.info("Warning -- high number of grid points: " + nGrid);
}} else if (resolution == 3.4028235E38) {
nGrid = gridMax;
}}if (nGrid == 1) nGrid = 2;
ptsPerAngstrom = (nGrid - 1) / range;
if (ptsPerAngstrom < minPointsPerAngstrom) {
ptsPerAngstrom = minPointsPerAngstrom;
nGrid = Clazz.doubleToInt(Math.floor(ptsPerAngstrom * range + 1));
ptsPerAngstrom = (nGrid - 1) / range;
}d = this.volumeData.volumetricVectorLengths[index] = 1 / ptsPerAngstrom;
this.voxelCounts[index] = nGrid;
if (this.params.sbOut != null) this.params.sbOut.append("isosurface resolution for axis " + (index + 1) + " set to " + ptsPerAngstrom + " points/Angstrom; " + this.voxelCounts[index] + " voxels\n");
switch (index) {
case 0:
this.volumetricVectors[0].set(d, 0, 0);
this.volumetricOrigin.x = min;
break;
case 1:
this.volumetricVectors[1].set(0, d, 0);
this.volumetricOrigin.y = min;
break;
case 2:
this.volumetricVectors[2].set(0, 0, d);
this.volumetricOrigin.z = min;
if (this.isEccentric) this.eccentricityMatrix.rotate(this.volumetricOrigin);
if (this.center != null && !Float.isNaN(this.center.x)) this.volumetricOrigin.add(this.center);
if (this.params.sbOut != null) this.params.sbOut.append((this.voxelCounts[0] * this.voxelCounts[1] * this.voxelCounts[2]) + " voxels total\n");
}
if (this.isEccentric) this.eccentricityMatrix.rotate(this.volumetricVectors[index]);
return this.voxelCounts[index];
}, "~N,~N,~N,~N,~N,~N");
Clazz.overrideMethod(c$, "readSurfaceData", 
function(isMapData){
this.readSurfaceDataVDR(isMapData);
}, "~B");
Clazz.defineMethod(c$, "readSurfaceDataVDR", 
function(isMapData){
if (this.isProgressive && !isMapData) {
this.nDataPoints = this.volumeData.setVoxelCounts(this.nPointsX, this.nPointsY, this.nPointsZ);
this.voxelData = null;
return;
}if (this.precalculateVoxelData) this.generateCube();
 else this.readVoxelDataIndividually(isMapData);
}, "~B");
Clazz.defineMethod(c$, "generateCube", 
function(){
JU.Logger.info("data type: user volumeData");
JU.Logger.info("voxel grid origin:" + this.volumetricOrigin);
for (var i = 0; i < 3; ++i) JU.Logger.info("voxel grid vector:" + this.volumetricVectors[i]);

JU.Logger.info("Read " + this.nPointsX + " x " + this.nPointsY + " x " + this.nPointsZ + " data points");
});
Clazz.overrideMethod(c$, "closeReader", 
function(){
});
});
;//5.0.1-v2 Tue Jul 23 17:25:20 CDT 2024
