#define NULL ((void *)0)
typedef unsigned long size_t; // Customize by platform.
typedef long intptr_t;
typedef unsigned long uintptr_t;
typedef long scalar_t__; // Either arithmetic or pointer type.
/* By default, we understand bool (as a convenience). */
typedef int bool;
#define false 0
#define true 1

/* Forward declarations */

/* Type definitions */

/* Variables and functions */
int /*<<< orphan*/ exit(int);
int /*<<< orphan*/ fprintf(int /*<<< orphan*/, char *);
int /*<<< orphan*/ stderr;
int usage() {
  fprintf(stderr, "usage: ms nsam howmany \n");
  fprintf(stderr, "  Options: \n");
  fprintf(stderr, "\t -t theta   (this option and/or the next must be used. "
                  "Theta = 4*N0*u )\n");
  fprintf(stderr, "\t -s segsites   ( fixed number of segregating sites)\n");
  fprintf(stderr, "\t -T          (Output gene tree.)\n");
  fprintf(stderr, "\t -F minfreq     Output only sites with freq of minor "
                  "allele >= minfreq.\n");
  fprintf(stderr, "\t -r rho nsites     (rho here is 4Nc)\n");
  fprintf(stderr, "\t\t -c f track_len   (f = ratio of conversion rate to rec "
                  "rate. tracklen is mean length.) \n");
  fprintf(
      stderr,
      "\t\t\t if rho = 0.,  f = 4*N0*g, with g the gene conversion rate.\n");
  fprintf(stderr,
          "\t -G alpha  ( N(t) = N0*exp(-alpha*t) .  alpha = -log(Np/Nr)/t\n");
  fprintf(stderr, "\t -I npop n1 n2 ... [mig_rate] (all elements of mig matrix "
                  "set to mig_rate/(npop-1) \n");
  fprintf(stderr,
          "\t\t -m i j m_ij    (i,j-th element of mig matrix set to m_ij.)\n");
  fprintf(stderr, "\t\t -ma m_11 m_12 m_13 m_21 m_22 m_23 ...(Assign values to "
                  "elements of migration matrix.)\n");
  fprintf(stderr, "\t\t -n i size_i   (popi has size set to size_i*N0 \n");
  fprintf(stderr,
          "\t\t -g i alpha_i  (If used must appear after -M option.)\n");
  fprintf(stderr, "\t   The following options modify parameters at the time "
                  "'t' specified as the first argument:\n");
  fprintf(stderr, "\t -eG t alpha  (Modify growth rate of all pop's.)\n");
  fprintf(stderr, "\t -eg t i alpha_i  (Modify growth rate of pop i.) \n");
  fprintf(stderr, "\t -eM t mig_rate   (Modify the mig matrix so all elements "
                  "are mig_rate/(npop-1)\n");
  fprintf(stderr, "\t -em t i j m_ij    (i,j-th element of mig matrix set to "
                  "m_ij at time t )\n");
  fprintf(stderr, "\t -ema t npop  m_11 m_12 m_13 m_21 m_22 m_23 ...(Assign "
                  "values to elements of migration matrix.)\n");
  fprintf(stderr, "\t -eN t size  (Modify pop sizes. New sizes = size*N0 ) \n");
  fprintf(stderr, "\t -en t i size_i  (Modify pop size of pop i.  New size of "
                  "popi = size_i*N0 .)\n");
  fprintf(stderr, "\t -es t i proportion  (Split: pop i -> pop-i + pop-npop, "
                  "npop increases by 1.\n");
  fprintf(stderr, "\t\t proportion is probability that each lineage stays in "
                  "pop-i. (p, 1-p are admixt. proport.\n");
  fprintf(stderr, "\t\t Size of pop npop is set to N0 and alpha = 0.0 , size "
                  "and alpha of pop i are unchanged.\n");
  fprintf(stderr,
          "\t -ej t i j   ( Join lineages in pop i and pop j into pop j\n");
  fprintf(stderr, "\t\t  size, alpha and M are unchanged.\n");
  fprintf(stderr, "\t -v J a_1 b_1 I_1 ... a_J b_J I_J  ( Incorporate J "
                  "crossover hotspots, each of intensity I_j on [a_j,b_j] )\n");
  fprintf(stderr,
          "\t -V K a_1 b_1 I_1 ... a_K b_K I_K  ( Incorporate K gene "
          "conversion hotspots, each of intensity I_k on [a_k,b_k] )\n");
  fprintf(stderr, "\t  -f filename     ( Read command line arguments from file "
                  "filename.)\n");
  fprintf(stderr, " See msdoc.pdf for explanation of these parameters.\n");
  exit(1);
}


