Subversion Repositories gelsvn

Rev

Rev 357 | Go to most recent revision | Show entire file | Regard whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 357 Rev 595
Line -... Line 1...
- 
 
 1 
/* ----------------------------------------------------------------------- *
- 
 
 2 
 * This file is part of GEL, http://www.imm.dtu.dk/GEL
- 
 
 3 
 * Copyright (C) the authors and DTU Informatics
- 
 
 4 
 * For license and list of authors, see ../../doc/intro.pdf
- 
 
 5 
 * ----------------------------------------------------------------------- */
- 
 
 6 
 
1 
#include <iostream>
 7 
#include <iostream>
2 
 
 8 
 
3 
#include "QEM.h"
 9 
#include "QEM.h"
4 
#include "LinAlg/LapackFunc.h"
 10 
#include "LinAlg/LapackFunc.h"
5 
 
 11 
 
Line 34... Line 40...
34 
			}
 40 
        }
35 
 
 41
36 
			CMatrix Ap = V * Sp * U.Transposed();
 42 
        CMatrix Ap = V * Sp * U.Transposed();
37 
 
 43
38 
			CVector x(3);
 44 
        CVector x(3);
39 
			x = Ap * (b*0.5);
 45 
        x = Ap * CVector(b*0.5);
40
46
41 
			return diff-Vec3d(x[0],x[1],x[2]);
 47 
        return diff-Vec3d(x[0],x[1],x[2]);
42 
	}
 48 
	}
43 
 
 49
44 
}
 50 
}