"

 COPYRIGHT (c) 1993 by Claus Gittinger

	      All Rights Reserved

 This software is furnished under a license and may be used

 only in accordance with the terms of that license and with the

 inclusion of the above copyright notice.   This software may not

 be provided or otherwise made available to, or used by, any

 other person.  No title to or ownership of the software is

 hereby transferred.

"

"{ Package: 'stx:libbasic' }"

AbstractNumberVector variableFloatSubclass:#FloatArray

	instanceVariableNames:''

	classVariableNames:''

	poolDictionaries:''

	category:'Collections-Arrayed'

!

!FloatArray class methodsFor:'documentation'!

copyright

"

 COPYRIGHT (c) 1993 by Claus Gittinger

	      All Rights Reserved

 This software is furnished under a license and may be used

 only in accordance with the terms of that license and with the

 inclusion of the above copyright notice.   This software may not

 be provided or otherwise made available to, or used by, any

 other person.  No title to or ownership of the software is

 hereby transferred.

"

!

documentation

"

    FloatArrays store floats (and nothing else).

    They have been added to support heavy duty number crunching and

    data exchange with openGL frameworks and other mass data libraries. 

    See documentation in DoubleArray for more information.

    [memory requirements:]

        OBJ-HEADER + (size * float-size)

    [See also:]

        DoubleArray Array

    [author:]

        Claus Gittinger

"

! !

!FloatArray class methodsFor:'queries'!

elementByteSize

    "for bit-like containers, return the number of bytes stored per element.

     Here, 4 is returned"

    ^ 4

    "Created: / 15-09-2011 / 14:12:39 / cg"

! !

!FloatArray methodsFor:'arithmetic'!

* anObject

    "return the product of the receiver and the argument.

     The argument may either be a scalar or another vector"

    ^ self clone *= anObject

    "

     #(1 2 3 4) asFloatArray * 3

     #(1 2 3 4) asFloatArray * #(1 2 3 4) asFloatArray

    "

!

+ anObject

    "return the sum of the receiver and the argument.

     The argument may either be a scalar or another vector"

    ^ self clone += anObject

    "

     #(1 2 3 4) asFloatArray + 3

     #(1 2 3 4) asFloatArray + #(1 2 3 4) asFloatArray

    "

!

- anObject

    "return the difference of the receiver and the argument.

     The argument may either be a scalar or another vector"

    ^ self clone -= anObject

    "

     #(1 2 3 4) asFloatArray - 3

     #(1 2 3 4) asFloatArray - #(1 2 3 4) asFloatArray

    "

    "Created: / 30-05-2007 / 17:41:46 / cg"

!

/ anObject

    "return the division of the receiver and the argument.

     The argument may either be a scalar or another vector"

    ^ self clone /= anObject

    "

     #(1 2 3 4) asFloatArray / 3

     #(1 2 3 4) asFloatArray / #(1 2 3 4) asFloatArray

    "

    "Created: / 30-05-2007 / 17:46:05 / cg"

!

abs

    "return a vector containing absolute values"

    ^ self clone primAbs

    "

     #( -1 2 -3 4 -5 6 -7 8) abs.

     #( -1 2 -3 4 -5 6 -7 8) asFloatArray abs. 

    "

!

negated

    "return a vector containing negated values"

    ^ self clone primNegated

    "

     #( -1 2 -3 4 -5 6 -7 8) negated.     

     #( -1 2 -3 4 -5 6 -7 8) asFloatArray negated. 

    "

    "Modified: / 30-05-2007 / 17:51:47 / cg"

!

primAbs

    "destructive absolute value of each element"

%{

    if (__ClassInstPtr(__qClass(self))->c_ninstvars == __mkSmallInteger(0)) {

	INT _sz = __floatArraySize(self);

	INT i;

	float *_p = __FloatArrayInstPtr(self)->f_element;

	/* how about inline-mmx-asm for this ... */

	for (i=0; i<_sz; i++) {

	    float f = _p[i];

	    if (f < 0) {

		_p[i] = -f;

	    }

	}

	RETURN (self);

    }

%}.

    1 to: self size do:[:i| self at: i put: (self at: i) abs].

    "

     |f|

     f := FloatArray withAll:#(-1 2 -3 4 -5).

     f abs.

     f

    "

    "Created: / 30-05-2007 / 17:50:17 / cg"

!

primAddArray: floatArray

    "add the vector argument into the receiver (destructive).

     The argument must be another vector"

