3. Getting started

Have a quick look at the test program tracytest.c ( Source):

#include <math.h>

#include "tracy.h"
#include "datatyp.h"
#include "physlib.h"

main() {
  boolean chroma;
  double dP;
  globvalrec globv;
  long lastpos;

  t2init();
  
  fi=fopen("sls.lat","r"); /* open the lattice input file  */
  fo=fopen("sls.out","w"); /* open the lattice output file */

  if (Lattice_Read(&fi, &fo)) { /* read lattice */
    Cell_Init(); /* initialize cell structure */

    getglobv_(&globv);

    /* define x/y aperture limits */
    for (i = 0; i <= globval.Cell_nLoc; i++)
      globv.maxampl[i][0] = globv.maxampl[i][1] = 1e0;

    globv.MatMeth = false; globv.Cavity_on = false;
    globv.radiation = false; globv.emittance = false;
    globv.pathlength = false; globv.CODimax = 15;
    globv.bpm = ElemIndex("mon");

    putglobv_(&globv);

    chroma=true;
    dP=0.0;

    getglobv_(&globv);
    printglob();

    Ring_GetTwiss(chroma, dP); /* get Twiss parameters */

    printglob(); /* print parameter list */
    printlatt(); /* dump linear lattice functions into "linlat.dat" */

    getcod(0.0, &lastpos); /* determine closed orbit */
    printcod(); /* dump closed orbit into "cod.dat" */
  }
}