<!doctype html public "-//w3c//dtd html 4.0 transitional//en">

<html>

Hy Robert and neutron-mc users,

<p>that's right, I think the buffer line in Scilab is limited, as well

as the 'sack' for user variables. This stack is automatically inreased

by McStas, but there is indeed an intrinsic limit in Scilab for the size

of the data and length of lines in files.

<p>The only solution to this 'known bug' is to limit the size of monitors

toreasonable values.

Emmanuel.

<p>"Dalgliesh, RM (Robert)" wrote:

<blockquote TYPE=CITE>Has anyone else run into a problem with the length

of the data line that

scilab appears to be able to read in?

<br>I find that if you run with 1000 channels in a linear detector and

gather

<br>enough counts that scilab issues this error when you try to run the

mcstas.sci file.

<br>(1000 channels may seem a bit greedy but I am investigating a few instrument

<br>specific things that require fine wavelength resultion over a wide

range.)

<p> !--error     6

<br>inconsistent row/column dimensions

<br>at line      34 of function get_lmon1_dat           

called by :

<br>mc_lmon1_dat=get_lmon1_dat();

<br>line    67 of exec file called by :

<p>line     6 of function mcload_inline           

called by :

<br>line     3 of function mcplot_inline           

called by :

<br>line    60 of function get_mcstas              

called by :

<br>mc_mcstas=get_mcstas();

<br>line    97 of exec file called by :

<br>exec('mcstas.sci',-1);s=get_mcstas();

<p>if I run fewer neutrons through the simulation everthing is ok beacuse

0 is

put into the file as 0 not 0.00000e0.

running the same instrument and reading into matlab is not a problem.

<p>The simple test instrument below illustrates the problem. if you run

1e4

neutrons through scilab is happy but 1e5 and it won't like it.

<br>I'm going to have a look through the scilab documentation to see if

there is

<br>a limit to the length of line it can read in but just wondered if anyone

else had experienced the same problem.

<p>DEFINE INSTRUMENT test()

<br>DECLARE

<br>%{

<br>%}

<br>INITIALIZE

<br>%{

<br>%}

<br>TRACE

<br>COMPONENT Origin = Arm()

<br>  AT (0,0,0) ABSOLUTE

<p>COMPONENT source1 = Source_flat(

<br>    radius = 0.1, dist = 2.0, xw = 0.1, yh = 0.1, E0

= 100.0,

<br>    dE = 90.0)

<br>  AT (0.0, 0.0, 0.00001) RELATIVE Origin

<p>COMPONENT lmon1 = L_monitor(

<br>    nchan = 1000, filename = "lmon1.dat", xmin = -0.1,

<br>    xmax = 0.1, ymin = -0.1, ymax = 0.1, Lmin = 0.0,

Lmax = 10.0)

<br>  AT (0.0, 0.0, 10.1) RELATIVE Origin

<br>FINALLY

<br>%{

<br>%}

<br>/* The END token marks the instrument definition end */

<br>END

<p>_______________________________________________

<br>neutron-mc mailing list

<br>neutron-mc@risoe.dk

<br><a href="http://neutron.risoe.dk/mailman/listinfo/neutron-mc">http://neutron.risoe.dk/mailman/listinfo/neutron-mc</a></blockquote>

<pre>-- 

What's up Doc ?

--------------------------------------------

Emmanuel FARHI, <A HREF="http://www.ill.fr/Computing/people/Farhi">http://www.ill.fr/Computing/people/Farhi</A>  \|/ ____ \|/

CS-Group ILL4/156, Institut Laue-Langevin (ILL) Grenoble  ~@-/ oO \-@~

6 rue J. Horowitz, BP 156, 38042 Grenoble Cedex 9,France  /_( \__/ )_\

Work :Tel (33/0) 4 76 20 71 35. Fax (33/0) 4 76 20 76 48     \__U_/</pre>

 </html>