x86 assembly (clang, -O0):
<usage>:
push   %rbp
mov    %rsp,%rbp
mov    0x2eb5(%rip),%rax
mov    (%rax),%edi
lea    0xecc(%rip),%rsi
call   1030 <fprintf@plt>
mov    0x2ea0(%rip),%rax
mov    (%rax),%edi
lea    0xed0(%rip),%rsi
call   1030 <fprintf@plt>
mov    0x2e8b(%rip),%rax
mov    (%rax),%edi
lea    0xec8(%rip),%rsi
call   1030 <fprintf@plt>
mov    0x2e76(%rip),%rax
mov    (%rax),%edi
lea    0xefd(%rip),%rsi
call   1030 <fprintf@plt>
mov    0x2e61(%rip),%rax
mov    (%rax),%edi
lea    0xf1e(%rip),%rsi
call   1030 <fprintf@plt>
mov    0x2e4c(%rip),%rax
mov    (%rax),%edi
lea    0xf2c(%rip),%rsi
call   1030 <fprintf@plt>
mov    0x2e37(%rip),%rax
mov    (%rax),%edi
lea    0xf61(%rip),%rsi
call   1030 <fprintf@plt>
mov    0x2e22(%rip),%rax
mov    (%rax),%edi
lea    0xf73(%rip),%rsi
call   1030 <fprintf@plt>
mov    0x2e0d(%rip),%rax
mov    (%rax),%edi
lea    0xfb9(%rip),%rsi
call   1030 <fprintf@plt>
mov    0x2df8(%rip),%rax
mov    (%rax),%edi
lea    0xfe4(%rip),%rsi
call   1030 <fprintf@plt>
mov    0x2de3(%rip),%rax
mov    (%rax),%edi
lea    0x100f(%rip),%rsi
call   1030 <fprintf@plt>
mov    0x2dce(%rip),%rax
mov    (%rax),%edi
lea    0x1050(%rip),%rsi
call   1030 <fprintf@plt>
mov    0x2db9(%rip),%rax
mov    (%rax),%edi
lea    0x107a(%rip),%rsi
call   1030 <fprintf@plt>
mov    0x2da4(%rip),%rax
mov    (%rax),%edi
lea    0x10bf(%rip),%rsi
call   1030 <fprintf@plt>
mov    0x2d8f(%rip),%rax
mov    (%rax),%edi
lea    0x10dd(%rip),%rsi
call   1030 <fprintf@plt>
mov    0x2d7a(%rip),%rax
mov    (%rax),%edi
lea    0x1101(%rip),%rsi
call   1030 <fprintf@plt>
mov    0x2d65(%rip),%rax
mov    (%rax),%edi
lea    0x114a(%rip),%rsi
call   1030 <fprintf@plt>
mov    0x2d50(%rip),%rax
mov    (%rax),%edi
lea    0x1168(%rip),%rsi
call   1030 <fprintf@plt>
mov    0x2d3b(%rip),%rax
mov    (%rax),%edi
lea    0x1187(%rip),%rsi
call   1030 <fprintf@plt>
mov    0x2d26(%rip),%rax
mov    (%rax),%edi
lea    0x11c3(%rip),%rsi
call   1030 <fprintf@plt>
mov    0x2d11(%rip),%rax
mov    (%rax),%edi
lea    0x11f9(%rip),%rsi
call   1030 <fprintf@plt>
mov    0x2cfc(%rip),%rax
mov    (%rax),%edi
lea    0x1246(%rip),%rsi
call   1030 <fprintf@plt>
mov    0x2ce7(%rip),%rax
mov    (%rax),%edi
lea    0x126a(%rip),%rsi
call   1030 <fprintf@plt>
mov    0x2cd2(%rip),%rax
mov    (%rax),%edi
lea    0x12a4(%rip),%rsi
call   1030 <fprintf@plt>
mov    0x2cbd(%rip),%rax
mov    (%rax),%edi
lea    0x12de(%rip),%rsi
call   1030 <fprintf@plt>
mov    0x2ca8(%rip),%rax
mov    (%rax),%edi
lea    0x1326(%rip),%rsi
call   1030 <fprintf@plt>
mov    0x2c93(%rip),%rax
mov    (%rax),%edi
lea    0x136c(%rip),%rsi
call   1030 <fprintf@plt>
mov    0x2c7e(%rip),%rax
mov    (%rax),%edi
lea    0x1394(%rip),%rsi
call   1030 <fprintf@plt>
mov    0x2c69(%rip),%rax
mov    (%rax),%edi
lea    0x13a5(%rip),%rsi
call   1030 <fprintf@plt>
mov    0x2c54(%rip),%rax
mov    (%rax),%edi
lea    0x13fe(%rip),%rsi
call   1030 <fprintf@plt>
mov    0x2c3f(%rip),%rax
mov    (%rax),%edi
lea    0x145d(%rip),%rsi
call   1030 <fprintf@plt>
mov    0x2c2a(%rip),%rax
mov    (%rax),%edi
lea    0x148f(%rip),%rsi
call   1030 <fprintf@plt>
mov    $0x1,%edi
call   1040 <exit@plt>

