ENH: scipy.spatial.distance: use BLAS and waste less memory

This module was using np.linalg.norm instead of the SciPy version,
which switches to BLAS when appropriate.

The change to sqeuclidean makes it four times faster for length-1e5
random vectors using NumPy 1.6.2 linked against ATLAS 3.8.4:

>>> %timeit ((u - v) ** 2).sum()
1000 loops, best of 3: 2.01 ms per loop
>>> %timeit u_v = u - v; np.dot(u_v, u_v)
1000 loops, best of 3: 552 us per loop
DOC: rename the shapes for randint distribution.

Make them match the numpy random_integer generator's.
DOC: rename the shape parameter for Von Mises distribution.

Now it matches the numpy one.
DOC: change the shape name for vonmises in stats tutorial
MAINT: remove LyX files from stats tutorial docs
reduced number of test failure cases for hyperu
BUG: spatial.distance: always convert to float in sqeuclidean

Integers are subject to silent wraparound in the subtraction and
squaring operations, as pointed out by @jakevdp. Some unnecessary
copying may occur, but any call where both arguments are float
arrays still only allocates a single temporary array.
BUG: sparse: fix crashes/wrong results in broadcasting multiply

The dia_matrices used to broadcast column/row vector elementwise
multiplications should be square.  The matrices are passed directly to
the internal _mul_* methods, which make no data size or shape checks, so
wrong size gives wrong results or segfaults with no warnings.
BUG: sparse/arpack: build with second_NONE

Arpack doesn't need the timing routine for anything.  arpack-ng also
builds using second_NONE
BUG: optimize: stub out lbfgsb timer

Also the L-BFGS-B code doesn't really need to know the timings,
as they are not exposed by scipy.
TST: special: fix inf test to be compatible with old Numpy
Merge pull request #2825 from larsmans/blas-spatial-distance

ENH: spatial: use BLAS in scipy.spatial.distance.sqeuclidean
Guarantee oriented simplices for 2D delaunay

This commit updates the existing Delaunay code in order to ensure
that 2D triangulations maintain a consistent ordering. This
is similar to specifying the 'i' flag in QHull. This brings
the behaviour inline with Matlab, and is the most sensible output
for 2D triangulations.

This is achieved by swapping the first and second indices if
the orientation of the facet is clockwise. The logic is written
to try and minimise the impact on higher order triangulations.
Remove the 0 when setting qh_ORIENTclock

The default value for an enum is 0 so the setter is redundant
Fix tests to ensure counterclockwise orientation
Merge pull request #2949 from pv/arscnd

Stub out timer routines in Fortran code
DOC: pep8
Merge pull request #2960 from sciunto/doc_pep8

DOC: fix pep8 issues in tutorials
Merge pull request #2937 from pv/sparse-mul-bugfix

BUG: sparse: fix crashes/wrong results in broadcasting multiply
Merge pull request #2921 from EvgeniBurovski/shapes

Stats: renaming shape parameters for several distributions