Dear colleagues.
A new version has been checked in into svn and GitHub devel repositories. Several fixes pointed out by A. Morita, S. Yamada and others have been included. Besides them major changes are going on in tracking (classical spin tracking and radiation calculation) and emittance (a kinematical algorithm). This version will be much buggier than usual, so please do not use this seriously yet. We need more time to debug.
Fixes:
1. A bug in qsol to have used qcopymat instead of qcopymatg (pointed out by A. Morita). An acceleration in SOL by MULT has been wrong.
Changes:
1. Use S. Yamada's patch on utils.c to avoid piling up of temp files. He saids that using a pipe instead of a file was better to receive results of shell commands, but it has not been included so far, due to a bad experience of mine on pipe in an old system, probably HP-UX.
2. DROTATE for BEND has been implemented in a more precise way without approximation, in both tracking and emittance/optics.
3. Tracking of classical spin vectors has been implemented. It is usable by TrackParticles with 8 variables under the flag POL. The 7th variable represent the phase angle of (s_x, s_z) and the 8th is s_z. s_{x,y,z} must be in the range [-1, 1]. These new variables are symplectic ones. Currently no algorithm like SLIM has not been implemented yet.
4. Added Nominal spin tune and Nominal polarization time to emitttance output. They are accessible by Emittance[], too.
5. Changed the algorithm to calculate radiation effects in tracking and emittance calculation by a kinematical method in all elements including fringes and multipoles. This algorithm can calculate radiations without knowing the detailed information on the field, and has been already applied to the tracking. However, several effects such as F1 of a dipole are not yet nicely included. It is under development and we shoud wait for the next year.
> Changes:
> 1. Use S. Yamada's patch on utils.c to avoid piling up of temp files. He saids that using a pipe instead of a file was better to receive results of shell commands, but it has not been included so far, due to a bad experience of mine on pipe in an old system, probably HP-UX.
>
この"bad experience"というのは、backgroundな sub-processでは無く、foreground sub-processの出力先にpipeを使ったケースではないでしょうか？
pipeへのwrite(2)は、kernelのI/O bufferが埋まるとreader側が読みにくるまでblocking waitに入るので、foreground sub-processの出力をpipeへ向けた場合、write(2)のblockingによりsub-processが終了しないために、pipeのreaderである main-processも進行せずにデッドロックします
MAIN trunk系の実装は、sub-processがforeground時はunlink済みのファイル記述子、sub-processがbackground時はpipe(2)が生成するファイル記述子を使って実装しています(参考までに)
Dear Users,
Again, this is an experimental version, and do not use seriously.
Changes:
1. Added a linear spin depolarization calculation in emit/Emittance[]. However, it does not look right yet, and not consistent with tracking.
2. The number of parallel processes is determined by NPARA, number of particles and the length of the beamline.
Dear Users,
Fixes:
1. Corrected handling of spin in emittance calculation, esp. in the slicing of tsolque and texchg.
2. Fixed a mistake on rotation in tchge.
Changes:
1. Added spin depolarization calculatioin with resonance perturbation (preliminary). Currently coded up to 4th order, but routines are extendable to any order.
Fixes:
1. Bugs in optics, geometry, tracking in SOL with BEND and QUAD, after spin was included.
Changes:
The version number has jumped.
Fixes:
1. Tracking with RAD was not properly handled at the end of QUAD and MULT to bring slightly larger energy loss, since "kinematical evaluation" was introduced.
Changes:
1. Calculate polarization rate at each element.
2. A few modification in the polarization calculation, but the results are not affected.
3. The radiation fluctuations per turn are halved for depolarization, since it is generated in the middle of one turn in average.
4. Rectified depolarization calculation.
5. Reduced the number of slices of a BEND from 0.1 photons to 1 photon per slice in emittance calculation. Also 1 photon/slice for 1 cm amplitude for MULT.
7. Set number of slices in tracking for BEND, QUAD, and MULT more consistently. Now the number of photons per slice is less than 1 (RING) ot 0.1 (TRPT), assuming r=0.02 m.
Fixed:
1. Confusion of INF and NaN for output.
2. A glitch when set-element-value-command is applied to multiple types.
Added:
1. A new function Put (similar to Mathematica):
Put[x__,file_String]:=With[{f=OpenWrite[file]},
Scan[Write[f,#]&,{x}];
Close[f]];
2. New tags Polarization2, Polarization4, Polarization6 in the results of Emittance[], corresponding to 2nd, 4th, 6th order calculation, respectively.
Polarization returs the same value as Polarization.
3. A new function FindClosedOrbit[z0_:{0,0,0,0}]: try to find out the closed orbit with CALC4D. z0 is the initial orbit at the entrance. CALC4D only.
Changes:
1. Eleminate a list as the arg of FFS USE .
2. The internal indices of idval, idtype, pname, lpname now start from 0.
3. Added cross terms into 4th and 6th order terms in the depolarization matrix.
4. Now element-value-command such as "Q* K1 0" allows a wildcard for the keyword. It is especially useful to reset all keywords before GetMAIN, etc.
5. Many math functions are moved to a new module mathfun, internally.
7. MAIN now accpets INF and NaN as the value of input parameters.
Fixed:
1. Imprementation of InverseErf to have caused seg fault, since the last version.
2. tsolqu had a large local varibles for spins.
3. A bung in LinearSolve (tsolvm), recently introduced.
Changes:
1. Flags NOINTRES and NOHALFRES now ignore unstableness with trace > 2 and trace < -2, respectively.
2. SADHelp.html starts to use LaTeX output for some texts with equations.
Added:
1. latex.n to provide library to include LaTeX output, now used in SADHelp.
Fixed:
1. A bug, probably very old, was found in OpenRead to read the results of a command (file name starting with "!"), when the result is null.
2. Compatibility of tfCrypt_.c (see Yamamoto-san's posts).
Changes:
1. The convergence algorithm for GO now always terminates at a certain number of iterations.