C++********************************************************************* C C APPLYW.F C C ********************************************************************** C=* FROM: SPIDER - MODULAR IMAGE PROCESSING SYSTEM. AUTHOR: J.FRANK * C=* Copyright (C) 1985-2005 Health Research Inc. * C=* * C=* HEALTH RESEARCH INCORPORATED (HRI), * C=* ONE UNIVERSITY PLACE, RENSSELAER, NY 12144-3455. * C=* * C=* Email: spider@wadsworth.org * C=* * C=* This program is free software; you can redistribute it and/or * C=* modify it under the terms of the GNU General Public License as * C=* published by the Free Software Foundation; either version 2 of the * C=* License, or (at your option) any later version. * C=* * C=* This program is distributed in the hope that it will be useful, * C=* but WITHOUT ANY WARRANTY; without even the implied warranty of * C=* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU * C=* General Public License for more details. * C=* * C=* You should have received a copy of the GNU General Public License * C=* along with this program; if not, write to the * C=* Free Software Foundation, Inc., * C=* 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA. * C=* * C ********************************************************************** C C IMAGE_PROCESSING_ROUTINE C C 1 2 3 4 5 6 7 C--********************************************************************* SUBROUTINE APPLYW(CIRC,LCIRC,NUMR,WR,NRING,MAXRIN) INTEGER NUMR(3,NRING),MAXRIN,NUMR3I,NUMR2I DIMENSION CIRC(LCIRC),WR(NRING) c$omp parallel do private(i,j,jc,NUMR3I,NUMR2I,W) DO I=1,NRING NUMR3I=NUMR(3,I) NUMR2I=NUMR(2,I) W=WR(I) CIRC(NUMR2I)=CIRC(NUMR2I)*W IF (NUMR3I.EQ.MAXRIN) THEN CIRC(NUMR2I+1)=CIRC(NUMR2I+1)*W ELSE CIRC(NUMR2I+1)=CIRC(NUMR2I+1)*0.5*W ENDIF DO J=3,NUMR3I JC=J+NUMR2I-1 CIRC(JC)=CIRC(JC)*W ENDDO ENDDO END