Type checked invariants

[vbarrielle]

Use phantom data to ensure at compile time that algorithms expecting e.g. sorted indices get matrices satisfying that property .