NiHu
2.0
helmholtz_bem_3d_surf.mex.cpp
1
// This file is a part of NiHu, a C++ BEM template library.
2
//
3
// Copyright (C) 2012-2014 Peter Fiala <fiala@hit.bme.hu>
4
// Copyright (C) 2012-2014 Peter Rucz <rucz@hit.bme.hu>
5
//
6
// This program is free software: you can redistribute it and/or modify
7
// it under the terms of the GNU General Public License as published by
8
// the Free Software Foundation, either version 3 of the License, or
9
// (at your option) any later version.
10
//
11
// This program is distributed in the hope that it will be useful,
12
// but WITHOUT ANY WARRANTY; without even the implied warranty of
13
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14
// GNU General Public License for more details.
15
//
16
// You should have received a copy of the GNU General Public License
17
// along with this program. If not, see <http://www.gnu.org/licenses/>.
18
20
#include "
core/weighted_residual.hpp
"
21
#include "
util/mex_matrix.hpp
"
22
#include "library/helmholtz_3d.hpp"
23
#include "
library/lib_element.hpp
"
24
25
typedef
NiHu::mex::real_matrix<double>
dMatrix
;
26
typedef
NiHu::mex::complex_matrix<double>
cMatrix;
28
30
void
mexFunction(
int
nlhs, mxArray *lhs[],
int
nrhs, mxArray
const
*rhs[])
31
{
32
dMatrix
surf_nodes(rhs[0]), surf_elem(rhs[1]);
33
auto
surf_mesh =
NiHu::create_mesh
(surf_nodes, surf_elem,
NiHu::quad_1_tag
(),
NiHu::tria_1_tag
());
35
37
auto
const
&surf_sp =
NiHu::constant_view
(surf_mesh);
39
41
size_t
n = surf_sp.get_num_dofs();
42
cMatrix Ls(n, n, lhs[0]), Ms(n, n, lhs[1]);
44
46
double
k = *mxGetPr(rhs[2]);
47
auto
I =
NiHu::identity_integral_operator
();
48
auto
L =
NiHu::create_integral_operator
(
NiHu::helmholtz_3d_SLP_kernel<double>
(k));
49
auto
M =
NiHu::create_integral_operator
(
NiHu::helmholtz_3d_DLP_kernel<double>
(k));
51
53
Ls << surf_sp * L[surf_sp];
54
Ms << surf_sp * M[surf_sp] + surf_sp * (-.5*I)[surf_sp];
55
}
57
NiHu::create_integral_operator
integral_operator< Kernel > create_integral_operator(Kernel &&kernel)
factory function of an integral operator
Definition:
integral_operator.hpp:421
NiHu::type2tag
Metafunction assigning a tag to a type.
Definition:
type2tag.hpp:17
NiHu::create_mesh
mesh< tmp::vector< typename tag2type< Args >::type... > > create_mesh(nodes_t const &nodes, elements_t const &elements, Args...)
factory function to create a mesh from nodes and elements matrices
Definition:
mesh.hpp:298
NiHu::identity_integral_operator
The identity integral operator .
Definition:
integral_operator.hpp:310
NiHu::constant_view
const function_space_view< Mesh, field_option::constant, Dimension > & constant_view(mesh_base< Mesh > const &msh, Dimension dim=Dimension())
factory function to transform a mesh into a constant function space
Definition:
function_space.hpp:304
weighted_residual.hpp
declaration of class NiHu::weighted_residual
NiHu::mex::complex_matrix
Container class of a complex matrix stored in Matlab format.
Definition:
mex_matrix_separate.hpp:256
mex_matrix.hpp
A Matlab mex matrix interface.
NiHu::mex::real_matrix
Definition:
mex_matrix_separate.hpp:133
NiHu::normal_derivative_kernel
Normal derivative of a distance dependent kernel.
Definition:
normal_derivative_kernel.hpp:26
lib_element.hpp
src
tutorial
helmholtz_bem_3d_surf.mex.cpp
Generated on Sun May 5 2024 01:00:17 for NiHu by
1.8.18