x86 assembly (clang, -O2):
<usage>:
push   %rbx
mov    0x2eb8(%rip),%rbx
mov    (%rbx),%edi
lea    0xecf(%rip),%rsi
call   1030 <fprintf@plt>
mov    (%rbx),%edi
lea    0xeda(%rip),%rsi
call   1030 <fprintf@plt>
mov    (%rbx),%edi
lea    0xed9(%rip),%rsi
call   1030 <fprintf@plt>
mov    (%rbx),%edi
lea    0xf15(%rip),%rsi
call   1030 <fprintf@plt>
mov    (%rbx),%edi
lea    0xf3d(%rip),%rsi
call   1030 <fprintf@plt>
mov    (%rbx),%edi
lea    0xf52(%rip),%rsi
call   1030 <fprintf@plt>
mov    (%rbx),%edi
lea    0xf8e(%rip),%rsi
call   1030 <fprintf@plt>
mov    (%rbx),%edi
lea    0xfa7(%rip),%rsi
call   1030 <fprintf@plt>
mov    (%rbx),%edi
lea    0xff4(%rip),%rsi
call   1030 <fprintf@plt>
mov    (%rbx),%edi
lea    0x1026(%rip),%rsi
call   1030 <fprintf@plt>
mov    (%rbx),%edi
lea    0x1058(%rip),%rsi
call   1030 <fprintf@plt>
mov    (%rbx),%edi
lea    0x10a0(%rip),%rsi
call   1030 <fprintf@plt>
mov    (%rbx),%edi
lea    0x10d1(%rip),%rsi
call   1030 <fprintf@plt>
mov    (%rbx),%edi
lea    0x111d(%rip),%rsi
call   1030 <fprintf@plt>
mov    (%rbx),%edi
lea    0x1142(%rip),%rsi
call   1030 <fprintf@plt>
mov    (%rbx),%edi
lea    0x116d(%rip),%rsi
call   1030 <fprintf@plt>
mov    (%rbx),%edi
lea    0x11bd(%rip),%rsi
call   1030 <fprintf@plt>
mov    (%rbx),%edi
lea    0x11e2(%rip),%rsi
call   1030 <fprintf@plt>
mov    (%rbx),%edi
lea    0x1208(%rip),%rsi
call   1030 <fprintf@plt>
mov    (%rbx),%edi
lea    0x124b(%rip),%rsi
call   1030 <fprintf@plt>
mov    (%rbx),%edi
lea    0x1288(%rip),%rsi
call   1030 <fprintf@plt>
mov    (%rbx),%edi
lea    0x12dc(%rip),%rsi
call   1030 <fprintf@plt>
mov    (%rbx),%edi
lea    0x1307(%rip),%rsi
call   1030 <fprintf@plt>
mov    (%rbx),%edi
lea    0x1348(%rip),%rsi
call   1030 <fprintf@plt>
mov    (%rbx),%edi
lea    0x1389(%rip),%rsi
call   1030 <fprintf@plt>
mov    (%rbx),%edi
lea    0x13d8(%rip),%rsi
call   1030 <fprintf@plt>
mov    (%rbx),%edi
lea    0x1425(%rip),%rsi
call   1030 <fprintf@plt>
mov    (%rbx),%edi
lea    0x1454(%rip),%rsi
call   1030 <fprintf@plt>
mov    (%rbx),%edi
lea    0x146c(%rip),%rsi
call   1030 <fprintf@plt>
mov    (%rbx),%edi
lea    0x14cc(%rip),%rsi
call   1030 <fprintf@plt>
mov    (%rbx),%edi
lea    0x1532(%rip),%rsi
call   1030 <fprintf@plt>
mov    (%rbx),%edi
lea    0x156b(%rip),%rsi
call   1030 <fprintf@plt>
mov    $0x1,%edi
call   1040 <exit@plt>

