Clipper
Main Page
Related Pages
Namespaces
Classes
Files
File List
clipper
core
nxmap_operator.h
1
4
//C Copyright (C) 2000-2006 Kevin Cowtan and University of York
5
//L
6
//L This library is free software and is distributed under the terms
7
//L and conditions of version 2.1 of the GNU Lesser General Public
8
//L Licence (LGPL) with the following additional clause:
9
//L
10
//L `You may also combine or link a "work that uses the Library" to
11
//L produce a work containing portions of the Library, and distribute
12
//L that work under terms of your choice, provided that you give
13
//L prominent notice with each copy of the work that the specified
14
//L version of the Library is used in it, and that you include or
15
//L provide public access to the complete corresponding
16
//L machine-readable source code for the Library including whatever
17
//L changes were used in the work. (i.e. If you make changes to the
18
//L Library you must distribute those, but you do not need to
19
//L distribute source or object code to those portions of the work
20
//L not covered by this licence.)'
21
//L
22
//L Note that this clause grants an additional right and does not impose
23
//L any additional restriction, and so does not affect compatibility
24
//L with the GNU General Public Licence (GPL). If you wish to negotiate
25
//L other terms, please contact the maintainer.
26
//L
27
//L You can redistribute it and/or modify the library under the terms of
28
//L the GNU Lesser General Public License as published by the Free Software
29
//L Foundation; either version 2.1 of the License, or (at your option) any
30
//L later version.
31
//L
32
//L This library is distributed in the hope that it will be useful, but
33
//L WITHOUT ANY WARRANTY; without even the implied warranty of
34
//L MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
35
//L Lesser General Public License for more details.
36
//L
37
//L You should have received a copy of the CCP4 licence and/or GNU
38
//L Lesser General Public License along with this library; if not, write
39
//L to the CCP4 Secretary, Daresbury Laboratory, Warrington WA4 4AD, UK.
40
//L The GNU Lesser General Public can also be obtained by writing to the
41
//L Free Software Foundation, Inc., 59 Temple Place, Suite 330, Boston,
42
//L MA 02111-1307 USA
43
44
45
#ifndef CLIPPER_NXMAP_OPERATOR
46
#define CLIPPER_NXMAP_OPERATOR
47
48
#include "xmap.h"
49
#include "nxmap.h"
50
51
52
namespace
clipper
53
{
54
56
67
class
NX_operator
68
{
69
public
:
71
NX_operator
();
73
NX_operator
(
const
Xmap_base
& xmap,
const
NXmap_base
& nxmap,
const
RTop_orth
& rtop );
75
NX_operator
(
const
Cell
& cell,
const
Grid_sampling
& grid,
const
NXmap_base
& nxmap,
const
RTop_orth
& rtop );
77
void
init
(
const
Xmap_base
& xmap,
const
NXmap_base
& nxmap,
const
RTop_orth
& rtop );
79
void
init
(
const
Cell
& cell,
const
Grid_sampling
& grid,
const
NXmap_base
& nxmap,
const
RTop_orth
& rtop );
80
82
inline
Coord_map
coord_map
(
const
Coord_frac
& c )
const
83
{
return
Coord_map
(
xfrac_nxgrid
* c ); }
85
inline
Coord_frac
coord_frac
(
const
Coord_map
& c )
const
86
{
return
Coord_frac
(
nxgrid_xfrac
* c ); }
88
template
<
class
I,
class
T,
class
M> T
nxmap_data
(
const
M& nxmap,
const
Coord_grid
& c )
const
;
90
template
<
class
I,
class
T,
class
M> T
xmap_data
(
const
M& xmap,
const
Coord_grid
& c )
const
;
91
93
bool
is_null
()
const
;
94
95
void
debug()
const
;
96
97
protected
:
98
RTop<>
xfrac_nxgrid
;
99
RTop<>
nxgrid_xfrac
;
100
RTop<>
xgrid_nxgrid
;
101
RTop<>
nxgrid_xgrid
;
102
RTop<int>
xgrid_nxgrid_int
;
103
RTop<int>
nxgrid_xgrid_int
;
104
bool
x_nx_is_int
;
105
bool
x_nx_is_trn
;
106
bool
nx_x_is_int
;
107
bool
nx_x_is_trn
;
108
};
109
110
112
127
template
<
class
T>
class
NXmap_operator
:
public
NX_operator
128
{
129
public
:
131
NXmap_operator
() {}
133
NXmap_operator
(
const
Xmap_base
& xmap,
const
NXmap<T>
&
nxmap
,
const
RTop_orth
& rtop ) {
init
( xmap, nxmap, rtop ); }
135
NXmap_operator
(
const
Cell
& cell,
const
Grid_sampling
& grid,
const
NXmap<T>
&
nxmap
,
const
RTop_orth
& rtop ) {
init
( cell, grid, nxmap, rtop ); }
137
void
init
(
const
Xmap_base
& xmap,
const
NXmap<T>
&
nxmap
,
const
RTop_orth
& rtop ) {
init
( xmap.
cell
(), xmap.
grid_sampling
(),
nxmap
, rtop ); }
139
void
init
(
const
Cell
& cell,
const
Grid_sampling
& grid,
const
NXmap<T>
&
nxmap
,
const
RTop_orth
& rtop ) { nxmap_ = &
nxmap
;
NX_operator::init
( cell, grid, nxmap, rtop ); }
140
142
template
<
class
I> T
nxmap_data
(
const
Coord_grid
& c )
const
143
{
return
NX_operator::nxmap_data<I,T>( *nxmap_, c ); }
144
146
const
NXmap<T>
&
nxmap
()
const
{
return
*nxmap_; }
147
148
private
:
149
const
NXmap<T>
* nxmap_;
150
};
151
152
153
154
// template implementations
155
165
template
<
class
I,
class
T,
class
M> T
NX_operator::nxmap_data
(
const
M& nxmap,
const
Coord_grid
& c )
const
166
{
167
if
(
x_nx_is_trn
) {
168
return
T( nxmap.get_data(
Coord_grid
( c +
xgrid_nxgrid_int
.
trn
() ) ) );
169
}
else
if
(
x_nx_is_int
) {
170
return
T( nxmap.get_data(
Coord_grid
(
xgrid_nxgrid_int
* c ) ) );
171
}
else
{
172
T val;
173
I::interp( nxmap,
Coord_map
(
xgrid_nxgrid
* c.
coord_map
() ), val );
174
return
val;
175
}
176
}
177
178
187
template
<
class
I,
class
T,
class
M> T
NX_operator::xmap_data
(
const
M& xmap,
const
Coord_grid
& c )
const
188
{
189
if
(
nx_x_is_trn
) {
190
return
T( xmap.get_data(
Coord_grid
( c +
nxgrid_xgrid_int
.
trn
() ) ) );
191
}
else
if
(
nx_x_is_int
) {
192
return
T( xmap.get_data(
Coord_grid
(
nxgrid_xgrid_int
* c ) ) );
193
}
else
{
194
T val;
195
I::interp( xmap,
Coord_map
(
nxgrid_xgrid
* c.
coord_map
() ), val );
196
return
val;
197
}
198
}
199
200
201
}
// namespace clipper
202
203
#endif
Generated on Fri Sep 6 2013 19:58:36 for Clipper by
1.8.4