write_toast_qm

0001 function write_toast_qm(fname, fmdl)
0002 fid = fopen(fname,'w');
0003 dims = size(fmdl.nodes,2);
0004 elecs= length(fmdl.electrode);
0005 elecpos= zeros(elecs,dims);
0006 for i=1:elecs
0007     elnodes = fmdl.electrode(i).nodes;
0008     elecpos(i,:) = mean( fmdl.nodes(elnodes,:),1);
0009 end
0010 
0011 fprintf(fid,'QM file %dD\n',dims);
0012 fprintf(fid,'Dimension %d\n',dims);
0013 fprintf(fid,'SourceList %d fixed\n',elecs);
0014 fprintf(fid,'%6.4g     %6.4g\n', elecpos');
0015 fprintf(fid,'MeasurementList %d\n',elecs);
0016 fprintf(fid,'%6.4g     %6.4g\n', elecpos');
0017 fprintf(fid,'LinkList\n');
0018 for i=1:elecs
0019   fprintf(fid,'%d:',elecs);
0020   fprintf(fid,'%d ',0:elecs-1);
0021   fprintf(fid,'\n');
0022 end
0023 
0024 fclose(fid);