GroovX: src/mtx/mtxops.cc Source File

00001 
00003 
00004 //
00005 // Copyright (c) 2002-2004 California Institute of Technology
00006 // Copyright (c) 2004-2007 University of Southern California
00007 // Rob Peters <rjpeters at usc dot edu>
00008 //
00009 // created: Mon Mar  4 11:13:05 2002
00010 // commit: $Id: mtxops.cc 10065 2007-04-12 05:54:56Z rjpeters $
00011 // $HeadURL: file:///lab/rjpeters/svnrepo/code/trunk/groovx/src/mtx/mtxops.cc $
00012 //
00013 // --------------------------------------------------------------------
00014 //
00015 // This file is part of GroovX.
00016 //   [http://ilab.usc.edu/rjpeters/groovx/]
00017 //
00018 // GroovX is free software; you can redistribute it and/or modify it
00019 // under the terms of the GNU General Public License as published by
00020 // the Free Software Foundation; either version 2 of the License, or
00021 // (at your option) any later version.
00022 //
00023 // GroovX is distributed in the hope that it will be useful, but
00024 // WITHOUT ANY WARRANTY; without even the implied warranty of
00025 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
00026 // General Public License for more details.
00027 //
00028 // You should have received a copy of the GNU General Public License
00029 // along with GroovX; if not, write to the Free Software Foundation,
00030 // Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA.
00031 //
00033 
00034 #ifndef GROOVX_PKGS_MTX_MTXOPS_CC_UTC20050626084022_DEFINED
00035 #define GROOVX_PKGS_MTX_MTXOPS_CC_UTC20050626084022_DEFINED
00036 
00037 #include "mtxops.h"
00038 
00039 #include "mtx/arithfunctor.h"
00040 #include "mtx/mtx.h"
00041 
00042 #include "rutz/error.h"
00043 #include "rutz/rand.h"
00044 
00045 #include "rutz/trace.h"
00046 
00047 mtx rand_mtx(int mrows, int ncols)
00048 {
00049 GVX_TRACE("rand_mtx");
00050 
00051   static rutz::urand generator(rutz::default_rand_seed);
00052 
00053   mtx result = mtx::uninitialized(mrows, ncols);
00054 
00055   for (mtx::colmaj_iter
00056          itr = result.colmaj_begin_nc(),
00057          stop = result.colmaj_end_nc();
00058        itr != stop;
00059        ++itr)
00060     {
00061       *itr = generator.fdraw();
00062     }
00063 
00064   return result;
00065 }
00066 
00067 mtx squared(const mtx& src)
00068 {
00069 GVX_TRACE("squared");
00070 
00071   mtx result(src);
00072   // FIXME would be better if mtx_base had a transform() template
00073   result.apply(dash::square());
00074 
00075   return result;
00076 }
00077 
00078 mtx zeropad(const mtx& src, int new_mrows, int new_ncols,
00079             int* ppadtop, int* ppadleft)
00080 {
00081   GVX_TRACE("zeropad");
00082 
00083   if (new_mrows < src.mrows() || new_ncols < src.ncols())
00084     throw rutz::error("zeropad(): new size must be >= old size",
00085                       SRC_POS);
00086 
00087   mtx result = mtx::zeros(new_mrows, new_ncols);
00088 
00089   const int padtop =  (new_mrows - src.mrows() + 1) / 2;
00090   const int padleft = (new_ncols - src.ncols() + 1) / 2;
00091 
00092   if (ppadtop) *ppadtop = padtop;
00093   if (ppadleft) *ppadleft = padleft;
00094 
00095   for (int col = 0; col < src.ncols(); ++col)
00096     {
00097       result.column(col + padleft)(range_n(padtop, src.mrows()))
00098         = src.column(col);
00099     }
00100 
00101   return result;
00102 }
00103 
00104 static const char __attribute__((used)) vcid_groovx_pkgs_mtx_mtxops_cc_utc20050626084022[] = "$Id: mtxops.cc 10065 2007-04-12 05:54:56Z rjpeters $ $HeadURL: file:
00105 #endif // !GROOVX_PKGS_MTX_MTXOPS_CC_UTC20050626084022_DEFINED