[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [ccp4bb]: amore problem



***  For details on how to be removed from this list visit the  ***
***          CCP4 home page http://www.ccp4.ac.uk         ***

Tommy Wang wrote:
> 
> ***  For details on how to be removed from this list visit the  ***
> ***          CCP4 home page http://www.ccp4.ac.uk         ***
> 
> Hi,
>         I got the following problem while running amore. Does anyone know what it
> means?
> 
>    98    0.11307E+06    0.87754E+05
>   100    0.41203E+05    0.56024E+05
>  Truncated Patterson Norms:    0.39794E+05    0.19799E+05
>  > S/R FUORI < STEP MISMATCH:     1.765     1.765
> <B><FONT COLOR="#FF0000"><!--SUMMARY_BEGIN-->
>  AMORE:   > s/r fuori < step mismatch:
>  Times: User:    1140.3s System:    0.3s Elapsed:   57:12
> </pre>
> <!--SUMMARY_END--></FONT></B>
> <B><FONT COLOR="#FF0000"><!--SUMMARY_BEGIN-->
> </html>
> <!--SUMMARY_END--></FONT></B>


 This seems to be a rounding error - some poor code here!

Line 5210  amore from ccp4-4.2
      IF (FIRST) THEN
        GMAX = 360.0/NSYM2
        AMAX = 360.0/NSYM1
        STEP = GMAX/NG
C        IF (STEP.NE.AMAX/NA) GO TO 80
        IF (ABS(STEP-AMAX/NA) .GE.0.01) GO TO 80
C
......
C
   80 WRITE (6,FMT=6008) GMAX/NG,AMAX/NA
C
....
 6008 FORMAT (' > S/R FUORI < STEP MISMATCH:',2F10.3)


Log file: 
 > S/R FUORI < STEP MISMATCH:     1.765     1.765

 It gives this error message if GMAX/NG is not equal to AMAX/NA.

However your log file shows they are to 3 decimal places anyway so try
the above change and see if that works.

Eleanor