Multiscale Universal Interface  2.0
A Concurrent Framework for Coupling Heterogeneous Solvers
sampler.h
Go to the documentation of this file.
1 /*****************************************************************************
2 * Multiscale Universal Interface Code Coupling Library *
3 * *
4 * Copyright (C) 2019 Y. H. Tang, S. Kudo, X. Bian, Z. Li, G. E. Karniadakis *
5 * *
6 * This software is jointly licensed under the Apache License, Version 2.0 *
7 * and the GNU General Public License version 3, you may use it according *
8 * to either. *
9 * *
10 * ** Apache License, version 2.0 ** *
11 * *
12 * Licensed under the Apache License, Version 2.0 (the "License"); *
13 * you may not use this file except in compliance with the License. *
14 * You may obtain a copy of the License at *
15 * *
16 * http://www.apache.org/licenses/LICENSE-2.0 *
17 * *
18 * Unless required by applicable law or agreed to in writing, software *
19 * distributed under the License is distributed on an "AS IS" BASIS, *
20 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. *
21 * See the License for the specific language governing permissions and *
22 * limitations under the License. *
23 * *
24 * ** GNU General Public License, version 3 ** *
25 * *
26 * This program is free software: you can redistribute it and/or modify *
27 * it under the terms of the GNU General Public License as published by *
28 * the Free Software Foundation, either version 3 of the License, or *
29 * (at your option) any later version. *
30 * *
31 * This program is distributed in the hope that it will be useful, *
32 * but WITHOUT ANY WARRANTY; without even the implied warranty of *
33 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
34 * GNU General Public License for more details. *
35 * *
36 * You should have received a copy of the GNU General Public License *
37 * along with this program. If not, see <http://www.gnu.org/licenses/>. *
38 *****************************************************************************/
39 
49 #ifndef MUI_SAMPLER_H_
50 #define MUI_SAMPLER_H_
51 
52 #include "../config.h"
53 #include "../geometry/geometry.h"
54 #include "../storage/virtual_container.h"
55 
56 namespace mui {
57 /*
58 class sampler {
59 public:
60  sampler() {
61  r = 1.0;
62  }
63 
64  double filter( point<double,3> focus, const virtual_container<double> &data_points ) const {
65  double sum = 0.;
66  int n = 0;
67  for( size_t i = 0 ; i < data_points.size() ; i++ ) {
68  if ( ( focus - data_points[i].first ).norm() < r ) {
69  sum += data_points[i].second;
70  n++;
71  }
72  }
73  if (n) return sum / n;
74  else return 0;
75  }
76  inline span<> support() const {
77  return span<>() || geometry::sphere<>( point<double,3>(0), r );
78  }
79 
80 protected:
81  double r;
82 };
83 */
84 }
85 
86 #endif /* MUI_SAMPLER_H_ */
Definition: comm.h:54