79 |
return *__A; |
return *__A; |
80 |
} |
} |
81 |
|
|
82 |
|
//! Returns a pointer on a SecondOrderOperator that was "multiplied" by u |
83 |
|
ReferenceCounting<PDEOperator> operator * (ConstReferenceCounting<UserFunction> u) const |
84 |
|
{ |
85 |
|
ReferenceCounting<SecondOrderOperator::Matrix> A2 |
86 |
|
= new SecondOrderOperator::Matrix; |
87 |
|
|
88 |
|
for (size_t i=0; i<3; ++i) |
89 |
|
for (size_t j=0; j<3; ++j) { |
90 |
|
if ((*__A)(i,j) != 0) |
91 |
|
(*A2)(i,j) = new MulUserFunction((*__A)(i,j), u); |
92 |
|
} |
93 |
|
|
94 |
|
return (new SecondOrderOperator(A2)); |
95 |
|
} |
96 |
|
|
97 |
//! Returns a pointer on the opposed SecondOrderOperator operator. |
//! Returns a pointer on the opposed SecondOrderOperator operator. |
98 |
ReferenceCounting<PDEOperator> operator - () const |
ReferenceCounting<PDEOperator> operator - () const |
99 |
{ |
{ |
|
|
|
100 |
ReferenceCounting<SecondOrderOperator::Matrix> A2 |
ReferenceCounting<SecondOrderOperator::Matrix> A2 |
101 |
= new SecondOrderOperator::Matrix; |
= new SecondOrderOperator::Matrix; |
102 |
|
|