Difference between revisions of "DIP58"

From D Wiki
Jump to: navigation, search
Line 116: Line 116:
 
plot(inputSpace, someFunction(inputSpace));
 
plot(inputSpace, someFunction(inputSpace));
 
</syntaxhighlight>
 
</syntaxhighlight>
 +
"meshgrid", as implemented in [http://www.mathworks.com/help/matlab/ref/meshgrid.html MATLAB] and [http://docs.scipy.org/doc/numpy/reference/generated/numpy.meshgrid.html NumPy], constructs a multidimensional range from a list of 1-dimensional ranges.
  
 
=== Simulating Physics ===
 
=== Simulating Physics ===

Revision as of 06:37, 17 March 2014

Title: ".." as a Binary Operator
DIP: 58
Version: 1
Status: Draft
Created: 2014-03-16
Last Modified: 2014-03-16
Author: Mason McGill

Abstract

Making ".." a binary operator allows multidimensional indexing, terse range construction, and other features especially useful for scientific and engineering applications, while reducing language complexity and preserving backwards compatibility.

Rationale

One of D's design goals is to "cater to the needs of numerical analysis programmers". To that end, D already has some syntax to support MATLAB/Python/Julia/R-style multidimensional array manipulation: IndexExpression, SliceExpression, and "$". However, the ability to tersely express numeric ranges arguably makes languages like Julia more usable for "vectorized" numerical programming than D. Making ".." a general-purpose binary operator would enable this functionality, and improve the usability of numerical APIs.

Proposal

Grammar

SliceExpression and ForeachRangeStatement are removed from the language. CaseRangeStatement remains unchanged.

Code in the form "Expression1 .. Expression2" is parsed as a ToExpression (a..b can be read "a to b"). ".." has the lowest operator precedence of any binary operator.

Resolving ToExpression

A ToExpression is delegated to the arguments if they implement "opBinary" or "opBinaryLeft" for the ".." operator. Otherwise, it yields a structure that can be used in an IndexExpression or ForeachStatement to mimic a SliceExpression or ForeachRangeExpression.

auto __evaluateToExpression(Left, Right)(Left left, Right right)
{
    static if (__traits(compiles, left.opBinary!".."(right)))
        return left.opBinary!".."(right);
    else static if (__traits(compiles, right.opBinaryRight!".."(left)))
        return right.opBinaryRight!".."(left);
    else
        return BoundedRange!(Left, Right)(left, right);
}

struct BoundedRange(Front, Back)
{
    Front front;
    Back back;

    static if (is(typeof(front >= back) : bool))
        bool empty() { return front >= back; }

    static if (__traits(compiles, front++))
        void popFront() { front++; }

    /* Other range operations, if supported. */
}

This behavior can be changed in a non-backwards-compatible release to make ".." behave analogously to other binary operators (refusing to compile if a matching "opBinary" or "opBinaryRight" isn't defined).

Changes to the resolution of IndexExpression

An IndexExpression that would previously have been parsed as a SliceExpression first attempts to delegate to "opSlice". Failing that, it delegates to "opIndex". Index-assignment, index-op-assignment, and index-unary operations are handled analogously. This is similar to the route Python took when changing its indexing/slicing semantics.

auto __evaluateIndexExpression(Base, Indices...)(Base base, Indices indices)
{
    enum is0ArgSlice = !indices.length;
    enum is2ArgSlice = is(indices[0] == BoundedRange!(F, B), F, B);

    static if (is2ArgSlice)
        auto front = indices[0].front, back = indices[0].back;

    static if (is0ArgSlice && __traits(compiles, base.opSlice())
        return base.opSlice();
    else static if (is2ArgSlice && __traits(compiles, base.opSlice(front, back))
        return base.opSlice(front, back);
    else
        return base.opIndex(indices);
}

These changes can be reverted in a non-backwards-compatible release to drop support for the "opSlice*" family of functions.

Usage

Multidimensonal Slicing

const submatrix = matrix[0..5, 1..3];
matrix[0..5, 1..3] *= 2;

Strided Slicing/Iteration

const oddEntries = vector[(0..10).by(2)]; // Either spelling could be
const evenEntries = vector[1..2..11];     // supported by a library.

writeln(array(0..2..10)); // Prints [0, 2, 4, 6, 8].

Terse Constant Declaration

enum size = [150.cm .. 200.cm, 25.cm .. 50.cm];
enum orientation = -30.deg .. 30.deg;

detectPedestrians(size, orientation);

Defining Multidimensional Grids

const inputSpace = meshgrid(0..100, 0..100);
plot(inputSpace, someFunction(inputSpace));

"meshgrid", as implemented in MATLAB and NumPy, constructs a multidimensional range from a list of 1-dimensional ranges.

Simulating Physics

const line = Point(1, 2, 3)..Point(3, 2, 1);
const bounds = Box(1..3, 2..4, 3..6);
const collision = bounds.contains(line);

Parsing Text

const letters = 'a'..'z';
const numbers = '0'..'9';

find(letters ~ numbers, userInput);

Handling Dates/Times

const vacation = july(10, 2014)..july(15, 2014);

foreach (day; monday..friday)
    work(day);

Copyright

This document has been placed in the Public Domain.