context: // Vector4f.cs
  //
  // Author:
  //   Rodrigo Kumpera (rkumpera@novell.com)
  //
  // (C) 2008 Novell, Inc. (http://www.novell.com)
  //
  // Permission is hereby granted, free of charge, to any person obtaining
  // a copy of this software and associated documentation files (the
  // "Software"), to deal in the Software without restriction, including
  // without limitation the rights to use, copy, modify, merge, publish,
  // distribute, sublicense, and/or sell copies of the Software, and to
  // permit persons to whom the Software is furnished to do so, subject to
  // the following conditions:
  //
  // The above copyright notice and this permission notice shall be
  // included in all copies or substantial portions of the Software.
  //
  // THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND,
  // EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF
  // MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND
  // NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE
  // LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION
  // OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION
  // WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
  //
  
  using System;
  using System.Runtime.InteropServices;
  
  namespace Mono.Simd
  {
  	public enum ShuffleSel
  	{
  		XFromX,
  		XFromY,
  		XFromZ,
  		XFromW,
  
  		YFromX = 0x00,
  		YFromY = 0x04,
  		YFromZ = 0x08,
  		YFromW = 0x0C,
  
  		ZFromX = 0x00,
  		ZFromY = 0x10,
  		ZFromZ = 0x20,
  		ZFromW = 0x30,
  
  		WFromX = 0x00,
  		WFromY = 0x40,
  		WFromZ = 0x80,
  		WFromW = 0xC0,
  
  		/*Expand a single element into all elements*/
  		ExpandX = XFromX | YFromX | ZFromX | WFromX,
  		ExpandY = XFromY | YFromY | ZFromY | WFromY,
  		ExpandZ = XFromZ | YFromZ | ZFromZ | WFromZ,
  		ExpandW = XFromW | YFromW | ZFromW | WFromW,
  
  		/*Expand a pair of elements (x,y,z,w) -> (x,x,y,y)*/
  		ExpandXY = XFromX | YFromX | ZFromY | WFromY,
  		ExpandZW = XFromZ | YFromZ | ZFromW | WFromW,
  
  		/*Expand interleaving elements (x,y,z,w) -> (x,y,x,y)*/
  		ExpandInterleavedXY = XFromX | YFromY | ZFromX | WFromY,
  		ExpandInterleavedZW = XFromZ | YFromW | ZFromZ | WFromW,
  
  		/*Rotate elements*/
  		RotateRight = XFromY | YFromZ | ZFromW | WFromX,
  		RotateLeft = XFromW | YFromX | ZFromY | WFromZ,
  
  		/*Swap order*/
  		Swap = XFromW | YFromZ | ZFromY | WFromX,
  	};
  
  /*
  	TODO:
          Unary - (implemented as mulps [-1,-1,-1,-1])
          Abs (implemented as pand [7fffffff,...] )
          Comparison functions
          Mask extraction function
          vector x float ops
  		Replace Shuffle with less bug prone methods
  */
  
  	[StructLayout(LayoutKind.Explicit, Pack = 0, Size = 16)]
  	public struct Vector4f
  	{
  		[ FieldOffset(0) ]
  		internal float x;
  		[ FieldOffset(4) ]
  		internal float y;
  		[ FieldOffset(8) ]
  		internal float z;
  		[ FieldOffset(12) ]
  		internal float w;
  
  		public float X { get { return x; } set { x = value; } }
  		public float Y { get { return y; } set { y = value; } }
  		public float Z { get { return z; } set { z = value; } }
  		public float W { get { return w; } set { w = value; } }
  
  		public static Vector4f Pi 
  		{
  			get { return new Vector4f ((float)System.Math.PI); }
  		}
  
  		public static Vector4f E 
  		{
  			get { return new Vector4f ((float)System.Math.E); }
  		}
  
  		public  static Vector4f One 
  		{
  			get { return new Vector4f (1); }
  		}
  
  		public  static Vector4f Zero 
  		{
  			get { return new Vector4f (0); }
  		}
  
  		public  static Vector4f MinusOne 
  		{
  			get { return new Vector4f (-1); }
  		}
  
  		[System.Runtime.CompilerServices.IndexerName ("Component")]
  		public unsafe float this [int index]
  		{
  			get {
  				if ((index | 0x3) != 0x3) //index < 0 || index > 3
  					throw new ArgumentOutOfRangeException ("index");
  				fixed (float *v = &x) {
  					return * (v + index);
  				}
  			}
  			set {
  				if ( (index | 0x3) != 0x3) //index < 0 || index > 3
  					throw new ArgumentOutOfRangeException ("index");
  				fixed (float *v = &x) {
  					* (v + index) = value;
  				}
  			}
  		}
  
  		public Vector4f (float x, float y, float z, float w)
  		{
  			this.x = x;
  			this.y = y;
  			this.z = z;
  			this.w = w;
  		}
  		
  		public Vector4f (float f)
  		{
  			this.x = f;
  			this.y = f;
  			this.z = f;
  			this.w = f;
  		}
  
  		[Acceleration (AccelMode.SSE1)]
  		public static unsafe Vector4f operator & (Vector4f v1, Vector4f v2)
  		{
  			Vector4f res = new Vector4f ();
  			int *a = (int*)&v1;
  			int *b = (int*)&v2;
  			int *c = (int*)&res;
  			*c++ = *a++ & *b++;
  			*c++ = *a++ & *b++;
  			*c++ = *a++ & *b++;
  			*c = *a & *b;
  			return res;
  		}
  
  		[Acceleration (AccelMode.SSE1)]
  		public static unsafe Vector4f operator | (Vector4f v1, Vector4f v2)
  		{
  			Vector4f res = new Vector4f ();
  			int *a = (int*)&v1;
  			int *b = (int*)&v2;
  			int *c = (int*)&res;
  			*c++ = *a++ | *b++;
  			*c++ = *a++ | *b++;
  			*c++ = *a++ | *b++;
  			*c = *a | *b;
  			return res;
  		}
  
  		[Acceleration (AccelMode.SSE1)]
  		public static unsafe Vector4f operator ^ (Vector4f v1, Vector4f v2)
  		{
  			Vector4f res = new Vector4f ();
  			int *a = (int*)&v1;
  			int *b = (int*)&v2;
  			int *c = (int*)&res;
  			*c++ = *a++ ^ *b++;
  			*c++ = *a++ ^ *b++;
  			*c++ = *a++ ^ *b++;
  			*c = *a ^ *b;
  			return res;
  		}
  
  		[Acceleration (AccelMode.SSE1)]
  		public static Vector4f operator + (Vector4f v1, Vector4f v2)
  		{
  			return new Vector4f (v1.x + v2.x, v1.y + v2.y, v1.z + v2.z, v1.w + v2.w);
  		}
  
  		[Acceleration (AccelMode.SSE1)]
  		public static Vector4f operator - (Vector4f v1, Vector4f v2)
  		{
  			return new Vector4f (v1.x - v2.x, v1.y - v2.y, v1.z - v2.z, v1.w - v2.w);
  		}
  
  		[Acceleration (AccelMode.SSE1)]
  		public static Vector4f operator * (Vector4f v1, Vector4f v2)
  		{
  			return new Vector4f (v1.x * v2.x, v1.y * v2.y, v1.z * v2.z, v1.w * v2.w);
  		}
  
  		[Acceleration (AccelMode.SSE1)]
  		public static Vector4f operator * (float scalar, Vector4f v)
  		{
  			return new Vector4f (scalar * v.x, scalar * v.y, scalar * v.z, scalar * v.w);
  		}
  
  		[Acceleration (AccelMode.SSE1)]
  		public static Vector4f operator * (Vector4f v, float scalar)
  		{
  			return new Vector4f (scalar * v.x, scalar * v.y, scalar * v.z, scalar * v.w);
  		}
  
  		[Acceleration (AccelMode.SSE1)]
  		public static Vector4f operator / (Vector4f v1, Vector4f v2)
  		{
  			return new Vector4f (v1.x / v2.x, v1.y / v2.y, v1.z / v2.z, v1.w / v2.w);
  		}
  
  		[Acceleration (AccelMode.SSE2)]
  		public static bool operator ==(Vector4f v1, Vector4f v2)
  		{
  			return v1.x == v2.x && v1.y == v2.y && v1.z == v2.z && v1.w == v2.w;
  		}
  
  		[Acceleration (AccelMode.SSE2)]
  		public static bool operator !=(Vector4f v1, Vector4f v2)
  		{
  			return v1.x != v2.x || v1.y != v2.y || v1.z != v2.z || v1.w != v2.w;
  		}
  
  		[Acceleration (AccelMode.SSE1)]
  		public static unsafe explicit operator Vector2d (Vector4f v)
  		{
  			Vector2d* p = (Vector2d*)&v;
  			return *p;
  		}
  
  		[Acceleration (AccelMode.SSE1)]
  		public static unsafe explicit operator Vector2l (Vector4f v)
  		{
  			Vector2l* p = (Vector2l*)&v;
  			return *p;
  		}
  
  		[Acceleration (AccelMode.SSE1)]
  		[CLSCompliant(false)]
  		public static unsafe explicit operator Vector2ul (Vector4f v)
  		{
  			Vector2ul* p = (Vector2ul*)&v;
  			return *p;
  		}
  
  		[Acceleration (AccelMode.SSE1)]
  		public static unsafe explicit operator Vector4i (Vector4f v)
  		{
  			Vector4i* p = (Vector4i*)&v;
  			return *p;
  		}
  
  		[Acceleration (AccelMode.SSE1)]
  		[CLSCompliant(false)]
  		public static unsafe explicit operator Vector4ui (Vector4f v)
  		{
  			Vector4ui* p = (Vector4ui*)&v;
  			return *p;
  		}
  
  		[Acceleration (AccelMode.SSE1)]
  		public static unsafe explicit operator Vector8s (Vector4f v)
  		{
  			Vector8s* p = (Vector8s*)&v;
  			return *p;
  		}
  
  		[Acceleration (AccelMode.SSE1)]
  		[CLSCompliant(false)]
  		public static unsafe explicit operator Vector8us (Vector4f v)
  		{
  			Vector8us* p = (Vector8us*)&v;
  			return *p;
  		}
  
  		[Acceleration (AccelMode.SSE1)]
  		[CLSCompliant(false)]
  		public static unsafe explicit operator Vector16sb (Vector4f v)
  		{
  			Vector16sb* p = (Vector16sb*)&v;
  			return *p;
  		}
  
  		[Acceleration (AccelMode.SSE1)]
  		public static unsafe explicit operator Vector16b (Vector4f v)
  		{
  
gt: 			Vector16b* p = (Vector16b*)&v;