/* $Id$ */
/* Version: $Name$ */
/* */
-/* Copyright (C) 2001, 2002, 2003, 2004, 2006, 2007 by the deal.II authors */
+/* Copyright (C) 2001, 2002, 2003, 2004, 2006, 2007, 2008 by the deal.II authors */
/* */
/* This file is subject to QPL and may not be distributed */
/* without copyright and license information. Please refer */
update_gradients | update_JxW_values);
const unsigned int dofs_per_cell = fe->dofs_per_cell;
- const unsigned int n_q_points = quadrature->n_quadrature_points;
+ const unsigned int n_q_points = quadrature->size();
FullMatrix<double> cell_matrix (dofs_per_cell, dofs_per_cell);
update_JxW_values);
const unsigned int dofs_per_cell = this->fe->dofs_per_cell;
- const unsigned int n_q_points = this->quadrature->n_quadrature_points;
+ const unsigned int n_q_points = this->quadrature->size();
Vector<double> cell_rhs (dofs_per_cell);
std::vector<double> rhs_values (n_q_points);
/* $Id$ */
/* Version: $Name$ */
/* */
-/* Copyright (C) 2002, 2003, 2004, 2005, 2006, 2007 by the deal.II authors */
+/* Copyright (C) 2002, 2003, 2004, 2005, 2006, 2007, 2008 by the deal.II authors */
/* */
/* This file is subject to QPL and may not be distributed */
/* without copyright and license information. Please refer */
vertex_quadrature,
update_gradients | update_quadrature_points);
std::vector<Tensor<1,dim> >
- solution_gradients (vertex_quadrature.n_quadrature_points);
+ solution_gradients (vertex_quadrature.size());
// ...and next loop over all cells
// and their vertices, and count
update_gradients | update_JxW_values);
const unsigned int dofs_per_cell = fe->dofs_per_cell;
- const unsigned int n_q_points = quadrature->n_quadrature_points;
+ const unsigned int n_q_points = quadrature->size();
FullMatrix<double> cell_matrix (dofs_per_cell, dofs_per_cell);
update_JxW_values);
const unsigned int dofs_per_cell = this->fe->dofs_per_cell;
- const unsigned int n_q_points = this->quadrature->n_quadrature_points;
+ const unsigned int n_q_points = this->quadrature->size();
Vector<double> cell_rhs (dofs_per_cell);
std::vector<double> rhs_values (n_q_points);
right_hand_side (&right_hand_side)
{
const unsigned int n_q_points
- = quadrature.n_quadrature_points;
+ = quadrature.size();
cell_residual.resize(n_q_points);
rhs_values.resize(n_q_points);
update_gradients)
{
const unsigned int n_face_q_points
- = face_quadrature.n_quadrature_points;
+ = face_quadrature.size();
jump_residual.resize(n_face_q_points);
dual_weights.resize(n_face_q_points);
/* $Id$ */
/* Version: $Name$ */
/* */
-/* Copyright (C) 2002, 2003, 2004, 2006, 2007 by the deal.II authors */
+/* Copyright (C) 2002, 2003, 2004, 2006, 2007, 2008 by the deal.II authors */
/* */
/* This file is subject to QPL and may not be distributed */
/* without copyright and license information. Please refer */
// hold the indices of the local degrees of
// freedom on each cell:
const unsigned int dofs_per_cell = fe.dofs_per_cell;
- const unsigned int n_q_points = quadrature_formula.n_quadrature_points;
+ const unsigned int n_q_points = quadrature_formula.size();
FullMatrix<double> cell_matrix (dofs_per_cell, dofs_per_cell);
Vector<double> cell_rhs (dofs_per_cell);
// at quadrature points. Here, we also need
// second derivatives, which is simple,
// however:
- std::vector<double> local_values (quadrature.n_quadrature_points);
- std::vector<Tensor<1,dim> > local_gradients (quadrature.n_quadrature_points);
- std::vector<Tensor<2,dim> > local_2nd_derivs (quadrature.n_quadrature_points);
+ std::vector<double> local_values (quadrature.size());
+ std::vector<Tensor<1,dim> > local_gradients (quadrature.size());
+ std::vector<Tensor<2,dim> > local_2nd_derivs (quadrature.size());
// With all this, we can start the loop
// over all cells. Since we need to write
// since we extract the gradient as a
// scalar value.
double cell_residual_norm = 0;
- for (unsigned int q=0; q<quadrature.n_quadrature_points; ++q)
+ for (unsigned int q=0; q<quadrature.size(); ++q)
{
const double x = fe_values.quadrature_point(q)[0];
const double u = local_values[q];
update_values | update_gradients |
update_quadrature_points | update_JxW_values);
- const unsigned int n_q_points = quadrature_formula.n_quadrature_points;
+ const unsigned int n_q_points = quadrature_formula.size();
// Then, just as when we integrated the
// linear system, we need two variables
/* $Id$ */
/* Version: $Name$ */
/* */
-/* Copyright (C) 2003, 2004, 2005, 2006, 2007 by the deal.II authors */
+/* Copyright (C) 2003, 2004, 2005, 2006, 2007, 2008 by the deal.II authors */
/* */
/* This file is subject to QPL and may not be distributed */
/* without copyright and license information. Please refer */
update_quadrature_points | update_JxW_values);
const unsigned int dofs_per_cell = fe.dofs_per_cell;
- const unsigned int n_q_points = quadrature_formula.n_quadrature_points;
+ const unsigned int n_q_points = quadrature_formula.size();
FullMatrix<double> cell_matrix (dofs_per_cell, dofs_per_cell);
Vector<double> cell_rhs (dofs_per_cell);
update_quadrature_points | update_JxW_values);
const unsigned int dofs_per_cell = fe.dofs_per_cell;
- const unsigned int n_q_points = quadrature_formula.n_quadrature_points;
+ const unsigned int n_q_points = quadrature_formula.size();
FullMatrix<double> cell_matrix (dofs_per_cell, dofs_per_cell);
/* $Id$ */
/* Version: $Name$ */
/* */
-/* Copyright (C) 2000, 2004, 2005, 2006, 2007 by the deal.II authors */
+/* Copyright (C) 2000, 2004, 2005, 2006, 2007, 2008 by the deal.II authors */
/* */
/* This file is subject to QPL and may not be distributed */
/* without copyright and license information. Please refer */
update_quadrature_points | update_JxW_values);
const unsigned int dofs_per_cell = fe.dofs_per_cell;
- const unsigned int n_q_points = quadrature_formula.n_quadrature_points;
+ const unsigned int n_q_points = quadrature_formula.size();
FullMatrix<double> cell_matrix (dofs_per_cell, dofs_per_cell);
Vector<double> cell_rhs (dofs_per_cell);
/* $Id$ */
/* Version: $Name$ */
/* */
-/* Copyright (C) 2000, 2004, 2005, 2006, 2007 by the deal.II authors */
+/* Copyright (C) 2000, 2004, 2005, 2006, 2007, 2008 by the deal.II authors */
/* */
/* This file is subject to QPL and may not be distributed */
/* without copyright and license information. Please refer */
update_quadrature_points | update_JxW_values);
const unsigned int dofs_per_cell = fe.dofs_per_cell;
- const unsigned int n_q_points = quadrature_formula.n_quadrature_points;
+ const unsigned int n_q_points = quadrature_formula.size();
FullMatrix<double> cell_matrix (dofs_per_cell, dofs_per_cell);
Vector<double> cell_rhs (dofs_per_cell);
// points...
SymmetricTensor<2,dim> accumulated_stress;
for (unsigned int q=0;
- q<quadrature_formula.n_quadrature_points;
+ q<quadrature_formula.size();
++q)
accumulated_stress +=
reinterpret_cast<PointHistory<dim>*>(cell->user_pointer())[q]
// average to their destination:
norm_of_stress(index)
= (accumulated_stress /
- quadrature_formula.n_quadrature_points).norm();
+ quadrature_formula.size()).norm();
}
// And on the cells that we are not
// interested in, set the respective
tmp.swap (quadrature_point_history);
}
quadrature_point_history.resize (our_cells *
- quadrature_formula.n_quadrature_points);
+ quadrature_formula.size());
// Finally loop over all cells again and
// set the user pointers from the cells
if (cell->subdomain_id() == this_mpi_process)
{
cell->set_user_pointer (&quadrature_point_history[history_index]);
- history_index += quadrature_formula.n_quadrature_points;
+ history_index += quadrature_formula.size();
}
// At the end, for good measure make sure
FEValues<dim> fe_values (fe, quadrature_formula,
update_values | update_gradients);
std::vector<std::vector<Tensor<1,dim> > >
- displacement_increment_grads (quadrature_formula.n_quadrature_points,
+ displacement_increment_grads (quadrature_formula.size(),
std::vector<Tensor<1,dim> >(dim));
// Then loop over all cells and do the
// Then loop over the quadrature
// points of this cell:
- for (unsigned int q=0; q<quadrature_formula.n_quadrature_points; ++q)
+ for (unsigned int q=0; q<quadrature_formula.size(); ++q)
{
// On each quadrature point,
// compute the strain increment
/* $Id$ */
/* Version: $Name$ */
/* */
-/* Copyright (C) 2005, 2006, 2007 by the deal.II authors */
+/* Copyright (C) 2005, 2006, 2007, 2008 by the deal.II authors */
/* */
/* This file is subject to QPL and may not be distributed */
/* without copyright and license information. Please refer */
update_quadrature_points | update_JxW_values);
const unsigned int dofs_per_cell = fe.dofs_per_cell;
- const unsigned int n_q_points = quadrature_formula.n_quadrature_points;
- const unsigned int n_face_q_points = face_quadrature_formula.n_quadrature_points;
+ const unsigned int n_q_points = quadrature_formula.size();
+ const unsigned int n_face_q_points = face_quadrature_formula.size();
FullMatrix<double> local_matrix (dofs_per_cell, dofs_per_cell);
Vector<double> local_rhs (dofs_per_cell);
/* $Id$ */
/* Version: $Name$ */
/* */
-/* Copyright (C) 2006, 2007 by the deal.II authors */
+/* Copyright (C) 2006, 2007, 2008 by the deal.II authors */
/* */
/* This file is subject to QPL and may not be distributed */
/* without copyright and license information. Please refer */
const unsigned int dofs_per_cell = fe.dofs_per_cell;
- const unsigned int n_q_points = quadrature_formula.n_quadrature_points;
- const unsigned int n_face_q_points = face_quadrature_formula.n_quadrature_points;
+ const unsigned int n_q_points = quadrature_formula.size();
+ const unsigned int n_face_q_points = face_quadrature_formula.size();
FullMatrix<double> local_matrix (dofs_per_cell, dofs_per_cell);
Vector<double> local_rhs (dofs_per_cell);
update_values);
const unsigned int dofs_per_cell = fe.dofs_per_cell;
- const unsigned int n_q_points = quadrature_formula.n_quadrature_points;
- const unsigned int n_face_q_points = face_quadrature_formula.n_quadrature_points;
+ const unsigned int n_q_points = quadrature_formula.size();
+ const unsigned int n_face_q_points = face_quadrature_formula.size();
Vector<double> local_rhs (dofs_per_cell);
{
QGauss<dim> quadrature_formula(degree+2);
const unsigned int n_q_points
- = quadrature_formula.n_quadrature_points;
+ = quadrature_formula.size();
FEValues<dim> fe_values (fe, quadrature_formula,
update_values);
/* $Id$ */
/* Version: $Name$ */
/* */
-/* Copyright (C) 2007 by the deal.II authors */
+/* Copyright (C) 2007, 2008 by the deal.II authors */
/* */
/* This file is subject to QPL and may not be distributed */
/* without copyright and license information. Please refer */
const unsigned int dofs_per_cell = fe.dofs_per_cell;
- const unsigned int n_q_points = quadrature_formula.n_quadrature_points;
- const unsigned int n_face_q_points = face_quadrature_formula.n_quadrature_points;
+ const unsigned int n_q_points = quadrature_formula.size();
+ const unsigned int n_face_q_points = face_quadrature_formula.size();
FullMatrix<double> local_matrix (dofs_per_cell, dofs_per_cell);
Vector<double> local_rhs (dofs_per_cell);
update_values);
const unsigned int dofs_per_cell = fe.dofs_per_cell;
- const unsigned int n_q_points = quadrature_formula.n_quadrature_points;
- const unsigned int n_face_q_points = face_quadrature_formula.n_quadrature_points;
+ const unsigned int n_q_points = quadrature_formula.size();
+ const unsigned int n_face_q_points = face_quadrature_formula.size();
Vector<double> local_rhs (dofs_per_cell);
{
QGauss<dim> quadrature_formula(degree+2);
const unsigned int n_q_points
- = quadrature_formula.n_quadrature_points;
+ = quadrature_formula.size();
FEValues<dim> fe_values (fe, quadrature_formula,
update_values);
/* $Id$ */
/* Version: $Name: $ */
/* */
-/* Copyright (C) 2006, 2007 by the deal.II authors */
+/* Copyright (C) 2006, 2007, 2008 by the deal.II authors */
/* Author: Xing Jin, Wolfgang Bangerth, Texas A&M University, 2006 */
/* */
/* This file is subject to QPL and may not be distributed */
update_values | update_JxW_values);
const unsigned int dofs_per_cell = fe.dofs_per_cell;
- const unsigned int n_q_points = quadrature_formula.n_quadrature_points;
+ const unsigned int n_q_points = quadrature_formula.size();
FullMatrix<double> cell_matrix (dofs_per_cell, dofs_per_cell);
/* $Id$ */
-/* Copyright (C) 2006, 2007 by the deal.II authors */
+/* Copyright (C) 2006, 2007, 2008 by the deal.II authors */
/* Author: Ivan Christov, Wolfgang Bangerth, Texas A&M University, 2006 */
/* */
/* This file is subject to QPL and may not be distributed */
update_quadrature_points);
const unsigned int dofs_per_cell = fe.dofs_per_cell;
- const unsigned int n_q_points = quadrature_formula.n_quadrature_points;
+ const unsigned int n_q_points = quadrature_formula.size();
Vector<double> local_nl_term (dofs_per_cell);
std::vector<unsigned int> local_dof_indices (dofs_per_cell);
update_values | update_JxW_values | update_quadrature_points);
const unsigned int dofs_per_cell = fe.dofs_per_cell;
- const unsigned int n_q_points = quadrature_formula.n_quadrature_points;
+ const unsigned int n_q_points = quadrature_formula.size();
FullMatrix<double> local_nl_matrix (dofs_per_cell, dofs_per_cell);
std::vector<unsigned int> local_dof_indices (dofs_per_cell);
/* $Id$ */
/* Version: $Name$ */
/* */
-/* Copyright (C) 2006, 2007 by the deal.II authors */
+/* Copyright (C) 2006, 2007, 2008 by the deal.II authors */
/* */
/* This file is subject to QPL and may not be distributed */
/* without copyright and license information. Please refer */
for (unsigned int j=0; j<fe_collection[fe].dofs_per_cell; ++j)
{
std::complex<double> sum = 0;
- for (unsigned int q=0; q<quadrature.n_quadrature_points; ++q)
+ for (unsigned int q=0; q<quadrature.size(); ++q)
{
const Point<dim> x_q = quadrature.point(q);
sum += std::exp(std::complex<double>(0,1) *
/* $Id$ */
/* Version: $Name: $ */
/* */
-/* Copyright (C) 2006, 2007 by the deal.II authors and Yaqi Wang */
+/* Copyright (C) 2006, 2007, 2008 by the deal.II authors and Yaqi Wang */
/* */
/* This file is subject to QPL and may not be distributed */
/* without copyright and license information. Please refer */
update_JxW_values);
const unsigned int dofs_per_cell = fe.dofs_per_cell;
- const unsigned int n_q_points = quadrature_formula.n_quadrature_points;
+ const unsigned int n_q_points = quadrature_formula.size();
FullMatrix<double> cell_matrix (dofs_per_cell, dofs_per_cell);
Vector<double> cell_rhs (dofs_per_cell);
const QGauss<dim> quadrature_formula (fe.degree + 1);
const unsigned int dofs_per_cell = fe.dofs_per_cell;
- const unsigned int n_q_points = quadrature_formula.n_quadrature_points;
+ const unsigned int n_q_points = quadrature_formula.size();
FEValues<dim> fe_values (fe, quadrature_formula,
update_values | update_quadrature_points |
if (!cell_g->has_children() && !cell_g_prime->has_children())
{
const QGauss<dim> quadrature_formula (fe.degree+1);
- const unsigned int n_q_points = quadrature_formula.n_quadrature_points;
+ const unsigned int n_q_points = quadrature_formula.size();
FEValues<dim> fe_values (fe, quadrature_formula,
update_values | update_JxW_values);
double EnergyGroup<dim>::get_fission_source () const
{
const QGauss<dim> quadrature_formula (fe.degree + 1);
- const unsigned int n_q_points = quadrature_formula.n_quadrature_points;
+ const unsigned int n_q_points = quadrature_formula.size();
FEValues<dim> fe_values (fe, quadrature_formula,
update_values | update_JxW_values);
/* Author: Moritz Allmaras, Texas A&M University, 2007 */
/* */
-/* Copyright (C) 2007 by the deal.II authors and M. Allmaras */
+/* Copyright (C) 2007, 2008 by the deal.II authors and M. Allmaras */
/* */
/* This file is subject to QPL and may not be distributed */
/* without copyright and license information. Please refer */
QGauss<dim> quadrature_formula(2);
QGauss<dim-1> face_quadrature_formula(2);
- const unsigned int n_q_points = quadrature_formula.n_quadrature_points,
- n_face_q_points = face_quadrature_formula.n_quadrature_points,
+ const unsigned int n_q_points = quadrature_formula.size(),
+ n_face_q_points = face_quadrature_formula.size(),
dofs_per_cell = fe.dofs_per_cell;
// The FEValues objects will
/* $Id$ */
/* Version: $Name$ */
/* */
-/* Copyright (C) 1999, 2000, 2001, 2002, 2003, 2004, 2006 by the deal.II authors */
+/* Copyright (C) 1999, 2000, 2001, 2002, 2003, 2004, 2006, 2007, 2008 by the deal.II authors */
/* */
/* This file is subject to QPL and may not be distributed */
/* without copyright and license information. Please refer */
// more or less by convention the standard
// names for these purposes:
const unsigned int dofs_per_cell = fe.dofs_per_cell;
- const unsigned int n_q_points = quadrature_formula.n_quadrature_points;
+ const unsigned int n_q_points = quadrature_formula.size();
// Now, we said that we wanted to assemble
// the global matrix and vector
/* $Id$ */
/* Version: $Name$ */
/* */
-/* Copyright (C) 1999, 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2007 by the deal.II authors */
+/* Copyright (C) 1999, 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2007, 2008 by the deal.II authors */
/* */
/* This file is subject to QPL and may not be distributed */
/* without copyright and license information. Please refer */
// you and you don't have to care about the
// dimension dependent parts:
const unsigned int dofs_per_cell = fe.dofs_per_cell;
- const unsigned int n_q_points = quadrature_formula.n_quadrature_points;
+ const unsigned int n_q_points = quadrature_formula.size();
FullMatrix<double> cell_matrix (dofs_per_cell, dofs_per_cell);
Vector<double> cell_rhs (dofs_per_cell);
/* $Id$ */
/* Version: $Name$ */
/* */
-/* Copyright (C) 1999, 2000, 2001, 2002, 2003, 2004, 2006, 2007 by the deal.II authors */
+/* Copyright (C) 1999, 2000, 2001, 2002, 2003, 2004, 2006, 2007, 2008 by the deal.II authors */
/* */
/* This file is subject to QPL and may not be distributed */
/* without copyright and license information. Please refer */
update_quadrature_points | update_JxW_values);
const unsigned int dofs_per_cell = fe.dofs_per_cell;
- const unsigned int n_q_points = quadrature_formula.n_quadrature_points;
+ const unsigned int n_q_points = quadrature_formula.size();
FullMatrix<double> cell_matrix (dofs_per_cell, dofs_per_cell);
Vector<double> cell_rhs (dofs_per_cell);
/* $Id$ */
/* Version: $Name$ */
/* */
-/* Copyright (C) 2000, 2001, 2002, 2003, 2004, 2006, 2007 by the deal.II authors */
+/* Copyright (C) 2000, 2001, 2002, 2003, 2004, 2006, 2007, 2008 by the deal.II authors */
/* */
/* This file is subject to QPL and may not be distributed */
/* without copyright and license information. Please refer */
update_quadrature_points | update_JxW_values);
const unsigned int dofs_per_cell = fe.dofs_per_cell;
- const unsigned int n_q_points = quadrature_formula.n_quadrature_points;
+ const unsigned int n_q_points = quadrature_formula.size();
FullMatrix<double> cell_matrix (dofs_per_cell, dofs_per_cell);
Vector<double> cell_rhs (dofs_per_cell);
/* $Id$ */
/* Version: $Name$ */
/* */
-/* Copyright (C) 2000, 2001, 2002, 2003, 2004, 2006, 2007 by the deal.II authors */
+/* Copyright (C) 2000, 2001, 2002, 2003, 2004, 2006, 2007, 2008 by the deal.II authors */
/* */
/* This file is subject to QPL and may not be distributed */
/* without copyright and license information. Please refer */
QGauss<dim> quadrature_formula(3);
QGauss<dim-1> face_quadrature_formula(3);
- const unsigned int n_q_points = quadrature_formula.n_quadrature_points;
- const unsigned int n_face_q_points = face_quadrature_formula.n_quadrature_points;
+ const unsigned int n_q_points = quadrature_formula.size();
+ const unsigned int n_face_q_points = face_quadrature_formula.size();
const unsigned int dofs_per_cell = fe->dofs_per_cell;
/* $Id$ */
/* Version: $Name$ */
/* */
-/* Copyright (C) 2000, 2001, 2002, 2003, 2004, 2006, 2007 by the deal.II authors */
+/* Copyright (C) 2000, 2001, 2002, 2003, 2004, 2006, 2007, 2008 by the deal.II authors */
/* */
/* This file is subject to QPL and may not be distributed */
/* without copyright and license information. Please refer */
update_quadrature_points | update_JxW_values);
const unsigned int dofs_per_cell = fe.dofs_per_cell;
- const unsigned int n_q_points = quadrature_formula.n_quadrature_points;
+ const unsigned int n_q_points = quadrature_formula.size();
FullMatrix<double> cell_matrix (dofs_per_cell, dofs_per_cell);
Vector<double> cell_rhs (dofs_per_cell);
/* $Id$ */
/* Version: $Name$ */
/* */
-/* Copyright (C) 2000, 2001, 2002, 2003, 2004, 2006, 2007 by the deal.II authors */
+/* Copyright (C) 2000, 2001, 2002, 2003, 2004, 2006, 2007, 2008 by the deal.II authors */
/* */
/* This file is subject to QPL and may not be distributed */
/* without copyright and license information. Please refer */
// abbreviations to avoid
// unnecessarily long lines:
const unsigned int dofs_per_cell = fe.dofs_per_cell;
- const unsigned int n_q_points = quadrature_formula.n_quadrature_points;
- const unsigned int n_face_q_points = face_quadrature_formula.n_quadrature_points;
+ const unsigned int n_q_points = quadrature_formula.size();
+ const unsigned int n_face_q_points = face_quadrature_formula.size();
// We declare cell matrix and cell
// right hand side...