/Users/jack/Code/basso_dev/inc/basic_functions.h

Go to the documentation of this file.

#ifndef _BASIC_FUNCTIONS_BASSO_H_
#define _BASIC_FUNCTIONS_BASSO_H_

#include "basso.h"

namespace basso {
        
        
        
        Numeric norm( const nArray &a, int type=2 )
        {
                switch ( type )
                {
                        case 0:
                        return vect_norminf(a);
                        
                        case 1:
                        return vect_norm1(a);
                        
                        case 2:
                        return vect_norm2(a);
                        
                        default:
                        WarningMessage("norm","unknown norm type, returning L2 norm");
                        return vect_norm2(a);
                }
        }
        
        Numeric normalize( nArray &a, int type=2 )
        {
                Numeric n=norm(a,type);
                for ( int i=0; i<a.size(); ++i )
                        a[i] = a[i] / n;
                return n;       
        }
        
        void cross_prod( const nArray &a, const nArray &b, nArray &c )
        {
#ifdef ALLOW_DYNAMIC_RESIZE
                if ( vect_size(c)!=3 )
                        resize( c, 3 );
#elif EXPLICIT_BOUNDS_CHECK
                if ( vect_size(c)<3 || vect_size(a)<3 || vect_size(b)<3 )
                        WarningMessage("cross_product","invalid size of vectors, must be of dimension 3");
#endif
                c[0]=a[1]*b[2]-a[2]*b[1];
                c[1]=a[2]*b[0]-a[0]*b[2];
                c[2]=a[0]*b[1]-a[1]*b[0];
        }
        
        void outer_prod ( const nArray &a, const nArray &b, nMatrix &c, Numeric alpha=1.0 )
        {
#ifdef ALLOW_DYNAMIC_RESIZE
                if ( mat_nrows(c)!=vect_size(a) || mat_ncols(c)!=vect_size(b) )
                        resize( c, a.size(), b.size() );
#elif EXPLICIT_BOUNDS_CHECK
                if ( mat_nrows(c)!=vect_size(a) || mat_ncols(c)!=vect_size(b) )
                        WarningMessage("outer_product","incompatible dimensions");
#endif
                for ( int i=0; i<vect_size(a); ++i )
                        for ( int j=0; j<vect_size(b); ++j )
                                c(i,j)=alpha*a[i]*b[j];
        }
        
        
} // end namespace

#endif

---