phynexis-deps/igl
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

compress in tutorial, better header in compress, ifndef IGL_NO_OPENGL and IGL_NO_ANTTWEAKBAR guards, compress compiles without these, undo svd

Browse Source 
Former-commit-id: 7acb909737
This commit is contained in:
Alec Jacobson (jalec
2013-09-12 15:31:20 +02:00
parent 9e0f5c7a15
commit 79341b5bd6
58 changed files with 508 additions and 539 deletions
Download Patch File Download Diff File Expand all files Collapse all files
examples/svd/example.cpp
+97 -96
View File
@@ -1,96 +1,97 @@
#include <igl/svd.h>
#include <cstdlib>
#include <Accelerate/Accelerate.h>
#include <cstdio>
/* Auxiliary routines prototypes */
extern void print_matrix( char* desc, int m, int n, double* a, int lda );
/* Parameters */
void print3x3(const char * s, double * a)
{
printf("%s =\n",s);
for(int i = 0;i<3;i++)
{
for(int j = 0;j<3;j++)
{
printf("%g ",a[j*3+i]);
}
printf("\n");
}
printf("\n");
}
int main(int argc, char * argv[])
{
//// List of rest positions
//// (0,1)
//// / \
//// / \
//// / \
//// / \
//// (-1,0)-----(1,0)
////
//double rest[3][3] = {
// {-1,0,0},
// {1,0,0},
// {0,1,0}};
//// List of pose positions
////
//// (0,1)
//// | \
//// | \
//// | (1,0)
//// | /
//// | /
//// (0,-1)
//double pose[3][3] = {
// {0,1,0},
// {0,-1,0},
// {1,0,0}};
//// Compute covariance matrix C
//double C[3*3];
//// Initialize to zero
//for(int i = 0;i<3*3;i++)
//{
// C[i] = 0;
//}
//// Loop over vertices
//for(int i = 0;i<3;i++)
//{
// // Compute outer product rest[i] * pose[i]
// // Loop over coordinates
// for(int j = 0;j<3;j++)
// {
// // Loop over coordinates
// for(int k = 0;k<3;k++)
// {
// C[k*3+j] = rest[i][j] * pose[i][k];
// }
// }
//}
//print3x3("C",C);
//
//double C[3*3] = {8,3,4,1,5,9,6,7,2};
double C[3*3] = {5242.55,3364,-0,-8170.15,-5242.56,0,-0,-0,0};
double u[3*3],s[3],vt[3*3];
print3x3("C",C);
// Compute SVD of C
igl::svd3x3(C,u,s,vt);
print3x3("u",u);
print3x3("vt",vt);
// Compute R = u*vt
double R[3*3];
const double _3 = 3;
const double _1 = 1;
cblas_dgemm(CblasColMajor, CblasNoTrans,CblasNoTrans,3,3,3,1,u,3,vt,3,1,R,3);
print3x3("RT (transposed to be row-major)",R);
return 0;
}
int main(){}
//#include <igl/svd.h>
//#include <cstdlib>
//#include <Accelerate/Accelerate.h>
//#include <cstdio>
//
///* Auxiliary routines prototypes */
//extern void print_matrix( char* desc, int m, int n, double* a, int lda );
//
///* Parameters */
//
//void print3x3(const char * s, double * a)
//{
// printf("%s =\n",s);
// for(int i = 0;i<3;i++)
// {
// for(int j = 0;j<3;j++)
// {
// printf("%g ",a[j*3+i]);
// }
// printf("\n");
// }
// printf("\n");
//}
//
//int main(int argc, char * argv[])
//{
// //// List of rest positions
// //// (0,1)
// //// / \
// //// / \
// //// / \
// //// / \
// //// (-1,0)-----(1,0)
// ////
// //double rest[3][3] = {
// // {-1,0,0},
// // {1,0,0},
// // {0,1,0}};
// //// List of pose positions
// ////
// //// (0,1)
// //// | \
// //// | \
// //// | (1,0)
// //// | /
// //// | /
// //// (0,-1)
// //double pose[3][3] = {
// // {0,1,0},
// // {0,-1,0},
// // {1,0,0}};
// //// Compute covariance matrix C
// //double C[3*3];
// //// Initialize to zero
// //for(int i = 0;i<3*3;i++)
// //{
// // C[i] = 0;
// //}
//
// //// Loop over vertices
// //for(int i = 0;i<3;i++)
// //{
// // // Compute outer product rest[i] * pose[i]
// // // Loop over coordinates
// // for(int j = 0;j<3;j++)
// // {
// // // Loop over coordinates
// // for(int k = 0;k<3;k++)
// // {
// // C[k*3+j] = rest[i][j] * pose[i][k];
// // }
// // }
// //}
// //print3x3("C",C);
//
//
// //
// //double C[3*3] = {8,3,4,1,5,9,6,7,2};
// double C[3*3] = {5242.55,3364,-0,-8170.15,-5242.56,0,-0,-0,0};
// double u[3*3],s[3],vt[3*3];
// print3x3("C",C);
// // Compute SVD of C
// igl::svd3x3(C,u,s,vt);
// print3x3("u",u);
// print3x3("vt",vt);
//
// // Compute R = u*vt
// double R[3*3];
// const double _3 = 3;
// const double _1 = 1;
// cblas_dgemm(CblasColMajor, CblasNoTrans,CblasNoTrans,3,3,3,1,u,3,vt,3,1,R,3);
// print3x3("RT (transposed to be row-major)",R);
//
//
// return 0;
//}