Clipper
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 
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:
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:
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