source/mathematicarepresentation.cc

   1 /* This file is part of Shapes.
   2  *
   3  * Shapes is free software: you can redistribute it and/or modify
   4  * it under the terms of the GNU General Public License as published by
   5  * the Free Software Foundation, either version 3 of the License, or
   6  * any later version.
   7  *
   8  * Shapes is distributed in the hope that it will be useful,
   9  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  10  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  11  * GNU General Public License for more details.
  12  *
  13  * You should have received a copy of the GNU General Public License
  14  * along with Shapes.  If not, see <http://www.gnu.org/licenses/>.
  15  *
  16  * Copyright 2008 Henrik Tidefelt
  17  */
  18
  19 #include "mathematicarepresentation.h"
  20
  21 #include <sstream>
  22
  23 using namespace Shapes;
  24
  25
  26 std::string
  27 Helpers::mathematicaFormat( size_t n1, size_t n2, size_t tda, const double * data )
  28 {
  29         std::ostringstream os;
  30         os << "{" ;
  31         for( const double * row = data; row != data + n1 * tda; row += tda )
  32                 {
  33                         if( row != data )
  34                                 {
  35                                         os << "," ;
  36                                 }
  37                         os << "{" ;
  38                         for( const double * src = row; src != row + n2; ++src )
  39                                 {
  40                                         if( src != row )
  41                                                 {
  42                                                         os << "," ;
  43                                                 }
  44                                         os << *src ;
  45                                 }
  46                         os << "}" ;
  47                 }
  48         os << "}" ;
  49         return os.str( );
  50 }
  51
  52 std::string
  53 Helpers::mathematicaFormatTransposed( size_t n1, size_t n2, size_t tda, const double * data )
  54 {
  55         std::ostringstream os;
  56         os << "{" ;
  57         for( const double * row = data; row != data + n2; ++row )
  58                 {
  59                         if( row != data )
  60                                 {
  61                                         os << "," ;
  62                                 }
  63                         os << "{" ;
  64                         for( const double * src = row; src != row + n1 * tda; src += tda )
  65                                 {
  66                                         if( src != row )
  67                                                 {
  68                                                         os << "," ;
  69                                                 }
  70                                         os << *src ;
  71                                 }
  72                         os << "}" ;
  73                 }
  74         os << "}" ;
  75         return os.str( );
  76 }
  77
  78 std::string
  79 Helpers::mathematicaFormat( size_t n, size_t stride, const double ** begin, const double ** end )
  80 {
  81         std::ostringstream os;
  82         os << "{" ;
  83         for( const double ** row = begin; row != end; ++row )
  84                 {
  85                         if( row != begin )
  86                                 {
  87                                         os << "," ;
  88                                 }
  89                         os << "{" ;
  90                         for( const double * src = *row; src != *row + n * stride; src += stride )
  91                                 {
  92                                         if( src != *row )
  93                                                 {
  94                                                         os << "," ;
  95                                                 }
  96                                         os << *src ;
  97                                 }
  98                         os << "}" ;
  99                 }
 100         os << "}" ;
 101         return os.str( );
 102 }
 103
 104 std::string
 105 Helpers::mathematicaFormat( size_t n, size_t stride, const double * data )
 106 {
 107         std::ostringstream os;
 108         os << "{" ;
 109         for( const double * src = data; src != data + n * stride; src += stride )
 110                 {
 111                         if( src != data )
 112                                 {
 113                                         os << "," ;
 114                                 }
 115                         os << *src ;
 116                 }
 117         os << "}" ;
 118         return os.str( );
 119 }
 120
 121 std::string
 122 Helpers::mathematicaFormat( const gsl_vector & x )
 123 {
 124         return Helpers::mathematicaFormat( x.size, x.stride, x.data );
 125 }
 126
 127 std::string
 128 Helpers::mathematicaFormat( const gsl_matrix & x, bool transpose )
 129 {
 130         if( transpose )
 131                 {
 132                         return Helpers::mathematicaFormatTransposed( x.size1, x.size2, x.tda, x.data );
 133                 }
 134         else
 135                 {
 136                         return Helpers::mathematicaFormat( x.size1, x.size2, x.tda, x.data );
 137                 }
 138 }