Mercurial > hg > kwantix
diff diff_op.cpp @ 15:5144dd3c5468
Almost completed adding all Doxygen docstrings.
| author | Jordi Gutiérrez Hermoso <jordigh@gmail.com> |
|---|---|
| date | Fri, 08 Aug 2008 00:08:57 -0500 |
| parents | e6d2cd3b6e77 |
| children |
line wrap: on
line diff
--- a/diff_op.cpp +++ b/diff_op.cpp @@ -10,39 +10,45 @@ //************** Differential operator functions ********************* - diff_op::diff_op(){ - } - - double diff_op::operator()(const realfunc &f, const point &p) const{ + double diff_op::operator()(const realfunc &f, const point &p) const + { return at(f,p); } - double linear_diff_op2::operator()(const realfunc &f, const point &p) const{ + double linear_diff_op2::operator()(const realfunc &f, const point + &p) const + { return diff_op::operator()(f,p); } double bdry_diff_op::operator()(const realfunc &f, - const point &p) const{ + const point &p) const + { return at(f,p); } double bdry_diff_op::operator()(const realfunc &f, - const point &p, const vector &n) const{ + const point &p, const vector &n) + const + { return at(f,p,n); } double dirichlet_op::at(const realfunc &f, - const point &p) const{ + const point &p) const + { return f(p); } double dirichlet_op::at(const realfunc &f, const point &p, - const vector &n) const{ + const vector &n) const + { n.size(); return f(p); } double neumann_op::at(const realfunc &f, const point &p, - const vector &n) const{ + const vector &n) const + { size_t dim = n.size(); vector grad(dim); for(size_t i = 1; i <= dim; i++) @@ -50,28 +56,34 @@ return grad%n/norm(n); } - double Id_op::at(const realfunc &f, const point &p) const{ + double Id_op::at(const realfunc &f, const point &p) const + { return f(p); } - double del1::at(const realfunc &f, const point &p) const{ + double del1::at(const realfunc &f, const point &p) const + { return f.d(p,1); } - double del1::at(const realfunc &f, const point &p, size_t i) const{ + double del1::at(const realfunc &f, const point &p, size_t i) const + { return f.d(p,i); } - double del2::at(const realfunc &f, const point &p) const{ + double del2::at(const realfunc &f, const point &p) const + { return f.d2(p,1,1); } double del2::at(const realfunc &f, const point &p, - size_t i, size_t j) const{ + size_t i, size_t j) const + { return f.d2(p,i,j); } - double Laplacian::at(const realfunc &f, const point &p) const{ + double Laplacian::at(const realfunc &f, const point &p) const + { size_t n = p.size(); double result = 0; del2 D2;