| 1 | 
+ | 
#define _LARGEFILE_SOURCE64 | 
| 2 | 
  | 
#define _FILE_OFFSET_BITS 64 | 
| 3 | 
  | 
 | 
| 4 | 
  | 
#include <string.h> | 
| 227 | 
  | 
  int nAtoms = entry_plug->n_atoms; | 
| 228 | 
  | 
#endif //is_mpi | 
| 229 | 
  | 
 | 
| 230 | 
< | 
  double q[4]; | 
| 230 | 
> | 
  double q[4], ji[3]; | 
| 231 | 
  | 
  DirectionalAtom* dAtom; | 
| 232 | 
  | 
  Atom** atoms = entry_plug->atoms; | 
| 233 | 
  | 
  double pos[3], vel[3]; | 
| 274 | 
  | 
 | 
| 275 | 
  | 
      dAtom = (DirectionalAtom *)atoms[i]; | 
| 276 | 
  | 
      dAtom->getQ( q ); | 
| 277 | 
+ | 
      dAtom->getJ( ji ); | 
| 278 | 
  | 
 | 
| 279 | 
  | 
      sprintf( tempBuffer, | 
| 280 | 
  | 
               "%lf\t%lf\t%lf\t%lf\t%lf\t%lf\t%lf\n", | 
| 282 | 
  | 
               q[1], | 
| 283 | 
  | 
               q[2], | 
| 284 | 
  | 
               q[3], | 
| 285 | 
< | 
               dAtom->getJx(), | 
| 286 | 
< | 
               dAtom->getJy(), | 
| 287 | 
< | 
               dAtom->getJz()); | 
| 285 | 
> | 
               ji[0], | 
| 286 | 
> | 
               ji[1], | 
| 287 | 
> | 
               ji[2]); | 
| 288 | 
  | 
      strcat( writeLine, tempBuffer ); | 
| 289 | 
  | 
    } | 
| 290 | 
  | 
    else | 
| 417 | 
  | 
             | 
| 418 | 
  | 
            dAtom = (DirectionalAtom *)atoms[local_index]; | 
| 419 | 
  | 
            dAtom->getQ( q ); | 
| 420 | 
+ | 
            dAtom->getJ( ji ); | 
| 421 | 
  | 
 | 
| 422 | 
  | 
            for (int j = 0; j < 6 ; j++) | 
| 423 | 
  | 
              atomData13[j] = atomData6[j];             | 
| 427 | 
  | 
            atomData13[8] = q[2]; | 
| 428 | 
  | 
            atomData13[9] = q[3]; | 
| 429 | 
  | 
             | 
| 430 | 
< | 
            atomData13[10] = dAtom->getJx(); | 
| 431 | 
< | 
            atomData13[11] = dAtom->getJy(); | 
| 432 | 
< | 
            atomData13[12] = dAtom->getJz(); | 
| 430 | 
> | 
            atomData13[10] = ji[0]; | 
| 431 | 
> | 
            atomData13[11] = ji[1]; | 
| 432 | 
> | 
            atomData13[12] = ji[2]; | 
| 433 | 
  | 
          } | 
| 434 | 
  | 
           | 
| 435 | 
  | 
        } else { | 
| 546 | 
  | 
             | 
| 547 | 
  | 
            dAtom = (DirectionalAtom *)atoms[local_index]; | 
| 548 | 
  | 
            dAtom->getQ( q ); | 
| 549 | 
+ | 
            dAtom->getJ( ji ); | 
| 550 | 
  | 
             | 
| 551 | 
  | 
            for (int j = 0; j < 6 ; j++) | 
| 552 | 
  | 
              atomData13[j] = atomData6[j]; | 
| 556 | 
  | 
            atomData13[8] = q[2]; | 
| 557 | 
  | 
            atomData13[9] = q[3]; | 
| 558 | 
  | 
   | 
| 559 | 
< | 
            atomData13[10] = dAtom->getJx(); | 
| 560 | 
< | 
            atomData13[11] = dAtom->getJy(); | 
| 561 | 
< | 
            atomData13[12] = dAtom->getJz(); | 
| 559 | 
> | 
            atomData13[10] = ji[0]; | 
| 560 | 
> | 
            atomData13[11] = ji[1]; | 
| 561 | 
> | 
            atomData13[12] = ji[2]; | 
| 562 | 
  | 
          } | 
| 563 | 
  | 
 | 
| 564 | 
  | 
        } else { |