00001 #ifndef SimTK_SIMMATRIX_SMALL_DEFS_THAT_NEED_BIG_H_ 00002 #define SimTK_SIMMATRIX_SMALL_DEFS_THAT_NEED_BIG_H_ 00003 00004 /* -------------------------------------------------------------------------- * 00005 * SimTK Core: SimTK Simmatrix(tm) * 00006 * -------------------------------------------------------------------------- * 00007 * This is part of the SimTK Core biosimulation toolkit originating from * 00008 * Simbios, the NIH National Center for Physics-Based Simulation of * 00009 * Biological Structures at Stanford, funded under the NIH Roadmap for * 00010 * Medical Research, grant U54 GM072970. See https://simtk.org. * 00011 * * 00012 * Portions copyright (c) 2005-7 Stanford University and the Authors. * 00013 * Authors: Michael Sherman * 00014 * Contributors: * 00015 * * 00016 * Permission is hereby granted, free of charge, to any person obtaining a * 00017 * copy of this software and associated documentation files (the "Software"), * 00018 * to deal in the Software without restriction, including without limitation * 00019 * the rights to use, copy, modify, merge, publish, distribute, sublicense, * 00020 * and/or sell copies of the Software, and to permit persons to whom the * 00021 * Software is furnished to do so, subject to the following conditions: * 00022 * * 00023 * The above copyright notice and this permission notice shall be included in * 00024 * all copies or substantial portions of the Software. * 00025 * * 00026 * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR * 00027 * IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, * 00028 * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL * 00029 * THE AUTHORS, CONTRIBUTORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, * 00030 * DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR * 00031 * OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE * 00032 * USE OR OTHER DEALINGS IN THE SOFTWARE. * 00033 * -------------------------------------------------------------------------- */ 00034 00041 namespace SimTK { 00042 00043 // TODO: This is really bad! Should share space instead of recopying. 00044 template <int M, int N, class ELT, int CS, int RS> 00045 typename Mat<M,N,ELT,CS,RS>::TInvert 00046 Mat<M,N,ELT,CS,RS>::invert() const { 00047 Matrix_< EStandard > bigm(M,N); 00048 for (int j=0; j<N; ++j) 00049 for (int i=0; i<M; ++i) 00050 bigm(i,j) = (*this)(i,j); 00051 Matrix_< EInvert > result = bigm.invert(); 00052 assert(result.nrow() == TInvert::NRows && result.ncol() == TInvert::NCols); 00053 TInvert out; 00054 for (int j=0; j<TInvert::NCols; ++j) 00055 for (int i=0; i<TInvert::NRows; ++i) 00056 out(i,j) = result(i,j); 00057 return out; 00058 } 00059 00060 } //namespace SimTK 00061 00062 00063 #endif // SimTK_SIMMATRIX_SMALL_DEFS_THAT_NEED_BIG_H_