From cca29d5604f37639a58fba5108ecc07732aa7316 Mon Sep 17 00:00:00 2001 From: Ralf Hartmann Date: Mon, 18 Dec 2000 12:44:25 +0000 Subject: [PATCH] Remove old implementation of shape_values_and_grads. git-svn-id: https://svn.dealii.org/trunk@3554 0785d39b-7218-0410-832d-ea1e28bc413d --- .../base/source/tensor_product_polynomials.cc | 171 +----------------- 1 file changed, 9 insertions(+), 162 deletions(-) diff --git a/deal.II/base/source/tensor_product_polynomials.cc b/deal.II/base/source/tensor_product_polynomials.cc index 46803d9325..024eb3aeeb 100644 --- a/deal.II/base/source/tensor_product_polynomials.cc +++ b/deal.II/base/source/tensor_product_polynomials.cc @@ -15,7 +15,6 @@ #include #include -#include template @@ -24,157 +23,7 @@ TensorProductPolynomials::TensorProductPolynomials( polynomials(pols) {} - // first version of this function -// template -// void TensorProductPolynomials::shape_values_and_grads( -// const Point &p, -// vector &values, -// vector > &grads, -// vector > &grad_grads) const -// { -// unsigned int n_pols=polynomials->size(); -// unsigned int n_pols2=n_pols*n_pols; -// unsigned int n_tensor_pols=pow(n_pols, dim); -// Assert(values.size()==n_tensor_pols || values.size()==0, -// ExcDimensionMismatch2(values.size(), n_tensor_pols, 0)); -// Assert(grads.size()==n_tensor_pols|| grads.size()==0, -// ExcDimensionMismatch2(grads.size(), n_tensor_pols, 0)); -// Assert(grad_grads.size()==n_tensor_pols|| grad_grads.size()==0, -// ExcDimensionMismatch2(grad_grads.size(), n_tensor_pols, 0)); -// unsigned int v_size=0; -// bool update_values=false, update_grads=false, update_grad_grads=false; -// if (values.size()==n_tensor_pols) -// { -// update_values=true; -// v_size=1; -// } -// if (grads.size()==n_tensor_pols) -// { -// update_grads=true; -// v_size=2; -// } -// if (grad_grads.size()==n_tensor_pols) -// { -// update_grad_grads=true; -// v_size=3; -// } - -// vector > > v( -// dim, vector > (n_pols, vector (v_size))); - -// for (unsigned int d=0; d > &v_d=v[d]; -// Assert(v_d.size()==n_pols, ExcInternalError()); -// for (unsigned int i=0; ioperator[](i)->value(p(d), v_d[i]); -// } - -// if (update_values) -// { -// // values[i]=p_j(x) -// for (unsigned i=0; i1) -// for (unsigned i=0; i2) -// for (unsigned int i=0; i1) -// for (unsigned i=0; i2) -// for (unsigned i=0; i -number pow(number x, unsigned int y) -{ - number value=1; - for (unsigned int i=0; i void TensorProductPolynomials::shape_values_and_grads( const Point &p, @@ -183,7 +32,12 @@ void TensorProductPolynomials::shape_values_and_grads( vector > &grad_grads) const { unsigned int n_pols=polynomials.size(); - unsigned int n_tensor_pols=pow(n_pols, dim); + vector n_pols_to(dim+1); + n_pols_to[0]=1; + for (unsigned int i=0; i::shape_values_and_grads( polynomials[i]->value(p(d), v_d[i]); } - vector n_pols_to(dim+1); - n_pols_to[0]=1; - for (unsigned int i=0; i