%{

    if ((__ClassInstPtr(__qClass(self))->c_ninstvars == __mkSmallInteger(0))

     && __isFloats(floatArray)

     && (__ClassInstPtr(__qClass(floatArray))->c_ninstvars == __mkSmallInteger(0))) {

	INT _sz1 = __floatArraySize(self);

	INT _sz2 = __floatArraySize(floatArray);

	INT i;

	float *_p1 = __FloatArrayInstPtr(self)->f_element;

	float *_p2 = __FloatArrayInstPtr(floatArray)->f_element;

	if (_sz2 >= _sz1) {

	    /* how about inline-mmx-asm for this ... */

	    for (i=0; i<_sz1; i++) {

		_p1[i] += _p2[i];

	    }

	}

	RETURN (self);

    }

%}.

    1 to: self size do:[:i| self at: i put: (self at: i) + (floatArray at: i)].

    "

     |f1 f2|

     f1 := FloatArray withAll:#(1 2 3 4 5).

     f2 := FloatArray withAll:#(2 2 2 3 3).

     f1 += f2.

     f1

    "

!

primAddScalar: aScalar

    "add the scalar argument into the receiver (destructive)."

%{

    if (__ClassInstPtr(__qClass(self))->c_ninstvars == __mkSmallInteger(0)) {

	INT _sz1 = __floatArraySize(self);

	INT i;

	float *_p1 = __FloatArrayInstPtr(self)->f_element;

	float v;

	if (__isFloat(aScalar)) {

	    v = (float)(__floatVal(aScalar));

	} else  if (__isShortFloat(aScalar)) {

	    v = __shortFloatVal(aScalar);

	} else if (__isSmallInteger(aScalar)) {

	    v = (float)(__intVal(aScalar));

	} else

	    goto badArg;

	/* how about inline-mmx-asm for this ... */

	for (i=0; i<_sz1; i++) {

	    _p1[i] += v;

	}

	RETURN (self);

    }

    badArg: ;

%}.

    1 to: self size do:[:i| self at: i put: (self at: i) + aScalar].

    "

     |f1 f2|

     f1 := FloatArray withAll:#(1 2 3 4 5).

     f1 += 2.0.

     Transcript showCR:f1.

     f1 += 2.0 asShortFloat.

     Transcript showCR:f1.

     f1 += 2.

     Transcript showCR:f1.

    "

!

primDivArray: floatArray

    "divide the vector argument into the receiver (destructive).

     The argument must be another vector"

%{

    if ((__ClassInstPtr(__qClass(self))->c_ninstvars == __mkSmallInteger(0))

     && __isFloats(floatArray)

     && (__ClassInstPtr(__qClass(floatArray))->c_ninstvars == __mkSmallInteger(0))) {

	INT _sz1 = __floatArraySize(self);

	INT _sz2 = __floatArraySize(floatArray);

	INT i;

	float *_p1 = __FloatArrayInstPtr(self)->f_element;

	float *_p2 = __FloatArrayInstPtr(floatArray)->f_element;

	if (_sz2 >= _sz1) {

	    /* how about inline-mmx-asm for this ... */

	    for (i=0; i<_sz1; i++) {

		_p1[i] /= _p2[i];

	    }

	}

	RETURN (self);

    }

%}.

    1 to: self size do:[:i| self at: i put: (self at: i) / (floatArray at: i)].

    "

     |f1 f2|

     f1 := FloatArray withAll:#(1 2 3 4 5).

     f2 := FloatArray withAll:#(2 2 2 3 3).

     f1 /= f2.

     f1

    "

    "Modified: / 29-05-2007 / 16:01:34 / cg"

!

primDivScalar: aScalar

    "divide the scalar argument into the receiver (destructive)."

%{

    if (__ClassInstPtr(__qClass(self))->c_ninstvars == __mkSmallInteger(0)) {

	INT _sz1 = __floatArraySize(self);

	INT i;

	float *_p1 = __FloatArrayInstPtr(self)->f_element;

	float v;

	if (__isFloat(aScalar)) {

	    v = (float)(__floatVal(aScalar));

	} else  if (__isShortFloat(aScalar)) {

	    v = __shortFloatVal(aScalar);

	} else if (__isSmallInteger(aScalar)) {

	    v = (float)(__intVal(aScalar));

	} else

	    goto badArg;

	/* how about inline-mmx-asm for this ... */

	for (i=0; i<_sz1; i++) {

	    _p1[i] /= v;

	}

	RETURN (self);

    }

    badArg: ;

%}.

    1 to: self size do:[:i| self at: i put: (self at: i) / aScalar].

    "

     |f1 f2|

     f1 := FloatArray withAll:#(1 2 3 4 5).

     f1 /= 2.0.

     Transcript showCR:f1.

     f1 /= 2.0 asShortFloat.

     Transcript showCR:f1.

     f1 /= 2.

     Transcript showCR:f1.

    "

    "Modified: / 29-05-2007 / 16:01:39 / cg"

!

