/usr/include/deal.II/algorithms/operator.templates.h is in libdeal.ii-dev 8.4.2-2+b1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 | // ---------------------------------------------------------------------
//
// Copyright (C) 2006 - 2015 by the deal.II authors
//
// This file is part of the deal.II library.
//
// The deal.II library is free software; you can use it, redistribute
// it, and/or modify it under the terms of the GNU Lesser General
// Public License as published by the Free Software Foundation; either
// version 2.1 of the License, or (at your option) any later version.
// The full text of the license can be found in the file LICENSE at
// the top level of the deal.II distribution.
//
// ---------------------------------------------------------------------
#include <deal.II/algorithms/operator.h>
#include <deal.II/base/logstream.h>
DEAL_II_NAMESPACE_OPEN
namespace Algorithms
{
template <typename VectorType>
OutputOperator<VectorType>::~OutputOperator()
{}
template <typename VectorType>
OutputOperator<VectorType>::OutputOperator()
:
os(0)
{}
template <typename VectorType>
void OutputOperator<VectorType>::initialize_stream(std::ostream &stream)
{
os =&stream;
}
template <typename VectorType>
OutputOperator<VectorType> &
OutputOperator<VectorType>::operator<< (const AnyData &vectors)
{
if (os == 0)
{
deallog << "Step " << step << std::endl;
for (unsigned int i=0; i<vectors.size(); ++i)
{
const VectorType *v = vectors.try_read_ptr<VectorType>(i);
if (v == 0) continue;
deallog << vectors.name(i);
for (unsigned int j=0; j<v->size(); ++j)
deallog << ' ' << (*v)(j);
deallog << std::endl;
}
deallog << std::endl;
}
else
{
(*os) << ' ' << step;
for (unsigned int i=0; i<vectors.size(); ++i)
{
const VectorType *v = vectors.try_read_ptr<VectorType>(i);
if (v == 0) continue;
for (unsigned int j=0; j<v->size(); ++j)
(*os) << ' ' << (*v)(j);
}
(*os) << std::endl;
}
return *this;
}
}
DEAL_II_NAMESPACE_CLOSE
|