Clipper
|
00001 00004 //C Copyright (C) 2000-2006 Kevin Cowtan and University of York 00005 //L 00006 //L This library is free software and is distributed under the terms 00007 //L and conditions of version 2.1 of the GNU Lesser General Public 00008 //L Licence (LGPL) with the following additional clause: 00009 //L 00010 //L `You may also combine or link a "work that uses the Library" to 00011 //L produce a work containing portions of the Library, and distribute 00012 //L that work under terms of your choice, provided that you give 00013 //L prominent notice with each copy of the work that the specified 00014 //L version of the Library is used in it, and that you include or 00015 //L provide public access to the complete corresponding 00016 //L machine-readable source code for the Library including whatever 00017 //L changes were used in the work. (i.e. If you make changes to the 00018 //L Library you must distribute those, but you do not need to 00019 //L distribute source or object code to those portions of the work 00020 //L not covered by this licence.)' 00021 //L 00022 //L Note that this clause grants an additional right and does not impose 00023 //L any additional restriction, and so does not affect compatibility 00024 //L with the GNU General Public Licence (GPL). If you wish to negotiate 00025 //L other terms, please contact the maintainer. 00026 //L 00027 //L You can redistribute it and/or modify the library under the terms of 00028 //L the GNU Lesser General Public License as published by the Free Software 00029 //L Foundation; either version 2.1 of the License, or (at your option) any 00030 //L later version. 00031 //L 00032 //L This library is distributed in the hope that it will be useful, but 00033 //L WITHOUT ANY WARRANTY; without even the implied warranty of 00034 //L MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 00035 //L Lesser General Public License for more details. 00036 //L 00037 //L You should have received a copy of the CCP4 licence and/or GNU 00038 //L Lesser General Public License along with this library; if not, write 00039 //L to the CCP4 Secretary, Daresbury Laboratory, Warrington WA4 4AD, UK. 00040 //L The GNU Lesser General Public can also be obtained by writing to the 00041 //L Free Software Foundation, Inc., 59 Temple Place, Suite 330, Boston, 00042 //L MA 02111-1307 USA 00043 00044 00045 #ifndef CLIPPER_HKL_OPERATORS 00046 #define CLIPPER_HKL_OPERATORS 00047 00048 #include "hkl_datatypes.h" 00049 00050 00051 namespace clipper 00052 { 00053 00054 namespace datatypes 00055 { 00056 00057 // logical operators 00058 HKL_data<Flag_bool> operator &( const HKL_data_base& d1, const HKL_data_base& d2 ); 00059 HKL_data<Flag_bool> operator |( const HKL_data_base& d1, const HKL_data_base& d2 ); 00060 HKL_data<Flag_bool> operator ^( const HKL_data_base& d1, const HKL_data_base& d2 ); 00061 HKL_data<Flag_bool> operator !( const HKL_data_base& d1 ); 00062 00063 HKL_data<Flag_bool> operator ==( const HKL_data<Flag>& d1, const int& n ); 00064 HKL_data<Flag_bool> operator !=( const HKL_data<Flag>& d1, const int& n ); 00065 HKL_data<Flag_bool> operator >=( const HKL_data<Flag>& d1, const int& n ); 00066 HKL_data<Flag_bool> operator <=( const HKL_data<Flag>& d1, const int& n ); 00067 HKL_data<Flag_bool> operator >( const HKL_data<Flag>& d1, const int& n ); 00068 HKL_data<Flag_bool> operator <( const HKL_data<Flag>& d1, const int& n ); 00069 00070 // individual data operators 00071 template<class dtype> F_phi<dtype> operator +( const F_phi<dtype>& d1, const F_phi<dtype>& d2 ); 00072 template<class dtype> F_phi<dtype> operator -( const F_phi<dtype>& d1, const F_phi<dtype>& d2 ); 00073 template<class dtype> F_phi<dtype> operator -( const F_phi<dtype>& d1 ); 00074 template<class dtype> ABCD<dtype> operator +( const ABCD<dtype>& d1, const ABCD<dtype>& d2 ); 00075 00076 // data list operators 00077 template<class dtype> HKL_data<F_phi<dtype> > operator +( const HKL_data<F_phi<dtype> >& d1, const HKL_data<F_phi<dtype> >& d2 ); 00078 template<class dtype> HKL_data<F_phi<dtype> > operator -( const HKL_data<F_phi<dtype> >& d1, const HKL_data<F_phi<dtype> >& d2 ); 00079 template<class dtype> HKL_data<F_phi<dtype> > operator *( const HKL_data<F_phi<dtype> >& d1, const ftype& s ); 00080 template<class dtype> HKL_data<F_phi<dtype> > operator -( const HKL_data<F_phi<dtype> >& d1 ); 00081 template<class dtype> HKL_data<ABCD<dtype> > operator +( const HKL_data<ABCD<dtype> >& d1, const HKL_data<ABCD<dtype> >& d2 ); 00082 template<class dtype> HKL_data<F_phi<dtype> > operator *( const ftype& s, const HKL_data<F_phi<dtype> >& d1 ) { return d1*s; } 00083 00084 } // namespace datatypes 00085 00086 00087 00089 00096 template<int N> class LogPhaseProb { 00097 public: 00099 LogPhaseProb( const HKL_class& hkl_class ); 00101 template<class dtype> void set_abcd( const datatypes::ABCD<dtype>& abcd ); 00103 template<class dtype> void get_abcd( datatypes::ABCD<dtype>& abcd ) const; 00105 template<class dtype> void set_phi_fom( const datatypes::Phi_fom<dtype>& phifom ); 00107 template<class dtype> void get_phi_fom( datatypes::Phi_fom<dtype>& phifom ) const; 00109 const ftype& operator[] ( const int& p ) const { return q[p]; } 00111 ftype& operator[] ( const int& p ) { return q[p]; } 00113 ftype phase( const int& p ) const 00114 { return Util::twopi()*ftype(p*pinc+pmin)/ftype(N); } 00115 int size() const { return q.size(); } 00116 static int sampling() { return N; } 00117 private: 00118 int pmin, pinc; 00119 std::vector<ftype> q; 00120 }; 00121 00122 00123 } // namespace clipper 00124 00125 #endif