// Construct grid surface representation of data
void Surface::constructGrid(PrimitiveList& primitiveList, const Axes& axes, const Array<double>& displayAbscissa, List<DisplayDataSet>& displayData, ColourScale colourScale)
{
// Forget all data in current primitives
primitiveList.forgetAll();
// Get extents of displayData to use based on current axes limits
Vec3<int> minIndex, maxIndex;
if (!calculateExtents(axes, displayAbscissa, displayData, minIndex, maxIndex)) return;
int nZ = (maxIndex.z - minIndex.z) + 1;
// Copy and transform abscissa values (still in data space) into axes coordinates
Array<double> x(displayAbscissa, minIndex.x, maxIndex.x);
axes.transformX(x);
int nX = x.nItems();
if (nX < 2) return;
// Decide how large to make VertexChunks in Primitives
// We will consider multiple slices at once in order to take most advantage of the post-transform cache
// forming (abscissa.nItems()/(cacheSize-1)+1) Primitives (we divide by (cacheSize-1) because we will force an overlap
// of one gridpoint between adjacent strips).
const int cacheSize = 12;
int nPrimitives = x.nItems()/(cacheSize-1) + 1;
// Get some values from axes so we can calculate colours properly
bool yLogarithmic = axes.logarithmic(1);
double yStretch = axes.stretch(1);
// Reinitialise primitive list
primitiveList.reinitialise(nPrimitives, true, GL_LINES, true);
// Temporary variables
int n, offset = 0, i, nLimit, nMax;
Vec4<GLfloat> colour;
Vec3<double> nrm(0.0, 1.0, 0.0);
Array<double> y;
Array<DisplayDataSet::DataPointType> yType;
DisplayDataSet** slices = displayData.array();
Primitive* currentPrimitive = primitiveList[0];
int verticesA[cacheSize], verticesB[cacheSize];
double z;
// Loop over abscissa indices
while (offset <= x.nItems())
{
// Set nLimit to ensure we don't go beyond the end of the data arrays
nLimit = std::min(cacheSize, x.nItems()-offset);
// Loop over remaining displayData
for (int slice = minIndex.z; slice <= maxIndex.z; ++slice)
{
// Grab arrays
y.copy(slices[slice]->y(), minIndex.x+offset, minIndex.x+offset+nLimit-1);
yType.copy(slices[slice]->yType(), minIndex.x+offset, minIndex.x+offset+nLimit-1);
axes.transformY(y, yType);
z = axes.transformZ(slices[slice]->z());
// Generate vertices for this row
for (n=0; n<nLimit; ++n)
{
i = offset+n;
if (yType.value(n) != DisplayDataSet::NoPoint)
{
// A value exists here, so define a vertex
colourScale.colour(yLogarithmic ? pow(10.0, y.value(n) / yStretch) : y.value(n) / yStretch, colour);
verticesB[n] = currentPrimitive->defineVertex(x.value(i), y.value(n), z, nrm, colour);
// If the previous vertex on this row also exists, draw a line here
if ((n != 0) && (verticesB[n-1] != -1)) currentPrimitive->defineIndices(verticesB[n-1], verticesB[n]);
}
else verticesB[n] = -1;
}
// Draw lines across slices (if slice != 0)
if (slice != 0)
{
nMax = (maxIndex.z-slice) > 1 ? nLimit-1 : nLimit;
for (n=0; n<nMax; ++n)
{
if ((verticesA[n] != -1) && (verticesB[n] != -1)) currentPrimitive->defineIndices(verticesA[n], verticesB[n]);
}
}
// Store current vertices for next pass
for (n=0; n<cacheSize; ++n) verticesA[n] = verticesB[n];
}
// Move to next primitive and increase offset
currentPrimitive = currentPrimitive->next;
offset += cacheSize-1;
}
}
// Write AxisBlock keywords
bool UChromaSession::writeAxisBlock(LineParser& parser, Axes& axes, int axis)
{
parser.writeLineF(" %s %i\n", UChromaSession::viewPaneKeyword(UChromaSession::AxisBlockKeyword), axis);
parser.writeLineF(" %s %s\n", UChromaSession::axisKeyword(UChromaSession::AutoScaleKeyword), Axes::autoScaleMethod(axes.autoScale(axis)));
parser.writeLineF(" %s %s\n", UChromaSession::axisKeyword(UChromaSession::AutoTicksKeyword), stringBool(axes.autoTicks(axis)));
parser.writeLineF(" %s %f\n", UChromaSession::axisKeyword(UChromaSession::FirstTickKeyword), axes.tickFirst(axis));
parser.writeLineF(" %s %s\n", UChromaSession::axisKeyword(UChromaSession::FractionalPositioningKeyword), stringBool(axes.positionIsFractional(axis)));
parser.writeLineF(" %s %s %s %s\n", UChromaSession::axisKeyword(UChromaSession::GridLinesKeyword), stringBool(axes.gridLinesMajor(axis)), stringBool(axes.gridLinesMinor(axis)), stringBool(axes.gridLinesFull(axis)));
LineStyle style = axes.gridLineMajorStyle(axis);
parser.writeLineF(" %s %f '%s' %f %f %f %f\n", UChromaSession::axisKeyword(UChromaSession::GridLineMajorStyleKeyword), style.width(), LineStipple::stipple[style.stipple()].name, style.colour().redF(), style.colour().greenF(), style.colour().blueF(), style.colour().alphaF());
style = axes.gridLineMinorStyle(axis);
parser.writeLineF(" %s %f '%s' %f %f %f %f\n", UChromaSession::axisKeyword(UChromaSession::GridLineMinorStyleKeyword), style.width(), LineStipple::stipple[style.stipple()].name, style.colour().redF(), style.colour().greenF(), style.colour().blueF(), style.colour().alphaF());
parser.writeLineF(" %s %s\n", UChromaSession::axisKeyword(UChromaSession::InvertKeyword), stringBool(axes.inverted(axis)));
parser.writeLineF(" %s %s\n", UChromaSession::axisKeyword(UChromaSession::LabelAnchorKeyword), TextPrimitive::textAnchor(axes.labelAnchor(axis)));
parser.writeLineF(" %s %f %f %f\n", UChromaSession::axisKeyword(UChromaSession::LabelOrientationKeyword), axes.labelOrientation(axis).x, axes.labelOrientation(axis).y, axes.labelOrientation(axis).z);
parser.writeLineF(" %s %f %f\n", UChromaSession::axisKeyword(UChromaSession::LimitsKeyword), axes.min(axis), axes.max(axis));
parser.writeLineF(" %s %s\n", UChromaSession::axisKeyword(UChromaSession::LogarithmicKeyword), stringBool(axes.logarithmic(axis)));
parser.writeLineF(" %s %i\n", UChromaSession::axisKeyword(UChromaSession::MinorTicksKeyword), axes.minorTicks(axis));
NumberFormat fmt = axes.numberFormat(axis);
parser.writeLineF(" %s '%s' %i %s %s\n", UChromaSession::axisKeyword(UChromaSession::NumberFormatKeyword), NumberFormat::formatType(fmt.type()), fmt.nDecimals(), stringBool(fmt.useUpperCaseExponent()), stringBool(fmt.forcePrecedingPlus()));
parser.writeLineF(" %s %f %f %f\n", UChromaSession::axisKeyword(UChromaSession::PositionFractionalKeyword), axes.positionFractional(axis).x, axes.positionFractional(axis).y, axes.positionFractional(axis).z);
parser.writeLineF(" %s %f %f %f\n", UChromaSession::axisKeyword(UChromaSession::PositionRealKeyword), axes.positionReal(axis).x, axes.positionReal(axis).y, axes.positionReal(axis).z);
parser.writeLineF(" %s %f\n", UChromaSession::axisKeyword(UChromaSession::StretchKeyword), axes.stretch(axis));
parser.writeLineF(" %s %f\n", UChromaSession::axisKeyword(UChromaSession::TickDeltaKeyword), axes.tickDelta(axis));
parser.writeLineF(" %s %f %f %f\n", UChromaSession::axisKeyword(UChromaSession::TickDirectionKeyword), axes.tickDirection(axis).x, axes.tickDirection(axis).y, axes.tickDirection(axis).z);
parser.writeLineF(" %s %s\n", UChromaSession::axisKeyword(UChromaSession::TitleAnchorKeyword), TextPrimitive::textAnchor(axes.titleAnchor(axis)));
parser.writeLine(QString(" ")+UChromaSession::axisKeyword(UChromaSession::TitleKeyword)+" '"+axes.title(axis)+"'\n");
parser.writeLineF(" %s %f %f %f %f\n", UChromaSession::axisKeyword(UChromaSession::TitleOrientationKeyword), axes.titleOrientation(axis).x, axes.titleOrientation(axis).y, axes.titleOrientation(axis).z, axes.titleOrientation(axis).w);
parser.writeLineF(" %s %s\n", UChromaSession::axisKeyword(UChromaSession::VisibleAxisKeyword), stringBool(axes.visible(axis)));
parser.writeLineF(" %s\n", UChromaSession::axisKeyword(UChromaSession::EndAxisKeyword));
return true;
}
