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
80
81
82
83
84
85
86
87
88
|
/*1:*/
#line 6 "./t_polynomial.cweb"
#include "t_polynomial.h"
#include "kron_prod.h"
/*2:*/
#line 20 "./t_polynomial.cweb"
const URSingleTensor&PowerProvider::getNext(const URSingleTensor*dummy)
{
if(ut){
URSingleTensor*ut_new= new URSingleTensor(nv,ut->dimen()+1);
KronProd::kronMult(ConstVector(origv),ConstVector(ut->getData()),ut_new->getData());
delete ut;
ut= ut_new;
}else{
ut= new URSingleTensor(nv,1);
ut->getData()= origv;
}
return*ut;
}
/*:2*/
#line 10 "./t_polynomial.cweb"
;
/*3:*/
#line 38 "./t_polynomial.cweb"
const FRSingleTensor&PowerProvider::getNext(const FRSingleTensor*dummy)
{
getNext(ut);
if(ft)
delete ft;
ft= new FRSingleTensor(*ut);
return*ft;
}
/*:3*/
#line 11 "./t_polynomial.cweb"
;
/*4:*/
#line 49 "./t_polynomial.cweb"
PowerProvider::~PowerProvider()
{
if(ut)
delete ut;
if(ft)
delete ft;
}
/*:4*/
#line 12 "./t_polynomial.cweb"
;
/*5:*/
#line 59 "./t_polynomial.cweb"
UTensorPolynomial::UTensorPolynomial(const FTensorPolynomial&fp)
:TensorPolynomial<UFSTensor,UGSTensor,URSingleTensor> (fp.nrows(),fp.nvars())
{
for(FTensorPolynomial::const_iterator it= fp.begin();
it!=fp.end();++it){
insert(new UFSTensor(*((*it).second)));
}
}
/*:5*/
#line 13 "./t_polynomial.cweb"
;
/*6:*/
#line 70 "./t_polynomial.cweb"
FTensorPolynomial::FTensorPolynomial(const UTensorPolynomial&up)
:TensorPolynomial<FFSTensor,FGSTensor,FRSingleTensor> (up.nrows(),up.nvars())
{
for(UTensorPolynomial::const_iterator it= up.begin();
it!=up.end();++it){
insert(new FFSTensor(*((*it).second)));
}
}
/*:6*/
#line 14 "./t_polynomial.cweb"
;
/*:1*/
|
