00001 /* 00002 * Author: Steven Ludtke, 04/10/2003 (sludtke@bcm.edu) 00003 * Copyright (c) 2000-2006 Baylor College of Medicine 00004 * 00005 * This software is issued under a joint BSD/GNU license. You may use the 00006 * source code in this file under either license. However, note that the 00007 * complete EMAN2 and SPARX software packages have some GPL dependencies, 00008 * so you are responsible for compliance with the licenses of these packages 00009 * if you opt to use BSD licensing. The warranty disclaimer below holds 00010 * in either instance. 00011 * 00012 * This complete copyright notice must be included in any revised version of the 00013 * source code. Additional authorship citations may be added, but existing 00014 * author citations must be preserved. 00015 * 00016 * This program is free software; you can redistribute it and/or modify 00017 * it under the terms of the GNU General Public License as published by 00018 * the Free Software Foundation; either version 2 of the License, or 00019 * (at your option) any later version. 00020 * 00021 * This program is distributed in the hope that it will be useful, 00022 * but WITHOUT ANY WARRANTY; without even the implied warranty of 00023 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 00024 * GNU General Public License for more details. 00025 * 00026 * You should have received a copy of the GNU General Public License 00027 * along with this program; if not, write to the Free Software 00028 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA 00029 * 00030 * */ 00031 00032 #include "cmp_template.h" 00033 #include "transform.h" 00034 00035 00036 using namespace EMAN; 00037 00038 const string XYZCmp::NAME = "middle"; 00039 00040 float XYZCmp::cmp(EMData * image, EMData * with) const 00041 { 00042 validate_input_args(image, with); 00043 00044 #if 0 00045 int param1 = params["param1"]; 00046 float param2 = params["param2"]; 00047 // do something 00048 #endif 00049 int nx=image->get_xsize(); 00050 int ny=image->get_ysize(); 00051 00052 double ret=0; 00053 for (int i=nx/4; i<nx*3/4; i++) { 00054 for (int j=ny/4; j<ny*3/4; j++) { 00055 ret+=image->get_value_at(i,j)*with->get_value_at(i,j); 00056 } 00057 } 00058 00059 return (float) ret; 00060 }