LLVM IR (-O0):
@stderr = global i32 0
@.str = [25 x i8] c"usage: ms nsam howmany \0A\00"
@.str.1 = [13 x i8] c" Options: \0A\00"
@.str.2 = [74 x i8] c"\09 -t theta (this option and/or the next must be used. Theta = 4*N0*u )\0A\00"
@.str.3 = [54 x i8] c"\09 -s segsites ( fixed number of segregating sites)\0A\00"
@.str.4 = [35 x i8] c"\09 -T (Output gene tree.)\0A\00"
@.str.5 = [74 x i8] c"\09 -F minfreq Output only sites with freq of minor allele >= minfreq.\0A\00"
@.str.6 = [39 x i8] c"\09 -r rho nsites (rho here is 4Nc)\0A\00"
@.str.7 = [91 x i8] c"\09\09 -c f track_len (f = ratio of conversion rate to rec rate. tracklen is mean length.) \0A\00"
@.str.8 = [64 x i8] c"\09\09\09 if rho = 0., f = 4*N0*g, with g the gene conversion rate.\0A\00"
@.str.9 = [64 x i8] c"\09 -G alpha ( N(t) = N0*exp(-alpha*t) . alpha = -log(Np/Nr)/t\0A\00"
@.str.10 = [86 x i8] c"\09 -I npop n1 n2 ... [mig_rate] (all elements of mig matrix set to mig_rate/(npop-1) \0A\00"
@.str.11 = [63 x i8] c"\09\09 -m i j m_ij (i,j-th element of mig matrix set to m_ij.)\0A\00"
@.str.12 = [90 x i8] c"\09\09 -ma m_11 m_12 m_13 m_21 m_22 m_23 ...(Assign values to elements of migration matrix.)\0A\00"
@.str.13 = [51 x i8] c"\09\09 -n i size_i (popi has size set to size_i*N0 \0A\00"
@.str.14 = [57 x i8] c"\09\09 -g i alpha_i (If used must appear after -M option.)\0A\00"
@.str.15 = [94 x i8] c"\09 The following options modify parameters at the time 't' specified as the first argument:\0A\00"
@.str.16 = [51 x i8] c"\09 -eG t alpha (Modify growth rate of all pop's.)\0A\00"
@.str.17 = [52 x i8] c"\09 -eg t i alpha_i (Modify growth rate of pop i.) \0A\00"
@.str.18 = [81 x i8] c"\09 -eM t mig_rate (Modify the mig matrix so all elements are mig_rate/(npop-1)\0A\00"
@.str.19 = [75 x i8] c"\09 -em t i j m_ij (i,j-th element of mig matrix set to m_ij at time t )\0A\00"
@.str.20 = [98 x i8] c"\09 -ema t npop m_11 m_12 m_13 m_21 m_22 m_23 ...(Assign values to elements of migration matrix.)\0A\00"
@.str.21 = [57 x i8] c"\09 -eN t size (Modify pop sizes. New sizes = size*N0 ) \0A\00"
@.str.22 = [79 x i8] c"\09 -en t i size_i (Modify pop size of pop i. New size of popi = size_i*N0 .)\0A\00"
@.str.23 = [79 x i8] c"\09 -es t i proportion (Split: pop i -> pop-i + pop-npop, npop increases by 1.\0A\00"
@.str.24 = [93 x i8] c"\09\09 proportion is probability that each lineage stays in pop-i. (p, 1-p are admixt. proport.\0A\00"
@.str.25 = [91 x i8] c"\09\09 Size of pop npop is set to N0 and alpha = 0.0 , size and alpha of pop i are unchanged.\0A\00"
@.str.26 = [61 x i8] c"\09 -ej t i j ( Join lineages in pop i and pop j into pop j\0A\00"
@.str.27 = [38 x i8] c"\09\09 size, alpha and M are unchanged.\0A\00"
@.str.28 = [110 x i8] c"\09 -v J a_1 b_1 I_1 ... a_J b_J I_J ( Incorporate J crossover hotspots, each of intensity I_j on [a_j,b_j] )\0A\00"
@.str.29 = [116 x i8] c"\09 -V K a_1 b_1 I_1 ... a_K b_K I_K ( Incorporate K gene conversion hotspots, each of intensity I_k on [a_k,b_k] )\0A\00"
@.str.30 = [71 x i8] c"\09 -f filename ( Read command line arguments from file filename.)\0A\00"
@.str.31 = [53 x i8] c" See msdoc.pdf for explanation of these parameters.\0A\00"
define i32 @usage() {
%1 = load i32, i32* @stderr
%2 = call i32 @fprintf(i32 %1, i8* getelementptr inbounds ([25 x i8], [25 x i8]* @.str, i64 0, i64 0))
%3 = load i32, i32* @stderr
%4 = call i32 @fprintf(i32 %3, i8* getelementptr inbounds ([13 x i8], [13 x i8]* @.str.1, i64 0, i64 0))
%5 = load i32, i32* @stderr
%6 = call i32 @fprintf(i32 %5, i8* getelementptr inbounds ([74 x i8], [74 x i8]* @.str.2, i64 0, i64 0))
%7 = load i32, i32* @stderr
%8 = call i32 @fprintf(i32 %7, i8* getelementptr inbounds ([54 x i8], [54 x i8]* @.str.3, i64 0, i64 0))
%9 = load i32, i32* @stderr
%10 = call i32 @fprintf(i32 %9, i8* getelementptr inbounds ([35 x i8], [35 x i8]* @.str.4, i64 0, i64 0))
%11 = load i32, i32* @stderr
%12 = call i32 @fprintf(i32 %11, i8* getelementptr inbounds ([74 x i8], [74 x i8]* @.str.5, i64 0, i64 0))
%13 = load i32, i32* @stderr
%14 = call i32 @fprintf(i32 %13, i8* getelementptr inbounds ([39 x i8], [39 x i8]* @.str.6, i64 0, i64 0))
%15 = load i32, i32* @stderr
%16 = call i32 @fprintf(i32 %15, i8* getelementptr inbounds ([91 x i8], [91 x i8]* @.str.7, i64 0, i64 0))
%17 = load i32, i32* @stderr
%18 = call i32 @fprintf(i32 %17, i8* getelementptr inbounds ([64 x i8], [64 x i8]* @.str.8, i64 0, i64 0))
%19 = load i32, i32* @stderr
%20 = call i32 @fprintf(i32 %19, i8* getelementptr inbounds ([64 x i8], [64 x i8]* @.str.9, i64 0, i64 0))
%21 = load i32, i32* @stderr
%22 = call i32 @fprintf(i32 %21, i8* getelementptr inbounds ([86 x i8], [86 x i8]* @.str.10, i64 0, i64 0))
%23 = load i32, i32* @stderr
%24 = call i32 @fprintf(i32 %23, i8* getelementptr inbounds ([63 x i8], [63 x i8]* @.str.11, i64 0, i64 0))
%25 = load i32, i32* @stderr
%26 = call i32 @fprintf(i32 %25, i8* getelementptr inbounds ([90 x i8], [90 x i8]* @.str.12, i64 0, i64 0))
%27 = load i32, i32* @stderr
%28 = call i32 @fprintf(i32 %27, i8* getelementptr inbounds ([51 x i8], [51 x i8]* @.str.13, i64 0, i64 0))
%29 = load i32, i32* @stderr
%30 = call i32 @fprintf(i32 %29, i8* getelementptr inbounds ([57 x i8], [57 x i8]* @.str.14, i64 0, i64 0))
%31 = load i32, i32* @stderr
%32 = call i32 @fprintf(i32 %31, i8* getelementptr inbounds ([94 x i8], [94 x i8]* @.str.15, i64 0, i64 0))
%33 = load i32, i32* @stderr
%34 = call i32 @fprintf(i32 %33, i8* getelementptr inbounds ([51 x i8], [51 x i8]* @.str.16, i64 0, i64 0))
%35 = load i32, i32* @stderr
%36 = call i32 @fprintf(i32 %35, i8* getelementptr inbounds ([52 x i8], [52 x i8]* @.str.17, i64 0, i64 0))
%37 = load i32, i32* @stderr
%38 = call i32 @fprintf(i32 %37, i8* getelementptr inbounds ([81 x i8], [81 x i8]* @.str.18, i64 0, i64 0))
%39 = load i32, i32* @stderr
%40 = call i32 @fprintf(i32 %39, i8* getelementptr inbounds ([75 x i8], [75 x i8]* @.str.19, i64 0, i64 0))
%41 = load i32, i32* @stderr
%42 = call i32 @fprintf(i32 %41, i8* getelementptr inbounds ([98 x i8], [98 x i8]* @.str.20, i64 0, i64 0))
%43 = load i32, i32* @stderr
%44 = call i32 @fprintf(i32 %43, i8* getelementptr inbounds ([57 x i8], [57 x i8]* @.str.21, i64 0, i64 0))
%45 = load i32, i32* @stderr
%46 = call i32 @fprintf(i32 %45, i8* getelementptr inbounds ([79 x i8], [79 x i8]* @.str.22, i64 0, i64 0))
%47 = load i32, i32* @stderr
%48 = call i32 @fprintf(i32 %47, i8* getelementptr inbounds ([79 x i8], [79 x i8]* @.str.23, i64 0, i64 0))
%49 = load i32, i32* @stderr
%50 = call i32 @fprintf(i32 %49, i8* getelementptr inbounds ([93 x i8], [93 x i8]* @.str.24, i64 0, i64 0))
%51 = load i32, i32* @stderr
%52 = call i32 @fprintf(i32 %51, i8* getelementptr inbounds ([91 x i8], [91 x i8]* @.str.25, i64 0, i64 0))
%53 = load i32, i32* @stderr
%54 = call i32 @fprintf(i32 %53, i8* getelementptr inbounds ([61 x i8], [61 x i8]* @.str.26, i64 0, i64 0))
%55 = load i32, i32* @stderr
%56 = call i32 @fprintf(i32 %55, i8* getelementptr inbounds ([38 x i8], [38 x i8]* @.str.27, i64 0, i64 0))
%57 = load i32, i32* @stderr
%58 = call i32 @fprintf(i32 %57, i8* getelementptr inbounds ([110 x i8], [110 x i8]* @.str.28, i64 0, i64 0))
%59 = load i32, i32* @stderr
%60 = call i32 @fprintf(i32 %59, i8* getelementptr inbounds ([116 x i8], [116 x i8]* @.str.29, i64 0, i64 0))
%61 = load i32, i32* @stderr
%62 = call i32 @fprintf(i32 %61, i8* getelementptr inbounds ([71 x i8], [71 x i8]* @.str.30, i64 0, i64 0))
%63 = load i32, i32* @stderr
%64 = call i32 @fprintf(i32 %63, i8* getelementptr inbounds ([53 x i8], [53 x i8]* @.str.31, i64 0, i64 0))
%65 = call i32 @exit(i32 1)
unreachable
}

