7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh#ifdef HAVE_CONFIG_H

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh# include <config.h>

b5b35fce2e3df933e5223ef6645d814eacf51cfamiklosh#endif

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh#include "odf.h"

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh#include <stdio.h>

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh#include <time.h>

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh#include <vector>

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh#include "clear-n_.h"

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh#include "inkscape.h"

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh#include <style.h>

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh#include "display/curve.h"

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh#include "libnr/n-art-bpath.h"

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh#include "extension/system.h"

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh#include "xml/repr.h"

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh#include "xml/attribute-record.h"

f9325af537ca5517eb50ef95f432a3204616f6b3apenner#include "sp-image.h"

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh#include "sp-gradient.h"

c53f16f52840e8c0f2be9c1cc3af633c0ba1552emiklosh#include "sp-stop.h"

c53f16f52840e8c0f2be9c1cc3af633c0ba1552emiklosh#include "gradient-chemistry.h"

c53f16f52840e8c0f2be9c1cc3af633c0ba1552emiklosh#include "sp-linear-gradient.h"

c53f16f52840e8c0f2be9c1cc3af633c0ba1552emiklosh#include "sp-radial-gradient.h"

c53f16f52840e8c0f2be9c1cc3af633c0ba1552emiklosh#include "sp-path.h"

c53f16f52840e8c0f2be9c1cc3af633c0ba1552emiklosh#include "sp-text.h"

c53f16f52840e8c0f2be9c1cc3af633c0ba1552emiklosh#include "sp-flowtext.h"

c53f16f52840e8c0f2be9c1cc3af633c0ba1552emiklosh#include "svg/svg.h"

c53f16f52840e8c0f2be9c1cc3af633c0ba1552emiklosh#include "text-editing.h"

c53f16f52840e8c0f2be9c1cc3af633c0ba1552emiklosh

c53f16f52840e8c0f2be9c1cc3af633c0ba1552emiklosh

c53f16f52840e8c0f2be9c1cc3af633c0ba1552emiklosh#include "dom/dom.h"

c53f16f52840e8c0f2be9c1cc3af633c0ba1552emiklosh#include "dom/util/ziptool.h"

c53f16f52840e8c0f2be9c1cc3af633c0ba1552emiklosh#include "dom/io/domstream.h"

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh#include "dom/io/bufferstream.h"

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh#include "dom/io/stringstream.h"

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

9dc68827cbd515262ecb8d5ae8547d9e82c72e00Jon A. Cruz

9dc68827cbd515262ecb8d5ae8547d9e82c72e00Jon A. Cruz

1cda9431ef400135f5e1bd899a94b921bdad0eafmiklosh

3711b3e25395437ee0a09dbbb2a76d999c4ef322mikloshnamespace Inkscape

