// $Id$
// Version: $Name$
//
-// Copyright (C) 2001, 2002, 2003, 2004, 2005, 2006, 2007, 2008, 2009 by the deal.II authors
+// Copyright (C) 2001, 2002, 2003, 2004, 2005, 2006, 2007, 2008, 2009, 2010 by the deal.II authors
//
// This file is subject to QPL and may not be distributed
// without copyright and license information. Please refer
// sure we go at most
// step_length into this
// direction
- const Point<spacedim> old_point = object_mid_point;
object_mid_point -= std::min(2 * current_value / (gradient*gradient),
step_length / gradient.norm()) *
gradient;
quad = triangulation.begin_quad(),
endq = triangulation.end_quad();
typename Triangulation<dim,spacedim>::raw_line_iterator
- next_unused_line = triangulation.begin_raw_line (),
- endl = triangulation.end_line();
+ next_unused_line = triangulation.begin_raw_line ();
typename Triangulation<dim,spacedim>::raw_quad_iterator
next_unused_quad = triangulation.begin_raw_quad ();
typename Triangulation<dim,spacedim>::active_hex_iterator
hex = triangulation.begin_active_hex(level),
endh = triangulation.begin_active_hex(level+1);
- typename Triangulation<dim,spacedim>::raw_line_iterator
- next_unused_line = triangulation.begin_raw_line ();
- typename Triangulation<dim,spacedim>::raw_quad_iterator
- next_unused_quad = triangulation.begin_raw_quad ();
typename Triangulation<dim,spacedim>::raw_hex_iterator
next_unused_hex = triangulation.begin_raw_hex (level+1);
// $Id$
// Version: $Name$
//
-// Copyright (C) 1998, 1999, 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2008, 2009 by the deal.II authors
+// Copyright (C) 1998, 1999, 2000, 2001, 2002, 2003, 2004, 2005, 2006, 2008, 2009, 2010 by the deal.II authors
//
// This file is subject to QPL and may not be distributed
// without copyright and license information. Please refer
template <int dim, int spacedim>
Point<spacedim>
Boundary<dim, spacedim>::
-get_new_point_on_quad (const typename Triangulation<dim, spacedim>::quad_iterator &) const
+get_new_point_on_quad (const typename Triangulation<dim, spacedim>::quad_iterator &) const
{
Assert (false, ExcPureFunctionCalled());
return Point<spacedim>();
case 3:
return get_new_point_on_quad (face);
}
-
+
return Point<spacedim>();
}
template <int dim, int spacedim>
Point<spacedim>
StraightBoundary<dim, spacedim>::
-get_new_point_on_line (const typename Triangulation<dim, spacedim>::line_iterator &line) const
+get_new_point_on_line (const typename Triangulation<dim, spacedim>::line_iterator &line) const
{
return (line->vertex(0) + line->vertex(1)) / 2;
}
template <int dim, int spacedim>
Point<spacedim>
StraightBoundary<dim, spacedim>::
-get_new_point_on_quad (const typename Triangulation<dim, spacedim>::quad_iterator &) const
+get_new_point_on_quad (const typename Triangulation<dim, spacedim>::quad_iterator &) const
{
Assert (false, ExcImpossibleInDim(dim));
return Point<spacedim>();
template <>
Point<3>
StraightBoundary<2,3>::
-get_new_point_on_quad (const Triangulation<2,3>::quad_iterator &quad) const
+get_new_point_on_quad (const Triangulation<2,3>::quad_iterator &quad) const
{
return compute_new_point_on_quad<2> (quad);
}
template <>
Point<3>
StraightBoundary<3>::
-get_new_point_on_quad (const Triangulation<3>::quad_iterator &quad) const
+get_new_point_on_quad (const Triangulation<3>::quad_iterator &quad) const
{
return compute_new_point_on_quad<3> (quad);
}
{
const unsigned int n=points.size();
Assert(n>0, ExcInternalError());
-
+
const double dx=1./(n+1);
double x=dx;
const Point<spacedim> vertices[2] = { line->vertex(0),
line->vertex(1) };
-
+
for (unsigned int i=0; i<n; ++i, x+=dx)
points[i] = (1-x)*vertices[0] + x*vertices[1];
}
std::vector<Point<3> > &points) const
{
const unsigned int spacedim = 3;
-
+
const unsigned int n=points.size(),
m=static_cast<unsigned int>(std::sqrt(static_cast<double>(n)));
// is n a square number
x_k += object->vertex(i) *
GeometryInfo<dim>::d_linear_shape_function (xi, i);
- do
+ do
{
Tensor<1,dim> F_k;
for (unsigned int i=0; i<GeometryInfo<dim>::vertices_per_cell; ++i)
const Point<dim> delta_xi = - invert(H_k) * F_k;
xi += delta_xi;
- Point<spacedim> old_xk = x_k;
x_k = Point<spacedim>();
for (unsigned int i=0; i<GeometryInfo<dim>::vertices_per_cell; ++i)
x_k += object->vertex(i) *
break;
}
while (true);
-
+
return x_k;
}
// $Id$
// Version: $Name$
//
-// Copyright (C) 2009 by Michael Rapson and the deal.II authors
+// Copyright (C) 2009, 2010 by Michael Rapson and the deal.II authors
//
// This file is subject to QPL and may not be distributed
// without copyright and license information. Please refer
to_gnuplot << "# Data independent of mesh location\n";
// write column headings
- std::map <std::string, std::vector <std::vector <double> > >::iterator
- data_store_begin = data_store.begin ();
to_gnuplot << "# <Key> ";
for (unsigned int component = 0; component < n_indep; component++)
// write general data stored
for (unsigned int key = 0; key < dataset_key.size (); key++)
{
- std::map <std::string, std::vector <std::vector <double> > >::iterator
- data_store_begin = data_store.begin ();
to_gnuplot << dataset_key[key];
for (unsigned int component = 0; component < n_indep; component++)