Discussion :

[sphere369]

Bonjour,

J'ai le message d'erreur suivant dans le code ci-dessous que je n'arrive pas a debugger.

1>TestPSOR.obj : error LNK2019: symbole externe non résolu "class Vector<double,int,class FullArray<double,class std::allocator<double> > > __cdecl operator+(class Vector<double,int,class FullArray<double,class std::allocator<double> > > const &,double const &)" (??H@YA?AV?$Vector@NHV?$FullArray@NV?$allocator@N@std@@@@@@ABV0@ABN@Z) référencé dans la fonction _main
1>TestPSOR.obj : error LNK2019: symbole externe non résolu "double __cdecl l2Norm<double,int>(class Vector<double,int,class FullArray<double,class std::allocator<double> > > const &,class Vector<double,int,class FullArray<double,class std::allocator<double> > > const &)" (??$l2Norm@NH@@YANABV?$Vector@NHV?$FullArray@NV?$allocator@N@std@@@@@@0@Z) référencé dans la fonction "public: virtual bool __thiscall MatrixIterativeSolver<double,int>::insideTolerance(void)const " (?insideTolerance@?$MatrixIterativeSolver@NH@@UBE_NXZ)
1>C:\MyProjects\Debug\test_1.exe : fatal error LNK1120: 2 externes non résolus

D'avance merci pour votr aide.

Sphere369.

TestPSOR.cpp

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
 
 
#include "matrixsolvermechanisms.cpp"
#include "vector.cpp"
#include <iostream>
#include "arraymechanisms.cpp"
#include "matrixmechanisms.cpp"
using namespace std;
 
 
int main()
{
	// Ingredient of problem, this is
	//
	//	Ax >= b, x >= c
	//	(x - c).(Ax - b) == 0 (inner product)
	//
 
	NumericMatrix<double, int> A(10, 10, 1, 1);		// The matrix to be inverted
	setDiagonal(A, 2.0);
	print(A);
 
	Vector<double, int> b(10, 1, 1.0);			// The right-hand side vector
	print(b);
 
	// Constraint stuff
	Vector<double, int> c(10, 1, 1.0);		
	Vector<double, int> v0 = c + 1.0;
	print(c); print(v0);
 
	PSORSolver<double, int> mySolver(A, b, c);
 
	/* virtual void initParameters(const V& tolerance, const V& SORFactor,
	const I& maxIter, NormType ntype,IterativeType itype);	// hook */
	mySolver.startVector(v0);	
	mySolver.initParameters(0.01, 1.0, 100, L2, Jacobi);
 
	mySolver.solve();
	Vector<double, int> res = mySolver.result();
	cout << "\nSolution for PSOR: \n"; print(res);
	cout << "\nNumber of iterations: " << mySolver.NumberIterations();
	if (mySolver.status() == false)
		cout << "\nStatus false";
	else
		cout << "\nStatus true";
 
	// Now test matrix iterative solvers using tridiagonal matrices
 
	int J = 10;
	Vector<double, int> a2(J-1,	2,	0.0);			// J-1 elements, start index == 2
	Vector<double, int> b2(J, 1,	2.0);
	Vector<double, int> c2(J-1,	1,	0.0);
	Vector<double, int> r2(J, 1, 1.0);				// Right-hand side
	Vector<double, int> constraintVector(J, 1, 1.0);			
 
	NumericMatrix<double, int> m2 = createMatrix(a2, b2, c2);		
	print(m2);
	print(r2);
 
	PSORSolver<double, int> mySolver2(m2, r2, constraintVector);
 
	v0 = r2 + 1.0;		// Linear offset
	mySolver2.startVector(v0);	
	mySolver2.initParameters(0.01, 1.0, 100, L2, Jacobi);
 
	mySolver2.solve();
	Vector<double, int> result3 = mySolver2.result();
 
	print(result3);
	cout << "\nNumber of iterations: " << mySolver2.NumberIterations();
	if (mySolver2.status() == false)
		cout << "\nStatus false";
	else
		cout << "\nStatus true";
 
 
	return 0;
 
}