68664e00e2372534b4df2fdc5f54f836bafece18miklosh{

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7mikloshnamespace Extension

9dc68827cbd515262ecb8d5ae8547d9e82c72e00Jon A. Cruz{

3711b3e25395437ee0a09dbbb2a76d999c4ef322mikloshnamespace Internal

a4d12a5147f3d1d6b568a326e39ef5dca384248dmiklosh{

1667116521643e2475184b048e0abb77a2aa9735miklosh

cb814cb0df20053ca3ef16ce55da474435daf045miklosh

cb814cb0df20053ca3ef16ce55da474435daf045miklosh

cb814cb0df20053ca3ef16ce55da474435daf045mikloshtypedef org::w3c::dom::DOMString DOMString;

75b857d473541532819bd791105cb352c9a43214buliabyaktypedef org::w3c::dom::XMLCh XMLCh;

7a7fa095a483e8b652af9f00e5169f62c84f09b9mikloshtypedef org::w3c::dom::io::OutputStreamWriter OutputStreamWriter;

7a7fa095a483e8b652af9f00e5169f62c84f09b9mikloshtypedef org::w3c::dom::io::BufferOutputStream BufferOutputStream;

7a7fa095a483e8b652af9f00e5169f62c84f09b9mikloshtypedef org::w3c::dom::io::StringOutputStream StringOutputStream;

1cda9431ef400135f5e1bd899a94b921bdad0eafmiklosh

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh#include <math.h>

cb814cb0df20053ca3ef16ce55da474435daf045miklosh

3711b3e25395437ee0a09dbbb2a76d999c4ef322mikloshclass SVDMatrix

1667116521643e2475184b048e0abb77a2aa9735miklosh{

7a7fa095a483e8b652af9f00e5169f62c84f09b9mikloshpublic:

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh SVDMatrix()

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh init();

1667116521643e2475184b048e0abb77a2aa9735miklosh }

1667116521643e2475184b048e0abb77a2aa9735miklosh

1667116521643e2475184b048e0abb77a2aa9735miklosh SVDMatrix(unsigned int rowSize, unsigned int colSize)

1667116521643e2475184b048e0abb77a2aa9735miklosh {

1667116521643e2475184b048e0abb77a2aa9735miklosh init();

1667116521643e2475184b048e0abb77a2aa9735miklosh rows = rowSize;

1667116521643e2475184b048e0abb77a2aa9735miklosh cols = colSize;

1667116521643e2475184b048e0abb77a2aa9735miklosh size = rows * cols;

75b857d473541532819bd791105cb352c9a43214buliabyak d = new double[size];

75b857d473541532819bd791105cb352c9a43214buliabyak for (unsigned int i=0 ; i<size ; i++)

75b857d473541532819bd791105cb352c9a43214buliabyak d[i] = 0.0;

75b857d473541532819bd791105cb352c9a43214buliabyak }

75b857d473541532819bd791105cb352c9a43214buliabyak

84d6d1f7365e49f2936df9df890ce49d2c000ce2Kris SVDMatrix(double *vals, unsigned int rowSize, unsigned int colSize)

75b857d473541532819bd791105cb352c9a43214buliabyak {

75b857d473541532819bd791105cb352c9a43214buliabyak init();

75b857d473541532819bd791105cb352c9a43214buliabyak rows = rowSize;

c53f16f52840e8c0f2be9c1cc3af633c0ba1552emiklosh cols = colSize;

c53f16f52840e8c0f2be9c1cc3af633c0ba1552emiklosh size = rows * cols;

c53f16f52840e8c0f2be9c1cc3af633c0ba1552emiklosh d = new double[size];

c53f16f52840e8c0f2be9c1cc3af633c0ba1552emiklosh for (unsigned int i=0 ; i<size ; i++)

c53f16f52840e8c0f2be9c1cc3af633c0ba1552emiklosh d[i] = vals[i];

1667116521643e2475184b048e0abb77a2aa9735miklosh }

1667116521643e2475184b048e0abb77a2aa9735miklosh

1667116521643e2475184b048e0abb77a2aa9735miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh SVDMatrix(const SVDMatrix &other)

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh {

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh init();

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh assign(other);

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh SVDMatrix &operator=(const SVDMatrix &other)

68664e00e2372534b4df2fdc5f54f836bafece18miklosh {

68664e00e2372534b4df2fdc5f54f836bafece18miklosh assign(other);

68664e00e2372534b4df2fdc5f54f836bafece18miklosh return *this;

68664e00e2372534b4df2fdc5f54f836bafece18miklosh }

68664e00e2372534b4df2fdc5f54f836bafece18miklosh

68664e00e2372534b4df2fdc5f54f836bafece18miklosh virtual ~SVDMatrix()

68664e00e2372534b4df2fdc5f54f836bafece18miklosh {

68664e00e2372534b4df2fdc5f54f836bafece18miklosh delete[] d;

68664e00e2372534b4df2fdc5f54f836bafece18miklosh }

68664e00e2372534b4df2fdc5f54f836bafece18miklosh

c53f16f52840e8c0f2be9c1cc3af633c0ba1552emiklosh double& operator() (unsigned int row, unsigned int col)

c53f16f52840e8c0f2be9c1cc3af633c0ba1552emiklosh {

c53f16f52840e8c0f2be9c1cc3af633c0ba1552emiklosh if (row >= rows || col >= cols)

c53f16f52840e8c0f2be9c1cc3af633c0ba1552emiklosh return badval;

c53f16f52840e8c0f2be9c1cc3af633c0ba1552emiklosh return d[cols*row + col];

c53f16f52840e8c0f2be9c1cc3af633c0ba1552emiklosh }

c53f16f52840e8c0f2be9c1cc3af633c0ba1552emiklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh double operator() (unsigned int row, unsigned int col) const

c53f16f52840e8c0f2be9c1cc3af633c0ba1552emiklosh {

c53f16f52840e8c0f2be9c1cc3af633c0ba1552emiklosh if (row >= rows || col >= cols)

c53f16f52840e8c0f2be9c1cc3af633c0ba1552emiklosh return badval;

c53f16f52840e8c0f2be9c1cc3af633c0ba1552emiklosh return d[cols*row + col];

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh unsigned int getRows()

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh {

c53f16f52840e8c0f2be9c1cc3af633c0ba1552emiklosh return rows;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

c53f16f52840e8c0f2be9c1cc3af633c0ba1552emiklosh unsigned int getCols()

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh return cols;

1667116521643e2475184b048e0abb77a2aa9735miklosh }

1667116521643e2475184b048e0abb77a2aa9735miklosh

1667116521643e2475184b048e0abb77a2aa9735miklosh SVDMatrix multiply(const SVDMatrix &other)

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh {

1667116521643e2475184b048e0abb77a2aa9735miklosh if (cols != other.rows)

1667116521643e2475184b048e0abb77a2aa9735miklosh {

1667116521643e2475184b048e0abb77a2aa9735miklosh SVDMatrix dummy;

1667116521643e2475184b048e0abb77a2aa9735miklosh return dummy;

1667116521643e2475184b048e0abb77a2aa9735miklosh }

1667116521643e2475184b048e0abb77a2aa9735miklosh SVDMatrix result(rows, other.cols);

1667116521643e2475184b048e0abb77a2aa9735miklosh for (unsigned int i=0 ; i<rows ; i++)

1667116521643e2475184b048e0abb77a2aa9735miklosh {

1667116521643e2475184b048e0abb77a2aa9735miklosh for (unsigned int j=0 ; j<other.cols ; j++)

1667116521643e2475184b048e0abb77a2aa9735miklosh {

1667116521643e2475184b048e0abb77a2aa9735miklosh double sum = 0.0;

1667116521643e2475184b048e0abb77a2aa9735miklosh for (unsigned int k=0 ; k<cols ; k++)

1667116521643e2475184b048e0abb77a2aa9735miklosh {

1667116521643e2475184b048e0abb77a2aa9735miklosh //sum += a[i][k] * b[k][j];

1667116521643e2475184b048e0abb77a2aa9735miklosh sum += d[i*cols +k] * other(k, j);

1667116521643e2475184b048e0abb77a2aa9735miklosh }

1667116521643e2475184b048e0abb77a2aa9735miklosh result(i, j) = sum;

1667116521643e2475184b048e0abb77a2aa9735miklosh }

1667116521643e2475184b048e0abb77a2aa9735miklosh

1667116521643e2475184b048e0abb77a2aa9735miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh return result;

c53f16f52840e8c0f2be9c1cc3af633c0ba1552emiklosh }

68664e00e2372534b4df2fdc5f54f836bafece18miklosh

1cda9431ef400135f5e1bd899a94b921bdad0eafmiklosh SVDMatrix transpose()

68664e00e2372534b4df2fdc5f54f836bafece18miklosh {

68664e00e2372534b4df2fdc5f54f836bafece18miklosh SVDMatrix result(cols, rows);

68664e00e2372534b4df2fdc5f54f836bafece18miklosh for (unsigned int i=0 ; i<rows ; i++)

68664e00e2372534b4df2fdc5f54f836bafece18miklosh for (unsigned int j=0 ; j<cols ; j++)

68664e00e2372534b4df2fdc5f54f836bafece18miklosh result(j, i) = d[i*cols + j];

68664e00e2372534b4df2fdc5f54f836bafece18miklosh return result;

68664e00e2372534b4df2fdc5f54f836bafece18miklosh }

68664e00e2372534b4df2fdc5f54f836bafece18miklosh

68664e00e2372534b4df2fdc5f54f836bafece18mikloshprivate:

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh virtual void init()

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh badval = 0.0;

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh d = NULL;

1667116521643e2475184b048e0abb77a2aa9735miklosh rows = 0;

c53f16f52840e8c0f2be9c1cc3af633c0ba1552emiklosh cols = 0;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh size = 0;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh void assign(const SVDMatrix &other)

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh if (d)

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh delete[] d;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh d = 0;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

d9a7c806ee7f408ddb61ff4f233c9d96111ee2b5johanengelen rows = other.rows;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh cols = other.cols;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh size = other.size;

fba63a357654d8b3e84c60007e40aa698cd45d19miklosh d = new double[size];

d9a7c806ee7f408ddb61ff4f233c9d96111ee2b5johanengelen for (unsigned int i=0 ; i<size ; i++)

d9a7c806ee7f408ddb61ff4f233c9d96111ee2b5johanengelen d[i] = other.d[i];

d9a7c806ee7f408ddb61ff4f233c9d96111ee2b5johanengelen }

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh double badval;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh double *d;

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh unsigned int rows;

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh unsigned int cols;

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh unsigned int size;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh};

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9mikloshclass SingularValueDecomposition

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh{

7a7fa095a483e8b652af9f00e5169f62c84f09b9mikloshpublic:

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh /** Construct the singular value decomposition

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh SingularValueDecomposition (const SVDMatrix &mat)

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh A = mat;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh s = NULL;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh s_size = 0;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh calculate();

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh virtual ~SingularValueDecomposition()

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh delete[] s;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh /**

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh SVDMatrix &getU();

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh /**

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh SVDMatrix &getV();

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh /**

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh */ double getS(unsigned int index);

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh /**

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh double norm2();

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh /**

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh double cond();

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh /**

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh int rank();

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9mikloshprivate:

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh void calculate();

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh SVDMatrix A;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh SVDMatrix U;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh double *s;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh unsigned int s_size;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh SVDMatrix V;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh};

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9mikloshstatic double svd_hypot(double a, double b)

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh{

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh double r;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh if (fabs(a) > fabs(b))

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh r = b/a;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh r = fabs(a) * sqrt(1+r*r);

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh else if (b != 0)

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh r = a/b;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh r = fabs(b) * sqrt(1+r*r);

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh else

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh r = 0.0;

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh return r;

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh}

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9mikloshvoid SingularValueDecomposition::calculate()

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh{

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh // Initialize.

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh int m = A.getRows();

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh int n = A.getCols();

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh int nu = (m > n) ? m : n;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh s_size = (m+1 < n) ? m+1 : n;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh s = new double[s_size];

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh U = SVDMatrix(m, nu);

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh V = SVDMatrix(n, n);

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh double *e = new double[n];

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh double *work = new double[m];

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh bool wantu = true;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh bool wantv = true;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh // Reduce A to bidiagonal form, storing the diagonal elements

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh // in s and the super-diagonal elements in e.

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh int nct = (m-1<n) ? m-1 : n;

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh int nrtx = (n-2<m) ? n-2 : m;

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh int nrt = (nrtx>0) ? nrtx : 0;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh for (int k = 0; k < 2; k++) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh if (k < nct) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh // Compute the transformation for the k-th column and

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh // place the k-th diagonal in s[k].

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh // Compute 2-norm of k-th column without under/overflow.

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh s[k] = 0;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh for (int i = k; i < m; i++) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh s[k] = svd_hypot(s[k],A(i, k));

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh if (s[k] != 0.0) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh if (A(k, k) < 0.0) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh s[k] = -s[k];

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh for (int i = k; i < m; i++) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh A(i, k) /= s[k];

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh A(k, k) += 1.0;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh s[k] = -s[k];

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh for (int j = k+1; j < n; j++) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh if ((k < nct) & (s[k] != 0.0)) {

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh // Apply the transformation.

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh double t = 0;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh for (int i = k; i < m; i++) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh t += A(i, k) * A(i, j);

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh t = -t/A(k, k);

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh for (int i = k; i < m; i++) {

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh A(i, j) += t*A(i, k);

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh // Place the k-th row of A into e for the

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh // subsequent calculation of the row transformation.

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh e[j] = A(k, j);

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh if (wantu & (k < nct)) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh // Place the transformation in U for subsequent back

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh // multiplication.

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh for (int i = k; i < m; i++) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh U(i, k) = A(i, k);

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh if (k < nrt) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh // Compute the k-th row transformation and place the

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh // k-th super-diagonal in e[k].

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh // Compute 2-norm without under/overflow.

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh e[k] = 0;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh for (int i = k+1; i < n; i++) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh e[k] = svd_hypot(e[k],e[i]);

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh if (e[k] != 0.0) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh if (e[k+1] < 0.0) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh e[k] = -e[k];

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh for (int i = k+1; i < n; i++) {

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh e[i] /= e[k];

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh }

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh e[k+1] += 1.0;

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh }

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh e[k] = -e[k];

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh if ((k+1 < m) & (e[k] != 0.0)) {

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh // Apply the transformation.

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh for (int i = k+1; i < m; i++) {

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh work[i] = 0.0;

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh }

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh for (int j = k+1; j < n; j++) {

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh for (int i = k+1; i < m; i++) {

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh work[i] += e[j]*A(i, j);

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh }

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh }

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh for (int j = k+1; j < n; j++) {

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh double t = -e[j]/e[k+1];

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh for (int i = k+1; i < m; i++) {

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh A(i, j) += t*work[i];

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh }

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh }

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh }

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh if (wantv) {

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh // Place the transformation in V for subsequent

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh // back multiplication.

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh for (int i = k+1; i < n; i++) {

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh V(i, k) = e[i];

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh }

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh }

1db439af43130c9695dbbb661e893d56006bb072miklosh }

1db439af43130c9695dbbb661e893d56006bb072miklosh }

1db439af43130c9695dbbb661e893d56006bb072miklosh

1db439af43130c9695dbbb661e893d56006bb072miklosh // Set up the final bidiagonal matrix or order p.

1db439af43130c9695dbbb661e893d56006bb072miklosh

1db439af43130c9695dbbb661e893d56006bb072miklosh int p = (n < m+1) ? n : m+1;

1db439af43130c9695dbbb661e893d56006bb072miklosh if (nct < n) {

1db439af43130c9695dbbb661e893d56006bb072miklosh s[nct] = A(nct, nct);

1db439af43130c9695dbbb661e893d56006bb072miklosh }

1db439af43130c9695dbbb661e893d56006bb072miklosh if (m < p) {

1db439af43130c9695dbbb661e893d56006bb072miklosh s[p-1] = 0.0;

1db439af43130c9695dbbb661e893d56006bb072miklosh }

1db439af43130c9695dbbb661e893d56006bb072miklosh if (nrt+1 < p) {

1db439af43130c9695dbbb661e893d56006bb072miklosh e[nrt] = A(nrt, p-1);

1db439af43130c9695dbbb661e893d56006bb072miklosh }

1db439af43130c9695dbbb661e893d56006bb072miklosh e[p-1] = 0.0;

1db439af43130c9695dbbb661e893d56006bb072miklosh

1db439af43130c9695dbbb661e893d56006bb072miklosh // If required, generate U.

1db439af43130c9695dbbb661e893d56006bb072miklosh

1db439af43130c9695dbbb661e893d56006bb072miklosh if (wantu) {

1db439af43130c9695dbbb661e893d56006bb072miklosh for (int j = nct; j < nu; j++) {

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh for (int i = 0; i < m; i++) {

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh U(i, j) = 0.0;

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh }

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh U(j, j) = 1.0;

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh }

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh for (int k = nct-1; k >= 0; k--) {

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh if (s[k] != 0.0) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh for (int j = k+1; j < nu; j++) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh double t = 0;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh for (int i = k; i < m; i++) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh t += U(i, k)*U(i, j);

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh t = -t/U(k, k);

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh for (int i = k; i < m; i++) {

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh U(i, j) += t*U(i, k);

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh for (int i = k; i < m; i++ ) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh U(i, k) = -U(i, k);

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh U(k, k) = 1.0 + U(k, k);

920fbbf5386a5d3b1d0a1304cb71eb48112abe0dmiklosh for (int i = 0; i < k-1; i++) {

920fbbf5386a5d3b1d0a1304cb71eb48112abe0dmiklosh U(i, k) = 0.0;

920fbbf5386a5d3b1d0a1304cb71eb48112abe0dmiklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh } else {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh for (int i = 0; i < m; i++) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh U(i, k) = 0.0;

07b7f1aaaf1087716e784a50cf574a059f7018d3Jon A. Cruz }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh U(k, k) = 1.0;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh // If required, generate V.

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh if (wantv) {

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh for (int k = n-1; k >= 0; k--) {

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh if ((k < nrt) & (e[k] != 0.0)) {

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh for (int j = k+1; j < nu; j++) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh double t = 0;

d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelen for (int i = k+1; i < n; i++) {

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh t += V(i, k)*V(i, j);

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh t = -t/V(k+1, k);

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh for (int i = k+1; i < n; i++) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh V(i, j) += t*V(i, k);

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh for (int i = 0; i < n; i++) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh V(i, k) = 0.0;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh V(k, k) = 1.0;

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh }

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh }

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh // Main iteration loop for the singular values.

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

f9325af537ca5517eb50ef95f432a3204616f6b3apenner int pp = p-1;

f9325af537ca5517eb50ef95f432a3204616f6b3apenner int iter = 0;

f9325af537ca5517eb50ef95f432a3204616f6b3apenner //double eps = pow(2.0,-52.0);

f9325af537ca5517eb50ef95f432a3204616f6b3apenner //double tiny = pow(2.0,-966.0);

f9325af537ca5517eb50ef95f432a3204616f6b3apenner //let's just calculate these now

f9325af537ca5517eb50ef95f432a3204616f6b3apenner //a double can be e � 308.25, so this is safe

f9325af537ca5517eb50ef95f432a3204616f6b3apenner double eps = 2.22e-16;

f9325af537ca5517eb50ef95f432a3204616f6b3apenner double tiny = 1.6e-291;

f9325af537ca5517eb50ef95f432a3204616f6b3apenner while (p > 0) {

f9325af537ca5517eb50ef95f432a3204616f6b3apenner int k,kase;

f9325af537ca5517eb50ef95f432a3204616f6b3apenner

f9325af537ca5517eb50ef95f432a3204616f6b3apenner // Here is where a test for too many iterations would go.

f9325af537ca5517eb50ef95f432a3204616f6b3apenner

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh // This section of the program inspects for

dd6d813ffad339352c39dc0645a792bdd9d8315cmiklosh // negligible elements in the s and e arrays. On

dd6d813ffad339352c39dc0645a792bdd9d8315cmiklosh // completion the variables kase and k are set as follows.

dd6d813ffad339352c39dc0645a792bdd9d8315cmiklosh

dd6d813ffad339352c39dc0645a792bdd9d8315cmiklosh // kase = 1 if s(p) and e[k-1] are negligible and k<p

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh // kase = 2 if s(k) is negligible and k<p

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh // kase = 3 if e[k-1] is negligible, k<p, and

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh // s(k), ..., s(p) are not negligible (qr step).

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh // kase = 4 if e(p-1) is negligible (convergence).

6d9c1b595cc989472a475874007760800d11e6b8pjrm

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh for (k = p-2; k >= -1; k--) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh if (k == -1) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh break;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh if (fabs(e[k]) <=

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh tiny + eps*(fabs(s[k]) + fabs(s[k+1]))) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh e[k] = 0.0;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh break;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh }

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh if (k == p-2) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh kase = 4;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh } else {

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh int ks;

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh for (ks = p-1; ks >= k; ks--) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh if (ks == k) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh break;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh double t = (ks != p ? fabs(e[ks]) : 0.) +

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh (ks != k+1 ? fabs(e[ks-1]) : 0.);

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh if (fabs(s[ks]) <= tiny + eps*t) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh s[ks] = 0.0;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh break;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh if (ks == k) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh kase = 3;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh } else if (ks == p-1) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh kase = 1;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh } else {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh kase = 2;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh k = ks;

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh k++;

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh // Perform the task indicated by kase.

f9325af537ca5517eb50ef95f432a3204616f6b3apenner

f9325af537ca5517eb50ef95f432a3204616f6b3apenner switch (kase) {

f9325af537ca5517eb50ef95f432a3204616f6b3apenner

f9325af537ca5517eb50ef95f432a3204616f6b3apenner // Deflate negligible s(p).

f9325af537ca5517eb50ef95f432a3204616f6b3apenner

f9325af537ca5517eb50ef95f432a3204616f6b3apenner case 1: {

f9325af537ca5517eb50ef95f432a3204616f6b3apenner double f = e[p-2];

f9325af537ca5517eb50ef95f432a3204616f6b3apenner e[p-2] = 0.0;

f9325af537ca5517eb50ef95f432a3204616f6b3apenner for (int j = p-2; j >= k; j--) {

f9325af537ca5517eb50ef95f432a3204616f6b3apenner double t = svd_hypot(s[j],f);

f9325af537ca5517eb50ef95f432a3204616f6b3apenner double cs = s[j]/t;

f9325af537ca5517eb50ef95f432a3204616f6b3apenner double sn = f/t;

f9325af537ca5517eb50ef95f432a3204616f6b3apenner s[j] = t;

f9325af537ca5517eb50ef95f432a3204616f6b3apenner if (j != k) {

f9325af537ca5517eb50ef95f432a3204616f6b3apenner f = -sn*e[j-1];

f9325af537ca5517eb50ef95f432a3204616f6b3apenner e[j-1] = cs*e[j-1];

f9325af537ca5517eb50ef95f432a3204616f6b3apenner }

f9325af537ca5517eb50ef95f432a3204616f6b3apenner if (wantv) {

f9325af537ca5517eb50ef95f432a3204616f6b3apenner for (int i = 0; i < n; i++) {

f9325af537ca5517eb50ef95f432a3204616f6b3apenner t = cs*V(i, j) + sn*V(i, p-1);

f9325af537ca5517eb50ef95f432a3204616f6b3apenner V(i, p-1) = -sn*V(i, j) + cs*V(i, p-1);

f9325af537ca5517eb50ef95f432a3204616f6b3apenner V(i, j) = t;

f9325af537ca5517eb50ef95f432a3204616f6b3apenner }

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh }

f9325af537ca5517eb50ef95f432a3204616f6b3apenner }

52df193b5d1e7bf4f7fbfbb457248c8d13f2007ebuliabyak }

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh break;

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh // Split at negligible s(k).

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh case 2: {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh double f = e[k-1];

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh e[k-1] = 0.0;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh for (int j = k; j < p; j++) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh double t = svd_hypot(s[j],f);

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh double cs = s[j]/t;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh double sn = f/t;

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh s[j] = t;

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh f = -sn*e[j];

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh e[j] = cs*e[j];

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh if (wantu) {

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh for (int i = 0; i < m; i++) {

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh t = cs*U(i, j) + sn*U(i, k-1);

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh U(i, k-1) = -sn*U(i, j) + cs*U(i, k-1);

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh U(i, j) = t;

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh }

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh }

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh }

d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelen }

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh break;

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh // Perform one qr step.

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh case 3: {

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh // Calculate the shift.

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh double scale = fabs(s[p-1]);

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh double d = fabs(s[p-2]);

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh if (d>scale) scale=d;

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh d = fabs(e[p-2]);

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh if (d>scale) scale=d;

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh d = fabs(s[k]);

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh if (d>scale) scale=d;

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh d = fabs(e[k]);

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh if (d>scale) scale=d;

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh double sp = s[p-1]/scale;

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh double spm1 = s[p-2]/scale;

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh double epm1 = e[p-2]/scale;

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh double sk = s[k]/scale;

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh double ek = e[k]/scale;

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh double b = ((spm1 + sp)*(spm1 - sp) + epm1*epm1)/2.0;

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh double c = (sp*epm1)*(sp*epm1);

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh double shift = 0.0;

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh if ((b != 0.0) | (c != 0.0)) {

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh shift = sqrt(b*b + c);

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh if (b < 0.0) {

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh shift = -shift;

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh }

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh shift = c/(b + shift);

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh }

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh double f = (sk + sp)*(sk - sp) + shift;

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh double g = sk*ek;

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh // Chase zeros.

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh for (int j = k; j < p-1; j++) {

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh double t = svd_hypot(f,g);

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh double cs = f/t;

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh double sn = g/t;

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh if (j != k) {

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh e[j-1] = t;

07b7f1aaaf1087716e784a50cf574a059f7018d3Jon A. Cruz }

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh f = cs*s[j] + sn*e[j];

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh e[j] = cs*e[j] - sn*s[j];

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh g = sn*s[j+1];

dc4f69a188c203f2fdc65f22d0d57904a8c52dd7miklosh s[j+1] = cs*s[j+1];

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh if (wantv) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh for (int i = 0; i < n; i++) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh t = cs*V(i, j) + sn*V(i, j+1);

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh V(i, j+1) = -sn*V(i, j) + cs*V(i, j+1);

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh V(i, j) = t;

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh }

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh t = svd_hypot(f,g);

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh cs = f/t;

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh sn = g/t;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh s[j] = t;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh f = cs*e[j] + sn*s[j+1];

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh s[j+1] = -sn*e[j] + cs*s[j+1];

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh g = sn*e[j+1];

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh e[j+1] = cs*e[j+1];

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh if (wantu && (j < m-1)) {

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh for (int i = 0; i < m; i++) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh t = cs*U(i, j) + sn*U(i, j+1);

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh U(i, j+1) = -sn*U(i, j) + cs*U(i, j+1);

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh U(i, j) = t;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh e[p-2] = f;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh iter = iter + 1;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh break;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh // Convergence.

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh case 4: {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh // Make the singular values positive.

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh if (s[k] <= 0.0) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh s[k] = (s[k] < 0.0 ? -s[k] : 0.0);

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh if (wantv) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh for (int i = 0; i <= pp; i++) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh V(i, k) = -V(i, k);

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh // Order the singular values.

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh while (k < pp) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh if (s[k] >= s[k+1]) {

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh break;

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh }

d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelen double t = s[k];

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh s[k] = s[k+1];

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh s[k+1] = t;

d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelen if (wantv && (k < n-1)) {

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh for (int i = 0; i < n; i++) {

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh t = V(i, k+1); V(i, k+1) = V(i, k); V(i, k) = t;

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh }

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh }

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh if (wantu && (k < m-1)) {

3686c32a570c3df738a09b34e85fc5d6bd50d020miklosh for (int i = 0; i < m; i++) {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh t = U(i, k+1); U(i, k+1) = U(i, k); U(i, k) = t;

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh k++;

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh }

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh iter = 0;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh p--;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh break;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

07b7f1aaaf1087716e784a50cf574a059f7018d3Jon A. Cruz }

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh delete e;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh delete work;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh}

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

d27f5758e12c3107ee69e66702043931e0756f6bmikloshSVDMatrix &SingularValueDecomposition::getU()

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh{

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh return U;

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh}

c53f16f52840e8c0f2be9c1cc3af633c0ba1552emiklosh

0e14f9e966c4b6012538d30cd0db7a775b879760JucaBlues

c53f16f52840e8c0f2be9c1cc3af633c0ba1552emikloshSVDMatrix &SingularValueDecomposition::getV()

c53f16f52840e8c0f2be9c1cc3af633c0ba1552emiklosh{

c53f16f52840e8c0f2be9c1cc3af633c0ba1552emiklosh return V;

c53f16f52840e8c0f2be9c1cc3af633c0ba1552emiklosh}

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh

d27f5758e12c3107ee69e66702043931e0756f6bmikloshdouble SingularValueDecomposition::getS(unsigned int index)

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh{

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh if (index >= s_size)

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh return 0.0;

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh return s[index];

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh}

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

d27f5758e12c3107ee69e66702043931e0756f6bmikloshdouble SingularValueDecomposition::norm2()

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh{

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh return s[0];

8ec52d39f409ecf67125cd4d878b844e9391e7e8miklosh}

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9mikloshdouble SingularValueDecomposition::cond()

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh{

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh return s[0]/s[2];

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh}

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh

d27f5758e12c3107ee69e66702043931e0756f6bmikloshint SingularValueDecomposition::rank()

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh{

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh double eps = pow(2.0,-52.0);

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh double tol = 3.0*s[0]*eps;

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh int r = 0;

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh for (int i = 0; i < 3; i++)

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh {

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh if (s[i] > tol)

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh r++;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh }

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh return r;

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh}

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh#define pi 3.14159

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh#define pxToCm 0.03

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh#define piToRad 0.0174532925

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh#define docHeightCm 22.86

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh

d27f5758e12c3107ee69e66702043931e0756f6bmiklosh

7a7fa095a483e8b652af9f00e5169f62c84f09b9mikloshstatic Glib::ustring getAttribute( Inkscape::XML::Node *node, char const *attrName)

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh{

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh Glib::ustring val;

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh char const *valstr = node->attribute(attrName);

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh if (valstr)

7a7fa095a483e8b652af9f00e5169f62c84f09b9miklosh val = valstr;

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh return val;

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh}

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh

3711b3e25395437ee0a09dbbb2a76d999c4ef322mikloshstatic Glib::ustring getExtension(const Glib::ustring &fname)

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh{

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh Glib::ustring ext;

75b857d473541532819bd791105cb352c9a43214buliabyak

75b857d473541532819bd791105cb352c9a43214buliabyak std::string::size_type pos = fname.rfind('.');

75b857d473541532819bd791105cb352c9a43214buliabyak if (pos == fname.npos)

75b857d473541532819bd791105cb352c9a43214buliabyak {

75b857d473541532819bd791105cb352c9a43214buliabyak ext = "";

75b857d473541532819bd791105cb352c9a43214buliabyak }

75b857d473541532819bd791105cb352c9a43214buliabyak else

75b857d473541532819bd791105cb352c9a43214buliabyak {

75b857d473541532819bd791105cb352c9a43214buliabyak ext = fname.substr(pos);

75b857d473541532819bd791105cb352c9a43214buliabyak }

75b857d473541532819bd791105cb352c9a43214buliabyak return ext;

75b857d473541532819bd791105cb352c9a43214buliabyak}

75b857d473541532819bd791105cb352c9a43214buliabyak

75b857d473541532819bd791105cb352c9a43214buliabyak

75b857d473541532819bd791105cb352c9a43214buliabyakstatic Glib::ustring formatTransform(NR::Matrix &tf)

75b857d473541532819bd791105cb352c9a43214buliabyak{

75b857d473541532819bd791105cb352c9a43214buliabyak Glib::ustring str;

75b857d473541532819bd791105cb352c9a43214buliabyak if (!tf.test_identity())

75b857d473541532819bd791105cb352c9a43214buliabyak {

75b857d473541532819bd791105cb352c9a43214buliabyak StringOutputStream outs;

75b857d473541532819bd791105cb352c9a43214buliabyak OutputStreamWriter out(outs);

75b857d473541532819bd791105cb352c9a43214buliabyak out.printf("matrix(%.3f %.3f %.3f %.3f %.3f %.3f)",

75b857d473541532819bd791105cb352c9a43214buliabyak tf[0], tf[1], tf[2], tf[3], tf[4], tf[5]);

75b857d473541532819bd791105cb352c9a43214buliabyak str = outs.getString();

75b857d473541532819bd791105cb352c9a43214buliabyak }

75b857d473541532819bd791105cb352c9a43214buliabyak return str;

75b857d473541532819bd791105cb352c9a43214buliabyak}

75b857d473541532819bd791105cb352c9a43214buliabyak

75b857d473541532819bd791105cb352c9a43214buliabyak

75b857d473541532819bd791105cb352c9a43214buliabyak

75b857d473541532819bd791105cb352c9a43214buliabyak

75b857d473541532819bd791105cb352c9a43214buliabyak

75b857d473541532819bd791105cb352c9a43214buliabyakstatic NR::Matrix getODFTransform(const SPItem *item)

75b857d473541532819bd791105cb352c9a43214buliabyak{

75b857d473541532819bd791105cb352c9a43214buliabyak //### Get SVG-to-ODF transform

75b857d473541532819bd791105cb352c9a43214buliabyak NR::Matrix tf;

75b857d473541532819bd791105cb352c9a43214buliabyak tf = sp_item_i2d_affine(item);

75b857d473541532819bd791105cb352c9a43214buliabyak //Flip Y into document coordinates

75b857d473541532819bd791105cb352c9a43214buliabyak double doc_height = sp_document_height(SP_ACTIVE_DOCUMENT);

75b857d473541532819bd791105cb352c9a43214buliabyak NR::Matrix doc2dt_tf = NR::Matrix(NR::scale(1.0, -1.0));

75b857d473541532819bd791105cb352c9a43214buliabyak doc2dt_tf = doc2dt_tf * NR::Matrix(NR::translate(0, doc_height));

75b857d473541532819bd791105cb352c9a43214buliabyak tf = tf * doc2dt_tf;

75b857d473541532819bd791105cb352c9a43214buliabyak tf = tf * NR::Matrix(NR::scale(pxToCm));

75b857d473541532819bd791105cb352c9a43214buliabyak return tf;

75b857d473541532819bd791105cb352c9a43214buliabyak}

75b857d473541532819bd791105cb352c9a43214buliabyak

75b857d473541532819bd791105cb352c9a43214buliabyak

75b857d473541532819bd791105cb352c9a43214buliabyak

75b857d473541532819bd791105cb352c9a43214buliabyak

75b857d473541532819bd791105cb352c9a43214buliabyakstatic NR::Maybe<NR::Rect> getODFBoundingBox(const SPItem *item)

75b857d473541532819bd791105cb352c9a43214buliabyak{

75b857d473541532819bd791105cb352c9a43214buliabyak NR::Maybe<NR::Rect> bbox = sp_item_bbox_desktop((SPItem *)item);

75b857d473541532819bd791105cb352c9a43214buliabyak if (bbox) {

c4723fe0caa2096d00cb31a7d1506351ba8102dbmiklosh double doc_height = sp_document_height(SP_ACTIVE_DOCUMENT);

c4723fe0caa2096d00cb31a7d1506351ba8102dbmiklosh NR::Matrix doc2dt_tf = NR::Matrix(NR::scale(1.0, -1.0));

c4723fe0caa2096d00cb31a7d1506351ba8102dbmiklosh doc2dt_tf = doc2dt_tf * NR::Matrix(NR::translate(0, doc_height));

c4723fe0caa2096d00cb31a7d1506351ba8102dbmiklosh bbox = *bbox * doc2dt_tf;

0e14f9e966c4b6012538d30cd0db7a775b879760JucaBlues bbox = *bbox * NR::Matrix(NR::scale(pxToCm));

0e14f9e966c4b6012538d30cd0db7a775b879760JucaBlues }

0e14f9e966c4b6012538d30cd0db7a775b879760JucaBlues return bbox;

0e14f9e966c4b6012538d30cd0db7a775b879760JucaBlues}

0e14f9e966c4b6012538d30cd0db7a775b879760JucaBlues

0e14f9e966c4b6012538d30cd0db7a775b879760JucaBlues

0e14f9e966c4b6012538d30cd0db7a775b879760JucaBlues

0e14f9e966c4b6012538d30cd0db7a775b879760JucaBluesstatic NR::Matrix getODFItemTransform(const SPItem *item)

0e14f9e966c4b6012538d30cd0db7a775b879760JucaBlues{

0e14f9e966c4b6012538d30cd0db7a775b879760JucaBlues NR::Matrix itemTransform = NR::Matrix(NR::scale(1, -1));

0e14f9e966c4b6012538d30cd0db7a775b879760JucaBlues itemTransform = itemTransform * item->transform;

0e14f9e966c4b6012538d30cd0db7a775b879760JucaBlues itemTransform = itemTransform * NR::Matrix(NR::scale(1, -1));

c4723fe0caa2096d00cb31a7d1506351ba8102dbmiklosh return itemTransform;

c4723fe0caa2096d00cb31a7d1506351ba8102dbmiklosh}

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh

3711b3e25395437ee0a09dbbb2a76d999c4ef322mikloshstatic void analyzeTransform(NR::Matrix &tf,

03e63790ef0fa2919fc5f9f3e0d018adf317919dmiklosh double &rotate, double &/*xskew*/, double &/*yskew*/,

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh double &xscale, double &yscale)

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh{

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh SVDMatrix mat(2, 2);

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh mat(0, 0) = tf[0];

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh mat(0, 1) = tf[1];

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh mat(1, 0) = tf[2];

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh mat(1, 1) = tf[3];

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh SingularValueDecomposition svd(mat);

082cac3fb4c7fc60e92bbd4de269550801dcea6bmiklosh

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh SVDMatrix U = svd.getU();

082cac3fb4c7fc60e92bbd4de269550801dcea6bmiklosh SVDMatrix V = svd.getV();

0e14f9e966c4b6012538d30cd0db7a775b879760JucaBlues SVDMatrix Vt = V.transpose();

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh SVDMatrix UVt = U.multiply(Vt);

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh double s0 = svd.getS(0);

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh double s1 = svd.getS(1);

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh xscale = s0;

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh yscale = s1;

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh //g_message("## s0:%.3f s1:%.3f", s0, s1);

c4723fe0caa2096d00cb31a7d1506351ba8102dbmiklosh //g_message("## u:%.3f %.3f %.3f %.3f", U(0,0), U(0,1), U(1,0), U(1,1));

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh //g_message("## v:%.3f %.3f %.3f %.3f", V(0,0), V(0,1), V(1,0), V(1,1));

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh //g_message("## vt:%.3f %.3f %.3f %.3f", Vt(0,0), Vt(0,1), Vt(1,0), Vt(1,1));

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh //g_message("## uvt:%.3f %.3f %.3f %.3f", UVt(0,0), UVt(0,1), UVt(1,0), UVt(1,1));

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh rotate = UVt(0,0);

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh}

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh

77364929ced3ec0bc5c9f47440606615c559084emiklosh

77364929ced3ec0bc5c9f47440606615c559084emikloshstatic void gatherText(Inkscape::XML::Node *node, Glib::ustring &buf)

77364929ced3ec0bc5c9f47440606615c559084emiklosh{

77364929ced3ec0bc5c9f47440606615c559084emiklosh if (node->type() == Inkscape::XML::TEXT_NODE)

c4723fe0caa2096d00cb31a7d1506351ba8102dbmiklosh {

77364929ced3ec0bc5c9f47440606615c559084emiklosh char *s = (char *)node->content();

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh if (s)

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh buf.append(s);

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh }

75b857d473541532819bd791105cb352c9a43214buliabyak

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh for (Inkscape::XML::Node *child = node->firstChild() ;

75b857d473541532819bd791105cb352c9a43214buliabyak child != NULL; child = child->next())

75b857d473541532819bd791105cb352c9a43214buliabyak {

75b857d473541532819bd791105cb352c9a43214buliabyak gatherText(child, buf);

75b857d473541532819bd791105cb352c9a43214buliabyak }

75b857d473541532819bd791105cb352c9a43214buliabyak

75b857d473541532819bd791105cb352c9a43214buliabyak}

75b857d473541532819bd791105cb352c9a43214buliabyak

3711b3e25395437ee0a09dbbb2a76d999c4ef322mikloshOdfOutput::preprocess(ZipFile &zf, Inkscape::XML::Node *node)

c4723fe0caa2096d00cb31a7d1506351ba8102dbmiklosh{

c4723fe0caa2096d00cb31a7d1506351ba8102dbmiklosh

c4723fe0caa2096d00cb31a7d1506351ba8102dbmiklosh Glib::ustring nodeName = node->name();

c4723fe0caa2096d00cb31a7d1506351ba8102dbmiklosh Glib::ustring id = getAttribute(node, "id");

c4723fe0caa2096d00cb31a7d1506351ba8102dbmiklosh

c4723fe0caa2096d00cb31a7d1506351ba8102dbmiklosh //### First, check for metadata

c4723fe0caa2096d00cb31a7d1506351ba8102dbmiklosh if (nodeName == "metadata" || nodeName == "svg:metadata")

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh {

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh Inkscape::XML::Node *mchild = node->firstChild() ;

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh if (!mchild || strcmp(mchild->name(), "rdf:RDF"))

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh return;

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh Inkscape::XML::Node *rchild = mchild->firstChild() ;

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh if (!rchild || strcmp(rchild->name(), "cc:Work"))

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh return;

7ec85862d9730e449ed5c2a86201bc9ca1daa0aamiklosh for (Inkscape::XML::Node *cchild = rchild->firstChild() ;

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh cchild ; cchild = cchild->next())

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh {

0e14f9e966c4b6012538d30cd0db7a775b879760JucaBlues Glib::ustring ccName = cchild->name();

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh Glib::ustring ccVal;

0e14f9e966c4b6012538d30cd0db7a775b879760JucaBlues gatherText(cchild, ccVal);

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh //g_message("ccName: %s ccVal:%s", ccName.c_str(), ccVal.c_str());

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh metadata[ccName] = ccVal;

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh }

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh return;

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh }

c4723fe0caa2096d00cb31a7d1506351ba8102dbmiklosh

c4723fe0caa2096d00cb31a7d1506351ba8102dbmiklosh //Now consider items.

c4723fe0caa2096d00cb31a7d1506351ba8102dbmiklosh SPObject *reprobj = SP_ACTIVE_DOCUMENT->getObjectByRepr(node);

c4723fe0caa2096d00cb31a7d1506351ba8102dbmiklosh if (!reprobj)

c4723fe0caa2096d00cb31a7d1506351ba8102dbmiklosh return;

7ec85862d9730e449ed5c2a86201bc9ca1daa0aamiklosh if (!SP_IS_ITEM(reprobj))

c4723fe0caa2096d00cb31a7d1506351ba8102dbmiklosh {

c4723fe0caa2096d00cb31a7d1506351ba8102dbmiklosh return;

7ec85862d9730e449ed5c2a86201bc9ca1daa0aamiklosh }

0e14f9e966c4b6012538d30cd0db7a775b879760JucaBlues SPItem *item = SP_ITEM(reprobj);

7ec85862d9730e449ed5c2a86201bc9ca1daa0aamiklosh //### Get SVG-to-ODF transform

7ec85862d9730e449ed5c2a86201bc9ca1daa0aamiklosh NR::Matrix tf = getODFTransform(item);

7ec85862d9730e449ed5c2a86201bc9ca1daa0aamiklosh

c4723fe0caa2096d00cb31a7d1506351ba8102dbmiklosh if (nodeName == "image" || nodeName == "svg:image")

c4723fe0caa2096d00cb31a7d1506351ba8102dbmiklosh {

c4723fe0caa2096d00cb31a7d1506351ba8102dbmiklosh //g_message("image");

c4723fe0caa2096d00cb31a7d1506351ba8102dbmiklosh Glib::ustring href = getAttribute(node, "xlink:href");

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh if (href.size() > 0)

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh {

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh Glib::ustring oldName = href;

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh Glib::ustring ext = getExtension(oldName);

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh if (ext == ".jpeg")

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh ext = ".jpg";

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh if (imageTable.find(oldName) == imageTable.end())

0e14f9e966c4b6012538d30cd0db7a775b879760JucaBlues {

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh char buf[64];

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh snprintf(buf, 63, "Pictures/image%d%s",

0e14f9e966c4b6012538d30cd0db7a775b879760JucaBlues (int)imageTable.size(), ext.c_str());

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh Glib::ustring newName = buf;

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh imageTable[oldName] = newName;

0e14f9e966c4b6012538d30cd0db7a775b879760JucaBlues Glib::ustring comment = "old name was: ";

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh comment.append(oldName);

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh URI oldUri(oldName);

0e14f9e966c4b6012538d30cd0db7a775b879760JucaBlues //g_message("oldpath:%s", oldUri.getNativePath().c_str());

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh //# if relative to the documentURI, get proper path

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh URI resUri = documentUri.resolve(oldUri);

0e14f9e966c4b6012538d30cd0db7a775b879760JucaBlues DOMString pathName = resUri.getNativePath();

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh //g_message("native path:%s", pathName.c_str());

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh ZipEntry *ze = zf.addFile(pathName, comment);

0e14f9e966c4b6012538d30cd0db7a775b879760JucaBlues if (ze)

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh {

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh ze->setFileName(newName);

0e14f9e966c4b6012538d30cd0db7a775b879760JucaBlues }

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh else

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh {

0e14f9e966c4b6012538d30cd0db7a775b879760JucaBlues g_warning("Could not load image file '%s'", pathName.c_str());

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh }

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh }

0e14f9e966c4b6012538d30cd0db7a775b879760JucaBlues }

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh }

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh for (Inkscape::XML::Node *child = node->firstChild() ;

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh child ; child = child->next())

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh preprocess(zf, child);

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh}

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh

fd39535b3a5276f8962a3f99072668f3e63421edmikloshbool OdfOutput::writeManifest(ZipFile &zf)

fd39535b3a5276f8962a3f99072668f3e63421edmiklosh{

fd39535b3a5276f8962a3f99072668f3e63421edmiklosh BufferOutputStream bouts;

fd39535b3a5276f8962a3f99072668f3e63421edmiklosh OutputStreamWriter outs(bouts);

77364929ced3ec0bc5c9f47440606615c559084emiklosh

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh time_t tim;

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh time(&tim);

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf("<?xml version=\"1.0\" encoding=\"UTF-8\"?>\n");

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf("<!DOCTYPE manifest:manifest PUBLIC \"-//OpenOffice.org//DTD Manifest 1.0//EN\" \"Manifest.dtd\">\n");

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf("\n");

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf("\n");

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf("<!--\n");

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf("*************************************************************************\n");

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf(" file: manifest.xml\n");

7ec85862d9730e449ed5c2a86201bc9ca1daa0aamiklosh outs.printf(" Generated by Inkscape: %s", ctime(&tim)); //ctime has its own <cr>

7ec85862d9730e449ed5c2a86201bc9ca1daa0aamiklosh outs.printf(" http://www.inkscape.org\n");

7ec85862d9730e449ed5c2a86201bc9ca1daa0aamiklosh outs.printf("*************************************************************************\n");

e45563a3c46261d8c32014f8e516857ba01bd7b7miklosh outs.printf("-->\n");

e45563a3c46261d8c32014f8e516857ba01bd7b7miklosh outs.printf("\n");

e45563a3c46261d8c32014f8e516857ba01bd7b7miklosh outs.printf("\n");

e45563a3c46261d8c32014f8e516857ba01bd7b7miklosh outs.printf("<manifest:manifest xmlns:manifest=\"urn:oasis:names:tc:opendocument:xmlns:manifest:1.0\">\n");

e45563a3c46261d8c32014f8e516857ba01bd7b7miklosh outs.printf(" <manifest:file-entry manifest:media-type=\"application/vnd.oasis.opendocument.graphics\" manifest:full-path=\"/\"/>\n");

e45563a3c46261d8c32014f8e516857ba01bd7b7miklosh outs.printf(" <manifest:file-entry manifest:media-type=\"text/xml\" manifest:full-path=\"content.xml\"/>\n");

e45563a3c46261d8c32014f8e516857ba01bd7b7miklosh outs.printf(" <manifest:file-entry manifest:media-type=\"text/xml\" manifest:full-path=\"styles.xml\"/>\n");

e45563a3c46261d8c32014f8e516857ba01bd7b7miklosh outs.printf(" <manifest:file-entry manifest:media-type=\"text/xml\" manifest:full-path=\"meta.xml\"/>\n");

e45563a3c46261d8c32014f8e516857ba01bd7b7miklosh outs.printf(" <!--List our images here-->\n");

e45563a3c46261d8c32014f8e516857ba01bd7b7miklosh std::map<Glib::ustring, Glib::ustring>::iterator iter;

e45563a3c46261d8c32014f8e516857ba01bd7b7miklosh for (iter = imageTable.begin() ; iter!=imageTable.end() ; iter++)

e45563a3c46261d8c32014f8e516857ba01bd7b7miklosh {

e45563a3c46261d8c32014f8e516857ba01bd7b7miklosh Glib::ustring oldName = iter->first;

e45563a3c46261d8c32014f8e516857ba01bd7b7miklosh Glib::ustring newName = iter->second;

e45563a3c46261d8c32014f8e516857ba01bd7b7miklosh

e45563a3c46261d8c32014f8e516857ba01bd7b7miklosh Glib::ustring ext = getExtension(oldName);

17d87f5698f5c2958d38c6a6207c7b322a7adaf9johanengelen if (ext == ".jpeg")

e45563a3c46261d8c32014f8e516857ba01bd7b7miklosh ext = ".jpg";

e45563a3c46261d8c32014f8e516857ba01bd7b7miklosh outs.printf(" <manifest:file-entry manifest:media-type=\"");

e45563a3c46261d8c32014f8e516857ba01bd7b7miklosh if (ext == ".gif")

7ec85862d9730e449ed5c2a86201bc9ca1daa0aamiklosh outs.printf("image/gif");

7ec85862d9730e449ed5c2a86201bc9ca1daa0aamiklosh else if (ext == ".png")

7ec85862d9730e449ed5c2a86201bc9ca1daa0aamiklosh outs.printf("image/png");

7ec85862d9730e449ed5c2a86201bc9ca1daa0aamiklosh else if (ext == ".jpg")

7ec85862d9730e449ed5c2a86201bc9ca1daa0aamiklosh outs.printf("image/jpeg");

03e63790ef0fa2919fc5f9f3e0d018adf317919dmiklosh outs.printf("\" manifest:full-path=\"");

03e63790ef0fa2919fc5f9f3e0d018adf317919dmiklosh outs.printf(newName.c_str());

03e63790ef0fa2919fc5f9f3e0d018adf317919dmiklosh outs.printf("\"/>\n");

03e63790ef0fa2919fc5f9f3e0d018adf317919dmiklosh }

03e63790ef0fa2919fc5f9f3e0d018adf317919dmiklosh outs.printf("</manifest:manifest>\n");

03e63790ef0fa2919fc5f9f3e0d018adf317919dmiklosh

03e63790ef0fa2919fc5f9f3e0d018adf317919dmiklosh outs.close();

03e63790ef0fa2919fc5f9f3e0d018adf317919dmiklosh

03e63790ef0fa2919fc5f9f3e0d018adf317919dmiklosh //Make our entry

03e63790ef0fa2919fc5f9f3e0d018adf317919dmiklosh ZipEntry *ze = zf.newEntry("META-INF/manifest.xml", "ODF file manifest");

03e63790ef0fa2919fc5f9f3e0d018adf317919dmiklosh ze->setUncompressedData(bouts.getBuffer());

d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelen ze->finish();

03e63790ef0fa2919fc5f9f3e0d018adf317919dmiklosh

03e63790ef0fa2919fc5f9f3e0d018adf317919dmiklosh return true;

03e63790ef0fa2919fc5f9f3e0d018adf317919dmiklosh}

03e63790ef0fa2919fc5f9f3e0d018adf317919dmiklosh

03e63790ef0fa2919fc5f9f3e0d018adf317919dmiklosh

03e63790ef0fa2919fc5f9f3e0d018adf317919dmikloshbool OdfOutput::writeMeta(ZipFile &zf)

03e63790ef0fa2919fc5f9f3e0d018adf317919dmiklosh{

03e63790ef0fa2919fc5f9f3e0d018adf317919dmiklosh BufferOutputStream bouts;

03e63790ef0fa2919fc5f9f3e0d018adf317919dmiklosh OutputStreamWriter outs(bouts);

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh time_t tim;

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh time(&tim);

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh std::map<Glib::ustring, Glib::ustring>::iterator iter;

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh Glib::ustring creator = "unknown";

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh iter = metadata.find("dc:creator");

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh if (iter != metadata.end())

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh creator = iter->second;

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh Glib::ustring date = "";

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh iter = metadata.find("dc:date");

7ec85862d9730e449ed5c2a86201bc9ca1daa0aamiklosh if (iter != metadata.end())

7ec85862d9730e449ed5c2a86201bc9ca1daa0aamiklosh date = iter->second;

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf("<?xml version=\"1.0\" encoding=\"UTF-8\"?>\n");

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf("\n");

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf("\n");

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf("<!--\n");

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf("*************************************************************************\n");

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf(" file: meta.xml\n");

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf(" Generated by Inkscape: %s", ctime(&tim)); //ctime has its own <cr>

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf(" http://www.inkscape.org\n");

d37634d73670180f99a3e0ea583621373d90ec4fJohan Engelen outs.printf("*************************************************************************\n");

7ec85862d9730e449ed5c2a86201bc9ca1daa0aamiklosh outs.printf("-->\n");

7ec85862d9730e449ed5c2a86201bc9ca1daa0aamiklosh outs.printf("\n");

7ec85862d9730e449ed5c2a86201bc9ca1daa0aamiklosh outs.printf("\n");

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf("<office:document-meta\n");

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf("xmlns:office=\"urn:oasis:names:tc:opendocument:xmlns:office:1.0\"\n");

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf("xmlns:xlink=\"http://www.w3.org/1999/xlink\"\n");

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf("xmlns:dc=\"http://purl.org/dc/elements/1.1/\"\n");

522aa9b8f493ba0c8e8b0bb536a563c96f5430a8miklosh outs.printf("xmlns:meta=\"urn:oasis:names:tc:opendocument:xmlns:meta:1.0\"\n");

17d87f5698f5c2958d38c6a6207c7b322a7adaf9johanengelen outs.printf("xmlns:presentation=\"urn:oasis:names:tc:opendocument:xmlns:presentation:1.0\"\n");

7ec85862d9730e449ed5c2a86201bc9ca1daa0aamiklosh outs.printf("xmlns:ooo=\"http://openoffice.org/2004/office\"\n");

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf("xmlns:smil=\"urn:oasis:names:tc:opendocument:xmlns:smil-compatible:1.0\"\n");

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf("xmlns:anim=\"urn:oasis:names:tc:opendocument:xmlns:animation:1.0\"\n");

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf("office:version=\"1.0\">\n");

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf("<office:meta>\n");

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf(" <meta:generator>Inkscape.org - 0.45</meta:generator>\n");

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf(" <meta:initial-creator>%#s</meta:initial-creator>\n",

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh creator.c_str());

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf(" <meta:creation-date>%#s</meta:creation-date>\n", date.c_str());

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh for (iter = metadata.begin() ; iter != metadata.end() ; iter++)

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh {

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh Glib::ustring name = iter->first;

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh Glib::ustring value = iter->second;

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh if (name.size() > 0 && value.size()>0)

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh {

7ec85862d9730e449ed5c2a86201bc9ca1daa0aamiklosh outs.printf(" <%#s>%#s</%#s>\n",

e45563a3c46261d8c32014f8e516857ba01bd7b7miklosh name.c_str(), value.c_str(), name.c_str());

7ec85862d9730e449ed5c2a86201bc9ca1daa0aamiklosh }

e45563a3c46261d8c32014f8e516857ba01bd7b7miklosh }

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf(" <meta:editing-cycles>2</meta:editing-cycles>\n");

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf(" <meta:editing-duration>PT56S</meta:editing-duration>\n");

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf(" <meta:user-defined meta:name=\"Info 1\"/>\n");

7ec85862d9730e449ed5c2a86201bc9ca1daa0aamiklosh outs.printf(" <meta:user-defined meta:name=\"Info 2\"/>\n");

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf(" <meta:user-defined meta:name=\"Info 3\"/>\n");

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf(" <meta:user-defined meta:name=\"Info 4\"/>\n");

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf(" <meta:document-statistic meta:object-count=\"2\"/>\n");

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf("</office:meta>\n");

522aa9b8f493ba0c8e8b0bb536a563c96f5430a8miklosh outs.printf("</office:document-meta>\n");

522aa9b8f493ba0c8e8b0bb536a563c96f5430a8miklosh outs.printf("\n");

522aa9b8f493ba0c8e8b0bb536a563c96f5430a8miklosh outs.printf("\n");

522aa9b8f493ba0c8e8b0bb536a563c96f5430a8miklosh

522aa9b8f493ba0c8e8b0bb536a563c96f5430a8miklosh

522aa9b8f493ba0c8e8b0bb536a563c96f5430a8miklosh outs.close();

522aa9b8f493ba0c8e8b0bb536a563c96f5430a8miklosh

522aa9b8f493ba0c8e8b0bb536a563c96f5430a8miklosh //Make our entry

522aa9b8f493ba0c8e8b0bb536a563c96f5430a8miklosh ZipEntry *ze = zf.newEntry("meta.xml", "ODF info file");

522aa9b8f493ba0c8e8b0bb536a563c96f5430a8miklosh ze->setUncompressedData(bouts.getBuffer());

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh ze->finish();

7ec85862d9730e449ed5c2a86201bc9ca1daa0aamiklosh

7ec85862d9730e449ed5c2a86201bc9ca1daa0aamiklosh return true;

7ec85862d9730e449ed5c2a86201bc9ca1daa0aamiklosh}

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh

3711b3e25395437ee0a09dbbb2a76d999c4ef322mikloshbool OdfOutput::writeStyle(ZipFile &zf)

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh{

7ec85862d9730e449ed5c2a86201bc9ca1daa0aamiklosh BufferOutputStream bouts;

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh OutputStreamWriter outs(bouts);

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh /*

7ec85862d9730e449ed5c2a86201bc9ca1daa0aamiklosh outs.printf("<!-- ####### Styles from Inkscape document ####### -->\n");

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh std::vector<StyleInfo>::iterator iter;

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh for (iter = styleTable.begin() ; iter != styleTable.end() ; iter++)

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh {

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf("<style:style style:name=\"%s\"", iter->name.c_str());

7ec85862d9730e449ed5c2a86201bc9ca1daa0aamiklosh StyleInfo s(*iter);

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf(" style:family=\"graphic\" style:parent-style-name=\"standard\">\n");

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf(" <style:graphic-properties");

522aa9b8f493ba0c8e8b0bb536a563c96f5430a8miklosh outs.printf(" draw:fill=\"%s\" ", s.fill.c_str());

522aa9b8f493ba0c8e8b0bb536a563c96f5430a8miklosh if (s.fill != "none")

522aa9b8f493ba0c8e8b0bb536a563c96f5430a8miklosh {

522aa9b8f493ba0c8e8b0bb536a563c96f5430a8miklosh outs.printf(" draw:fill-color=\"%s\" ", s.fillColor.c_str());

522aa9b8f493ba0c8e8b0bb536a563c96f5430a8miklosh outs.printf(" draw:fill-opacity=\"%s\" ", s.fillOpacity.c_str());

522aa9b8f493ba0c8e8b0bb536a563c96f5430a8miklosh }

522aa9b8f493ba0c8e8b0bb536a563c96f5430a8miklosh outs.printf(" draw:stroke=\"%s\" ", s.stroke.c_str());

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh if (s.stroke != "none")

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh {

17d87f5698f5c2958d38c6a6207c7b322a7adaf9johanengelen outs.printf(" svg:stroke-width=\"%s\" ", s.strokeWidth.c_str());

e45563a3c46261d8c32014f8e516857ba01bd7b7miklosh outs.printf(" svg:stroke-color=\"%s\" ", s.strokeColor.c_str());

e45563a3c46261d8c32014f8e516857ba01bd7b7miklosh outs.printf(" svg:stroke-opacity=\"%s\" ", s.strokeOpacity.c_str());

e45563a3c46261d8c32014f8e516857ba01bd7b7miklosh }

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf("/>\n");

7ec85862d9730e449ed5c2a86201bc9ca1daa0aamiklosh outs.printf("</style:style>\n");

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh }

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh

7ec85862d9730e449ed5c2a86201bc9ca1daa0aamiklosh //## Dump our gradient table

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh int gradientCount = 0;

522aa9b8f493ba0c8e8b0bb536a563c96f5430a8miklosh outs.printf("\n");

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf("<!-- ####### Gradients from Inkscape document ####### -->\n");

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh std::vector<GradientInfo>::iterator giter;

e454b92b3d16b0909892cddef064b745898c924dJon A. Cruz for (giter = gradientTable.begin() ; giter != gradientTable.end() ; giter++)

e454b92b3d16b0909892cddef064b745898c924dJon A. Cruz {

e9b6af083e34e2397a8ddbe9781920733d09d151Ted Gould GradientInfo gi(*giter);

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh if (gi.style == "linear")

7ec85862d9730e449ed5c2a86201bc9ca1daa0aamiklosh {

84d6d1f7365e49f2936df9df890ce49d2c000ce2Kris /*

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh if (gi.stops.size() < 2)

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh {

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh g_warning("Need at least 2 tops for a linear gradient");

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh continue;

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh }

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf("<svg:linearGradient ");

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf("id=\"%#s_g\" ", gi.name.c_str());

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf("draw:name=\"%#s_g\"\n", gi.name.c_str());

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf(" draw:display-name=\"imported linear %d\"\n",

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh gradientCount);

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf(" svg:x1=\"%05.3fcm\" svg:y1=\"%05.3fcm\"\n",

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh gi.x1, gi.y1);

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf(" svg:x2=\"%05.3fcm\" svg:y2=\"%05.3fcm\"\n",

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh gi.x2, gi.y2);

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf(" svg:gradientUnits=\"objectBoundingBox\">\n");

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf(" <svg:stop\n");

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf(" svg:stop-color=\"#%06lx\"\n",

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh gi.stops[0].rgb);

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf(" svg:stop-opacity=\"%f%%\"\n",

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh gi.stops[0].opacity * 100.0);

405079a6cf76f81a583a9b8556c1ca49840abc6eKris outs.printf(" svg:offset=\"0\"/>\n");

17d87f5698f5c2958d38c6a6207c7b322a7adaf9johanengelen outs.printf(" <svg:stop\n");

797add969622b3786436365fdd275aff9a9a6373miklosh outs.printf(" svg:stop-color=\"#%06lx\"\n",

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh gi.stops[1].rgb);

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh outs.printf(" svg:stop-opacity=\"%f%%\"\n",

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh gi.stops[1].opacity * 100.0);

e454b92b3d16b0909892cddef064b745898c924dJon A. Cruz outs.printf(" svg:offset=\"1\"/>\n");

e454b92b3d16b0909892cddef064b745898c924dJon A. Cruz outs.printf("</svg:linearGradient>\n");

17d87f5698f5c2958d38c6a6207c7b322a7adaf9johanengelen }

e45563a3c46261d8c32014f8e516857ba01bd7b7miklosh else if (gi.style == "radial")

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh {

3711b3e25395437ee0a09dbbb2a76d999c4ef322miklosh /*