The OPAL Discussion Forum

["Finn O'Shea" <finn.oshea AT nusano.com>]

Hi All,

In the course of trying to get Opal to agree with another code, I discovered that the normal component of the octupole strength is set using the formula for decapoles.  Line 114 of Multipole.cpp shows that the OCTUPOLE case is combined with the DECAPOLE case:

    case OCTUPOLE:
    case DECAPOLE:
        NormalComponents[n - 1] = (v + vError) / 24;
        NormalComponentErrors[n - 1] = vError / 24;
        break;

This can be fixed with the following changes:

    case OCTUPOLE:
        NormalComponents[n - 1] = (v + vError) / 6;
        NormalComponentErrors[n - 1] = vError / 6;
        break;
    case DECAPOLE:
        NormalComponents[n - 1] = (v + vError) / 24;
        NormalComponentErrors[n - 1] = vError / 24;
        break;

Finn

--

Senior Scientist

Nusano, Inc

28575 Livingston Ave

Valencia, CA 91355

Office: 1 (424) 293-3174