LLVM IR (-O2):
@stderr = global i32 0
@.str = [25 x i8] c"usage: ms nsam howmany \0A\00"
@.str.1 = [13 x i8] c" Options: \0A\00"
@.str.2 = [74 x i8] c"\09 -t theta (this option and/or the next must be used. Theta = 4*N0*u )\0A\00"
@.str.3 = [54 x i8] c"\09 -s segsites ( fixed number of segregating sites)\0A\00"
@.str.4 = [35 x i8] c"\09 -T (Output gene tree.)\0A\00"
@.str.5 = [74 x i8] c"\09 -F minfreq Output only sites with freq of minor allele >= minfreq.\0A\00"
@.str.6 = [39 x i8] c"\09 -r rho nsites (rho here is 4Nc)\0A\00"
@.str.7 = [91 x i8] c"\09\09 -c f track_len (f = ratio of conversion rate to rec rate. tracklen is mean length.) \0A\00"
@.str.8 = [64 x i8] c"\09\09\09 if rho = 0., f = 4*N0*g, with g the gene conversion rate.\0A\00"
@.str.9 = [64 x i8] c"\09 -G alpha ( N(t) = N0*exp(-alpha*t) . alpha = -log(Np/Nr)/t\0A\00"
@.str.10 = [86 x i8] c"\09 -I npop n1 n2 ... [mig_rate] (all elements of mig matrix set to mig_rate/(npop-1) \0A\00"
@.str.11 = [63 x i8] c"\09\09 -m i j m_ij (i,j-th element of mig matrix set to m_ij.)\0A\00"
@.str.12 = [90 x i8] c"\09\09 -ma m_11 m_12 m_13 m_21 m_22 m_23 ...(Assign values to elements of migration matrix.)\0A\00"
@.str.13 = [51 x i8] c"\09\09 -n i size_i (popi has size set to size_i*N0 \0A\00"
@.str.14 = [57 x i8] c"\09\09 -g i alpha_i (If used must appear after -M option.)\0A\00"
@.str.15 = [94 x i8] c"\09 The following options modify parameters at the time 't' specified as the first argument:\0A\00"
@.str.16 = [51 x i8] c"\09 -eG t alpha (Modify growth rate of all pop's.)\0A\00"
@.str.17 = [52 x i8] c"\09 -eg t i alpha_i (Modify growth rate of pop i.) \0A\00"
@.str.18 = [81 x i8] c"\09 -eM t mig_rate (Modify the mig matrix so all elements are mig_rate/(npop-1)\0A\00"
@.str.19 = [75 x i8] c"\09 -em t i j m_ij (i,j-th element of mig matrix set to m_ij at time t )\0A\00"
@.str.20 = [98 x i8] c"\09 -ema t npop m_11 m_12 m_13 m_21 m_22 m_23 ...(Assign values to elements of migration matrix.)\0A\00"
@.str.21 = [57 x i8] c"\09 -eN t size (Modify pop sizes. New sizes = size*N0 ) \0A\00"
@.str.22 = [79 x i8] c"\09 -en t i size_i (Modify pop size of pop i. New size of popi = size_i*N0 .)\0A\00"
@.str.23 = [79 x i8] c"\09 -es t i proportion (Split: pop i -> pop-i + pop-npop, npop increases by 1.\0A\00"
@.str.24 = [93 x i8] c"\09\09 proportion is probability that each lineage stays in pop-i. (p, 1-p are admixt. proport.\0A\00"
@.str.25 = [91 x i8] c"\09\09 Size of pop npop is set to N0 and alpha = 0.0 , size and alpha of pop i are unchanged.\0A\00"
@.str.26 = [61 x i8] c"\09 -ej t i j ( Join lineages in pop i and pop j into pop j\0A\00"
@.str.27 = [38 x i8] c"\09\09 size, alpha and M are unchanged.\0A\00"
@.str.28 = [110 x i8] c"\09 -v J a_1 b_1 I_1 ... a_J b_J I_J ( Incorporate J crossover hotspots, each of intensity I_j on [a_j,b_j] )\0A\00"
@.str.29 = [116 x i8] c"\09 -V K a_1 b_1 I_1 ... a_K b_K I_K ( Incorporate K gene conversion hotspots, each of intensity I_k on [a_k,b_k] )\0A\00"
@.str.30 = [71 x i8] c"\09 -f filename ( Read command line arguments from file filename.)\0A\00"
@.str.31 = [53 x i8] c" See msdoc.pdf for explanation of these parameters.\0A\00"
define i32 @usage() {
%1 = load i32, i32* @stderr
%2 = tail call i32 @fprintf(i32 %1, i8* getelementptr inbounds ([25 x i8], [25 x i8]* @.str, i64 0, i64 0))
%3 = load i32, i32* @stderr
%4 = tail call i32 @fprintf(i32 %3, i8* getelementptr inbounds ([13 x i8], [13 x i8]* @.str.1, i64 0, i64 0))
%5 = load i32, i32* @stderr
%6 = tail call i32 @fprintf(i32 %5, i8* getelementptr inbounds ([74 x i8], [74 x i8]* @.str.2, i64 0, i64 0))
%7 = load i32, i32* @stderr
%8 = tail call i32 @fprintf(i32 %7, i8* getelementptr inbounds ([54 x i8], [54 x i8]* @.str.3, i64 0, i64 0))
%9 = load i32, i32* @stderr
%10 = tail call i32 @fprintf(i32 %9, i8* getelementptr inbounds ([35 x i8], [35 x i8]* @.str.4, i64 0, i64 0))
%11 = load i32, i32* @stderr
%12 = tail call i32 @fprintf(i32 %11, i8* getelementptr inbounds ([74 x i8], [74 x i8]* @.str.5, i64 0, i64 0))
%13 = load i32, i32* @stderr
%14 = tail call i32 @fprintf(i32 %13, i8* getelementptr inbounds ([39 x i8], [39 x i8]* @.str.6, i64 0, i64 0))
%15 = load i32, i32* @stderr
%16 = tail call i32 @fprintf(i32 %15, i8* getelementptr inbounds ([91 x i8], [91 x i8]* @.str.7, i64 0, i64 0))
%17 = load i32, i32* @stderr
%18 = tail call i32 @fprintf(i32 %17, i8* getelementptr inbounds ([64 x i8], [64 x i8]* @.str.8, i64 0, i64 0))
%19 = load i32, i32* @stderr
%20 = tail call i32 @fprintf(i32 %19, i8* getelementptr inbounds ([64 x i8], [64 x i8]* @.str.9, i64 0, i64 0))
%21 = load i32, i32* @stderr
%22 = tail call i32 @fprintf(i32 %21, i8* getelementptr inbounds ([86 x i8], [86 x i8]* @.str.10, i64 0, i64 0))
%23 = load i32, i32* @stderr
%24 = tail call i32 @fprintf(i32 %23, i8* getelementptr inbounds ([63 x i8], [63 x i8]* @.str.11, i64 0, i64 0))
%25 = load i32, i32* @stderr
%26 = tail call i32 @fprintf(i32 %25, i8* getelementptr inbounds ([90 x i8], [90 x i8]* @.str.12, i64 0, i64 0))
%27 = load i32, i32* @stderr
%28 = tail call i32 @fprintf(i32 %27, i8* getelementptr inbounds ([51 x i8], [51 x i8]* @.str.13, i64 0, i64 0))
%29 = load i32, i32* @stderr
%30 = tail call i32 @fprintf(i32 %29, i8* getelementptr inbounds ([57 x i8], [57 x i8]* @.str.14, i64 0, i64 0))
%31 = load i32, i32* @stderr
%32 = tail call i32 @fprintf(i32 %31, i8* getelementptr inbounds ([94 x i8], [94 x i8]* @.str.15, i64 0, i64 0))
%33 = load i32, i32* @stderr
%34 = tail call i32 @fprintf(i32 %33, i8* getelementptr inbounds ([51 x i8], [51 x i8]* @.str.16, i64 0, i64 0))
%35 = load i32, i32* @stderr
%36 = tail call i32 @fprintf(i32 %35, i8* getelementptr inbounds ([52 x i8], [52 x i8]* @.str.17, i64 0, i64 0))
%37 = load i32, i32* @stderr
%38 = tail call i32 @fprintf(i32 %37, i8* getelementptr inbounds ([81 x i8], [81 x i8]* @.str.18, i64 0, i64 0))
%39 = load i32, i32* @stderr
%40 = tail call i32 @fprintf(i32 %39, i8* getelementptr inbounds ([75 x i8], [75 x i8]* @.str.19, i64 0, i64 0))
%41 = load i32, i32* @stderr
%42 = tail call i32 @fprintf(i32 %41, i8* getelementptr inbounds ([98 x i8], [98 x i8]* @.str.20, i64 0, i64 0))
%43 = load i32, i32* @stderr
%44 = tail call i32 @fprintf(i32 %43, i8* getelementptr inbounds ([57 x i8], [57 x i8]* @.str.21, i64 0, i64 0))
%45 = load i32, i32* @stderr
%46 = tail call i32 @fprintf(i32 %45, i8* getelementptr inbounds ([79 x i8], [79 x i8]* @.str.22, i64 0, i64 0))
%47 = load i32, i32* @stderr
%48 = tail call i32 @fprintf(i32 %47, i8* getelementptr inbounds ([79 x i8], [79 x i8]* @.str.23, i64 0, i64 0))
%49 = load i32, i32* @stderr
%50 = tail call i32 @fprintf(i32 %49, i8* getelementptr inbounds ([93 x i8], [93 x i8]* @.str.24, i64 0, i64 0))
%51 = load i32, i32* @stderr
%52 = tail call i32 @fprintf(i32 %51, i8* getelementptr inbounds ([91 x i8], [91 x i8]* @.str.25, i64 0, i64 0))
%53 = load i32, i32* @stderr
%54 = tail call i32 @fprintf(i32 %53, i8* getelementptr inbounds ([61 x i8], [61 x i8]* @.str.26, i64 0, i64 0))
%55 = load i32, i32* @stderr
%56 = tail call i32 @fprintf(i32 %55, i8* getelementptr inbounds ([38 x i8], [38 x i8]* @.str.27, i64 0, i64 0))
%57 = load i32, i32* @stderr
%58 = tail call i32 @fprintf(i32 %57, i8* getelementptr inbounds ([110 x i8], [110 x i8]* @.str.28, i64 0, i64 0))
%59 = load i32, i32* @stderr
%60 = tail call i32 @fprintf(i32 %59, i8* getelementptr inbounds ([116 x i8], [116 x i8]* @.str.29, i64 0, i64 0))
%61 = load i32, i32* @stderr
%62 = tail call i32 @fprintf(i32 %61, i8* getelementptr inbounds ([71 x i8], [71 x i8]* @.str.30, i64 0, i64 0))
%63 = load i32, i32* @stderr
%64 = tail call i32 @fprintf(i32 %63, i8* getelementptr inbounds ([53 x i8], [53 x i8]* @.str.31, i64 0, i64 0))
%65 = tail call i32 @exit(i32 1)
unreachable
}