/usr/include/dolfin/swig/fem_post.i is in libdolfin1.0-dev 1.0.0-1.
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 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 | /* -*- C -*- */
// Copyright (C) 2009 Johan Hake
//
// This file is part of DOLFIN.
//
// DOLFIN is free software: you can 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 3 of the License, or
// (at your option) any later version.
//
// DOLFIN is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
// GNU Lesser General Public License for more details.
//
// You should have received a copy of the GNU Lesser General Public License
// along with DOLFIN. If not, see <http://www.gnu.org/licenses/>.
//
// First added: 2007-10-22
// Last changed: 2011-08-18
// ===========================================================================
// SWIG directives for the DOLFIN fem kernel module (post)
//
// The directives in this file are applied _after_ the header files of the
// modules has been loaded.
// ===========================================================================
//-----------------------------------------------------------------------------
// Extend Function so f.function_space() return a dolfin.FunctionSpace
//-----------------------------------------------------------------------------
%extend dolfin::BoundaryCondition {
%pythoncode %{
def function_space(self):
"Return the FunctionSpace"
from dolfin.functions.functionspace import FunctionSpaceFromCpp
return FunctionSpaceFromCpp(self._function_space())
%}
}
//-----------------------------------------------------------------------------
// Extend GenericDofMap.tabulate_coordinates()
//-----------------------------------------------------------------------------
%extend dolfin::GenericDofMap {
void _tabulate_coordinates(PyObject* coordinates, const Cell& cell)
{
// NOTE: No NumPy array check. Assumed everything is coorect!
// Get NumPy array
PyArrayObject *xa = reinterpret_cast<PyArrayObject*>(coordinates);
// Initialize the boost::multi_array structure
boost::multi_array<double, 2>::extent_gen extents;
boost::multi_array<double, 2> tmparray;
tmparray.resize(extents[self->cell_dimension(cell.index())]\
[self->geometric_dimension()]);
// Tabulate the coordinates
dolfin::UFCCell ufc_cell(cell);
self->tabulate_coordinates(tmparray, ufc_cell);
// Copy data
double* data = static_cast<double*>(PyArray_DATA(xa));
for (dolfin::uint i=0; i<self->cell_dimension(cell.index()); i++)
for (dolfin::uint j=0; j<self->geometric_dimension(); j++)
data[i*self->geometric_dimension()+j] = tmparray[i][j];
}
%pythoncode %{
def tabulate_coordinates(self, cell, coordinates=None):
""" Tabulate the coordinates of all dofs on a cell
*Arguments*
cell (_Cell_)
The cell.
coordinates (NumPy array)
Optional argument: The coordinates of all dofs on a cell.
*Returns*
coordinates
The coordinates of all dofs on a cell.
"""
import numpy as np
# Check coordinate argument
shape = (self.max_cell_dimension(), self.geometric_dimension())
if coordinates is None:
coordinates = np.zeros(shape, 'd')
if not isinstance(coordinates, np.ndarray) or \
not (coordinates.flags.c_contiguous and \
coordinates.dtype == np.dtype('d') and \
coordinates.shape==shape):
raise TypeError, "expected a C-contiguous numpy array " \
"of 'double' (dtype='d') with shape %s"%str(shape)
# Call the extended method
self._tabulate_coordinates(coordinates, cell)
return coordinates
%}
}
|