Skip to content
Snippets Groups Projects
Select Git revision
  • 3e69309153ece95eef6a5b197393b70eafa9eb82
  • master default protected
  • develop-GA
  • timeFstatmap
  • add-higher-spindown-components
  • develop-DK
  • adds-header-to-grid-search
  • v1.2
  • v1.1.2
  • v1.1.0
  • v1.0.1
11 results

semicoherent_glitch_robust_directed_MCMC_search_on_1_glitch.py

Blame
  • Forked from Gregory Ashton / PyFstat
    Source project has a limited visibility.
    FT_example3.m 1.75 KiB
    
    % Example file for using SimTools with Finesse
    
    clear all;
    
    % name of kat file which contains 'blocks'
    inname='testblock.kat';
    
    % read in block from testblock.kat
    block=FT_read_blocks_from_file(inname);
    
    myblock=FT_copy_block(block,'cavity');
    
    %FT_print_block(myblock,'cavity');
    
    disp('----------------------------------------------------------');
    disp('Printing reflectivities of all mirrors:');
    
    % Printing all reflectivities of the mirrors:
    for i=1:myblock.n_lines
      mir=FT_parse_line(myblock.txt_lines(i),[1,2,3]);
      command=cell2mat(mir{1});
      name=cell2mat(mir{2});
      if (command(1)=='m')
        disp(sprintf('mirror %s reflectance: %g',name,mir{3}));
      end
    end
     
    % now we change the one of m2
    [mir_lines]=FT_find_element_in_active_block(myblock,1,'m');
    
    for i=mir_lines
      [mir_name,p1,p2]=FT_get_element_from_line(myblock.txt_lines(i),2);
      if (mir_name=='m2')
        myblock.txt_lines(i)={FT_replace_value_in_line(myblock.txt_lines(i),3,0.7)};
      end
    end
    
    disp('----------------------------------------------------------');
    disp('Printing block after m2 R has been changed');
    
    FT_print_block(myblock,'cavity');
    
    disp('----------------------------------------------------------');
    disp('Finding all blocks containg mirrors');
    [b_index,l1,l2,lines]=FT_find_element_in_all_blocks(block,1,'m');
    for i=1:length(b_index)
      disp(sprintf('Found %d mirrors in block %s\n',l2(i)-l1(i)+1,block(b_index(i)).name));
    end
    
    disp('----------------------------------------------------------');
    disp('Finding all blocks containg the text m1');
    [b_index,l1,l2,lines]=FT_find_text_in_all_blocks(block,'m1');
    for i=1:length(b_index)
      disp(sprintf('in block %s:\n',block(b_index(i)).name));
      for j=l1(i):l2(i)
        disp(sprintf('line %d: %s:\n',lines(j),cell2mat(block(b_index(i)).txt_lines(lines(j)))));
      end
    end