1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174
|
/* Copyright (C) 2000-2003 Damir Zucic */
/*=============================================================================
center.c
Purpose:
Execute center command: calculate the geometric center position
for currently selected atoms and use this position as a new
origin of the internal coordinate system. Only one macromol.
complex should be caught, otherwise this function will fail to
work.
Input:
(1) Pointer to MolComplexS structure.
(2) The number of macromolecular complexes.
(3) Pointer to RuntimeS structure.
(4) Pointer to ConfigS structure.
(5) Pointer to GUIS structure.
(6) Pointer to NearestAtomS structure.
(7) The number of pixels in the main window free area.
(8) Pointer to refreshI.
Output:
(1)
(2) Return value.
Return value:
(1) Positive (command) code on success.
(2) Negative (error) code on failure.
========includes:============================================================*/
#include <stdio.h>
#include <string.h>
#include <X11/Xlib.h>
#include <X11/Xutil.h>
#include <X11/Xos.h>
#include <X11/Xatom.h>
#include "defines.h"
#include "commands.h"
#include "typedefs.h"
/*======function prototypes:=================================================*/
size_t MainRefresh_ (MolComplexS *, int,
RuntimeS *, ConfigS *, GUIS *,
NearestAtomS *, size_t, unsigned int);
int ControlRefresh_ (MolComplexS *, ConfigS *, GUIS *);
/*======execute center command:==============================================*/
int Center_ (MolComplexS *mol_complexSP, int mol_complexesN,
RuntimeS *runtimeSP, ConfigS *configSP, GUIS *guiSP,
NearestAtomS *nearest_atomSP, size_t pixelsN,
unsigned int *refreshIP)
{
int mol_complexI;
int caught_complexesN = 0;
MolComplexS *curr_mol_complexSP;
size_t atomsN, atomI;
AtomS *curr_atomSP;
size_t selected_atomsN = 0;
double x = 0.0, y = 0.0, z = 0.0, denom;
/* Count the number of currently caught macromolecular complexes: */
for (mol_complexI = 0; mol_complexI < mol_complexesN; mol_complexI++)
{
if ((mol_complexSP + mol_complexI)->catchF == 1) caught_complexesN++;
}
/* Check the number of cought complexes (must be equal to one): */
if (caught_complexesN != 1)
{
strcpy (runtimeSP->messageA,
"One (and only one) complex should be caught!");
runtimeSP->message_length = strlen (runtimeSP->messageA);
return ERROR_CENTER;
}
/* Find the cought complex: */
for (mol_complexI = 0; mol_complexI < mol_complexesN; mol_complexI++)
{
/** Pointer to the current macromolecular complex: **/
curr_mol_complexSP = mol_complexSP + mol_complexI;
/** Check is the current macromolecular complex caught: **/
if (curr_mol_complexSP->catchF == 0) continue;
/** Number of atoms in a macromolecular complex: **/
atomsN = curr_mol_complexSP->atomsN;
if (atomsN == 0)
{
strcpy (runtimeSP->messageA,
"There are no atoms in this macromol. complex!");
runtimeSP->message_length = strlen (runtimeSP->messageA);
return ERROR_CENTER;
}
/** Scan all atoms in the current complex: **/
for (atomI = 0; atomI < atomsN; atomI++)
{
/*** Pointer to the current atom: ***/
curr_atomSP = curr_mol_complexSP->atomSP + atomI;
/*** Check is this atom selected: ***/
if (curr_atomSP->selectedF == 0) continue;
/*** Update the number of selected atoms: ***/
selected_atomsN++;
/*** Prepare data required to find geometric center: ***/
x += curr_atomSP->raw_atomS.x[0];
y += curr_atomSP->raw_atomS.y;
z += curr_atomSP->raw_atomS.z[0];
}
/** Check how many atoms are selected: **/
if (selected_atomsN == 0)
{
sprintf (runtimeSP->messageA, "There are no selected atoms!");
runtimeSP->message_length = strlen (runtimeSP->messageA);
return ERROR_CENTER;
}
/** Calculate the geometric center position: **/
denom = 1.0 / (double) selected_atomsN;
x *= denom;
y *= denom;
z *= denom;
/** Update the geometric center position: **/
curr_mol_complexSP->geometric_center_vectorS.x = x;
curr_mol_complexSP->geometric_center_vectorS.y = y;
curr_mol_complexSP->geometric_center_vectorS.z = z;
/** Update the rotation center position: **/
curr_mol_complexSP->rotation_center_vectorS.x = x;
curr_mol_complexSP->rotation_center_vectorS.y = y;
curr_mol_complexSP->rotation_center_vectorS.z = z;
/** Update the slab center: **/
curr_mol_complexSP->slab_center_vectorS.x = x;
curr_mol_complexSP->slab_center_vectorS.y = y;
curr_mol_complexSP->slab_center_vectorS.z = z;
/** Update the color fading center: **/
curr_mol_complexSP->fading_center_vectorS.x = x;
curr_mol_complexSP->fading_center_vectorS.y = y;
curr_mol_complexSP->fading_center_vectorS.z = z;
/* Update the position_changedF: */
curr_mol_complexSP->position_changedF = 1;
}
/* Refresh the main wind w: */
(*refreshIP)++;
MainRefresh_ (mol_complexSP, mol_complexesN,
runtimeSP, configSP, guiSP,
nearest_atomSP, pixelsN, *refreshIP);
/* Refresh the control window: */
ControlRefresh_ (mol_complexSP + runtimeSP->default_complexI, configSP, guiSP);
/* Return the command code: */
return COMMAND_CENTER;
}
/*===========================================================================*/
|