primMulArray: floatArray

    "multiply the vector argument into the receiver (destructive).

     The argument must be another vector"

%{

    if ((__ClassInstPtr(__qClass(self))->c_ninstvars == __mkSmallInteger(0))

     && __isFloats(floatArray)

     && (__ClassInstPtr(__qClass(floatArray))->c_ninstvars == __mkSmallInteger(0))) {

	INT _sz1 = __floatArraySize(self);

	INT _sz2 = __floatArraySize(floatArray);

	INT i;

	float *_p1 = __FloatArrayInstPtr(self)->f_element;

	float *_p2 = __FloatArrayInstPtr(floatArray)->f_element;

	if (_sz2 >= _sz1) {

	    /* how about inline-mmx-asm for this ... */

	    for (i=0; i<_sz1; i++) {

		_p1[i] *= _p2[i];

	    }

	}

	RETURN (self);

    }

%}.

    1 to: self size do:[:i| self at: i put: (self at: i) * (floatArray at: i)].

    "

     |f1 f2|

     f1 := FloatArray withAll:#(1 2 3 4 5).

     f2 := FloatArray withAll:#(2 2 2 3 3).

     f1 *= f2.

     f1

    "

!

primMulScalar: aScalar

    "multiply the scalar argument into the receiver (destructive)."

%{

    if (__ClassInstPtr(__qClass(self))->c_ninstvars == __mkSmallInteger(0)) {

	INT _sz1 = __floatArraySize(self);

	INT i;

	float *_p1 = __FloatArrayInstPtr(self)->f_element;

	float v;

	if (__isFloat(aScalar)) {

	    v = (float)(__floatVal(aScalar));

	} else  if (__isShortFloat(aScalar)) {

	    v = __shortFloatVal(aScalar);

	} else if (__isSmallInteger(aScalar)) {

	    v = (float)(__intVal(aScalar));

	} else

	    goto badArg;

	/* how about inline-mmx-asm for this ... */

	for (i=0; i<_sz1; i++) {

	    _p1[i] *= v;

	}

	RETURN (self);

    }

    badArg: ;

%}.

    1 to: self size do:[:i| self at: i put: (self at: i) * aScalar].

    "

     |f1 f2|

     f1 := FloatArray withAll:#(1 2 3 4 5).

     f1 *= 2.0.

     Transcript showCR:f1.

     f1 *= 2.0 asShortFloat.

     Transcript showCR:f1.

     f1 *= 2.

     Transcript showCR:f1.

    "

!

primNegated

    "destructive negative value of each element"

%{

    if (__ClassInstPtr(__qClass(self))->c_ninstvars == __mkSmallInteger(0)) {

	INT _sz = __floatArraySize(self);

	INT i;

	float *_p = __FloatArrayInstPtr(self)->f_element;

	/* how about inline-mmx-asm for this ... */

	for (i=0; i<_sz; i++) {

	    float f = _p[i];

	    _p[i] = -f;

	}

	RETURN (self);

    }

%}.

    1 to: self size do:[:i| self at: i put: (self at: i) negated].

    "

     |f|

     f := FloatArray withAll:#(-1 2 -3 4 -5).

     f negated.

     f

    "

    "Created: / 30-05-2007 / 17:51:29 / cg"

!

primSubtractArray: floatArray

    "subtract the vector argument from the receiver (destructive).

     The argument must be another vector"

%{

    if ((__ClassInstPtr(__qClass(self))->c_ninstvars == __mkSmallInteger(0))

     && __isFloats(floatArray)

     && (__ClassInstPtr(__qClass(floatArray))->c_ninstvars == __mkSmallInteger(0))) {

	INT _sz1 = __floatArraySize(self);

	INT _sz2 = __floatArraySize(floatArray);

	INT i;

	float *_p1 = __FloatArrayInstPtr(self)->f_element;

	float *_p2 = __FloatArrayInstPtr(floatArray)->f_element;

	if (_sz2 >= _sz1) {

	    /* how about inline-mmx-asm for this ... */

	    for (i=0; i<_sz1; i++) {

		_p1[i] -= _p2[i];

	    }

	}

	RETURN (self);

    }

%}.

    1 to: self size do:[:i| self at: i put: (self at: i) - (floatArray at: i)].

    "

     |f1 f2|

     f1 := FloatArray withAll:#(1 2 3 4 5).

     f2 := FloatArray withAll:#(2 2 2 3 3).

     f1 -= f2.

     f1

    "

    "Created: / 30-05-2007 / 17:42:41 / cg"

!

primSubtractScalar: aScalar

    "subtract the scalar argument from the receiver (destructive)."

%{

    if (__ClassInstPtr(__qClass(self))->c_ninstvars == __mkSmallInteger(0)) {