F64 SpanTime()
{
  if (run_bttn.state)
    return a.elapsed_t+tS-a.start_wall_t;
  else
    return a.elapsed_t;
}

F64 Cost(CMathODE *ode)
{
  MyMass   *tmpm;
  MySpring *tmps;
  F64 res=0;
  tmpm=ode->next_mass;
  while (tmpm!=&ode->next_mass)
  {
    res+=tmpm->cost;
    tmpm=tmpm->next;
  }
  tmps=ode->next_spring;
  while (tmps!=&ode->next_spring)
  {
    res+=tmps->cost;
    tmps=tmps->next;
  }
  return res;
}

U0 DrawIt(CTask *task,CDC *dc)
{
  MyMass   *tmpm;
  MySpring *tmps;
  tmps=ode->next_spring;
  while (tmps!=&ode->next_spring)
  {
    if (!(tmps->flags&SSF_INACTIVE))
    {
      dc->color=tmps->color;
      dc->thick=tmps->thick;
      GrLine3(dc,tmps->end1->x,tmps->end1->y,0,
              tmps->end2->x,tmps->end2->y,0);
    }
    tmps=tmps->next;
  }
  if (cursor_mass)
  {
    dc->color=RED;
    dc->thick=2;
    GrLine3(dc,ms.pos.x-task->pix_left-task->scroll_x,
            ms.pos.y-task->pix_top-task->scroll_y,0,
            cursor_mass->x,cursor_mass->y,0);
  }
  tmpm=ode->next_mass;
  while (tmpm!=&ode->next_mass)
  {
    if (!(tmpm->flags&MSF_INACTIVE))
    {
      dc->color=BLACK;
      GrCircle(dc,tmpm->x,tmpm->y,tmpm->radius);
      GrFloodFill(dc,tmpm->x,tmpm->y,TRUE);
      dc->color=tmpm->color;
      GrCircle(dc,tmpm->x,tmpm->y,tmpm->radius);
      GrFloodFill(dc,tmpm->x,tmpm->y,TRUE);
      dc->color=BLACK;
      GrCircle(dc,tmpm->x,tmpm->y,tmpm->radius);
    }
    tmpm=tmpm->next;
  }
  dc->color=BLACK;
  GrPrint(dc,90,0,"Cost:%12.2,f",Cost(ode));
  GrPrint(dc,90,FONT_HEIGHT,"Time:%12.2f",SpanTime);
}

MyMass *PlaceMass(I64 x, I64 y)
{
  MyMass *tmpm=CAlloc(sizeof(MyMass));
  tmpm->drag_profile_factor=1.0;
  tmpm->x=x;
  tmpm->y=y;
  tmpm->mass=MASS_MASS;
  tmpm->radius=MASS_RADIUS;
  tmpm->cost=25.0*COST_SCALE;
  tmpm->color=YELLOW;
  QueIns(tmpm,ode->last_mass);
  return tmpm;
}

U0 NullSpring(MySpring *tmps,F64 scale)
{
  F64 d=D3Dist(&tmps->end1->x,&tmps->end2->x);
  tmps->rest_len=d*scale;
  tmps->compression_strength=
  tmps->base_compression_strength/(tmps->rest_len+1.0);
  tmps->tensile_strength=tmps->base_tensile_strength/(tmps->rest_len+1.0);
  tmps->const=tmps->base_const/(tmps->rest_len+1.0);
  tmps->cost=tmps->base_cost*tmps->rest_len;
}

U0 MoveMass(MyMass *tmpm,I64 x, I64 y)
{
  MySpring *tmps;
  tmpm->x=x;
  tmpm->y=y;
  tmpm->DxDt=0;
  tmpm->DyDt=0;
  tmps=ode->next_spring;
  while (tmps!=&ode->next_spring)
  {
    if (tmps->end1==tmpm || tmps->end2==tmpm)
    {
      if (tmps->flags&SSF_NO_COMPRESSION)
        NullSpring(tmps,WIRE_PERCENT);
      else
        NullSpring(tmps,1.0);
    }
    tmps=tmps->next;
  }
}

U0 DelSpring(MySpring *tmps)
{
  QueRem(tmps);
  Free(tmps);
}

U0 DelMass(MyMass *tmpm)
{
  MySpring *tmps,*tmps1;
  tmps=ode->next_spring;
  while (tmps!=&ode->next_spring)
  {
    tmps1=tmps->next;
    if (tmps->end1==tmpm || tmps->end2==tmpm)
      DelSpring(tmps);
    tmps=tmps1;
  }
  QueRem(tmpm);
  Free(tmpm);
}

U0 DrawSpring(CDC *dc,MyMass *tmpm,I64 x,I64 y)
{
  switch (mode_bttn.state)
  {
  case MD_CONCRETE:
    dc->color=LTGRAY;
    dc->thick=2;
    break;
  case MD_STEEL:
    dc->color=DKGRAY;
    dc->thick=2;
    break;
  case MD_WIRE:
    dc->color=RED;
    dc->thick=1;
    break;
  }
  GrLine3(dc,tmpm->x,tmpm->y,0,x,y,0);
}

U0 PlaceSpring(MyMass *tmpm1,MyMass *tmpm2)
{
  MySpring *tmps=CAlloc(sizeof(MySpring));
  tmps->end1=tmpm1;
  tmps->end2=tmpm2;
  switch (mode_bttn.state)
  {
  case MD_CONCRETE:
    tmps->base_const          = 3.00*SPRING_SCALE;
    tmps->base_compression_strength=10.00*STRENGTH_SCALE;
    tmps->base_tensile_strength    = 0.35*STRENGTH_SCALE;
    tmps->base_cost                   = 0.30*COST_SCALE;
    NullSpring(tmps,1.0);
    tmps->color=LTGRAY;
    tmps->thick=2;
    break;
  case MD_STEEL:
    tmps->base_const          = 1.00*SPRING_SCALE;
    tmps->base_compression_strength= 1.00*STRENGTH_SCALE;
    tmps->base_tensile_strength    = 1.00*STRENGTH_SCALE;
    tmps->base_cost                   = 1.00*COST_SCALE;
    NullSpring(tmps,1.0);
    tmps->color=DKGRAY;
    tmps->thick=2;
    break;
  case MD_WIRE:
    tmps->base_const          = 0.25*SPRING_SCALE;
    tmps->base_compression_strength= 0.00;
    tmps->base_tensile_strength    = 0.50*STRENGTH_SCALE;
    tmps->base_cost                   = 0.10*COST_SCALE;
    NullSpring(tmps,WIRE_PERCENT);
    tmps->color=RED;
    tmps->thick=1;
    tmps->flags|=SSF_NO_COMPRESSION;
    break;
  }
  QueIns(tmps,ode->last_spring);
}

U0 AnimateTask(SpanAnimateStruct *a)
{
  MySpring *tmps,*tmps1;
  Bool old_run=FALSE;
  F64 f;
  while (TRUE)
  {
    tmps=ode->next_spring;
    while (tmps!=&ode->next_spring)
    {
      tmps1=tmps->next;
      f=tmps->f;
      if (f>0 &&  f>tmps->compression_strength &&
          !(tmps->flags&SSF_NO_COMPRESSION)||
          f<0 && -f>tmps->tensile_strength &&
          !(tmps->flags&SSF_NO_TENSION))
        tmps->flags|=SSF_INACTIVE;
      tmps=tmps1;
    }
    AdjustLoads(ode);
    Refresh; //CMathODE updated once per refresh.
    if (old_run!=run_bttn.state)
    {
      if (run_bttn.state)
      {
        if (!a->elapsed_t || !a->saved_ode)
        {
          Free(a->saved_ode);
          a->saved_ode=SpanSave(ode);
        }
        a->start_wall_t=tS;
        ODEPause(ode,OFF);
      }
      else
      {
        ODEPause(ode);
        a->elapsed_t+=tS-a->start_wall_t;
      }
      old_run=run_bttn.state;
    }
  }
}

U0 Init(SpanAnimateStruct *a)
{
  SpanDel(ode);
  ode=SpanNew;
  run_bttn.state=0;
  Refresh(2); //Allow stop to reg in animate task.
  if (a->saved_ode)
    SpanLoad(ode,a->saved_ode);
  else
    SpanBridge1Init(ode);
  a->elapsed_t=0;
  cursor_mass=NULL;
}

U0 SongTask(I64)
{//Song by Terry A. Davis
  Fs->task_end_cb=&SndTaskEndCB;
  MusicSettingsRst;
  music.tempo= 3.636;
  music.stacatto_factor= 0.902;
  while (TRUE)
  {
    Play("5q.EeDqED4G5DhE");
    Play("5q.EeDqED4G5DhE");
    Play("5q.FeEFEqF4G5EhF");
    Play("5q.FeEFEqF4G5EhF");
  }
}

U0 Span()
{
  I64 msg_code,arg1,arg2;
  MyMass   *tmpm1=NULL,*tmpm2=NULL;
  MySpring *tmps;
  CCtrl *bt_run,*bt_mode;
  U8 *src;
  CDC *dc=DCAlias;
  SettingsPush; //See SettingsPush
  Fs->text_attr=BROWN<<4+BLACK;
  AutoComplete;
  WinBorder;
  WinMax;
  DocCursor;
  Fs->song_task=Spawn(&SongTask,NULL,"Song",,Fs);
  bt_run =CtrlBttnNew(0,0,                80,,
                      2,"Stopped\0Running\0",run_colors,&run_bttn);
  bt_mode=CtrlBttnNew(0,3.0*FONT_HEIGHT,80,,
                      MD_MODES_NUM,Define("ST_SPAN_MODES"),mode_colors,&mode_bttn);
  a.saved_ode=NULL;
  Fs->win_inhibit|=WIG_DBL_CLICK;
  MenuPush(
  "File {"
  "  New(,CH_CTRLN);"
  "  Open(,CH_CTRLO);"
  "  SaveAs(,CH_CTRLA);"
  "  Abort(,CH_SHIFT_ESC);"
  "  Exit(,CH_ESC);"
  "}"
  "Play {"
  "  Restart(,'\n');"
  "  RunStop(,CH_SPACE);"
  "  Mass(,'m');"
  "  Concrete(,'c');"
  "  Steel(,'s');"
  "  Wire(,'w');"
  "  Move(,'v');"
  "  Delete(,'d');"
  "}"
  );
  ode=NULL;
  Init(&a);
  Fs->animate_task=Spawn(&AnimateTask,&a,"Animate",,Fs);
  Fs->draw_it=&DrawIt;
  PopUpOk(
  "Build a bridge to hold-up the\n"
  "red masses.  Test your design\n"
  "by pressing run/stop.\n\n"
  "The lowest cost bridge that\n"
  "stays standing wins.\n\n"
  "For a variation, try without\n"
  "using the center base point.\n"
  "\n"
  "Use\n"
  "\t$GREEN$'m'$FG$ass\n"
  "\t$GREEN$'c'$FG$oncrete\n"
  "\t$GREEN$'s'$FG$teel\n"
  "\t$GREEN$'w'$FG$ire\n"
  "\nto sel materials.\n");
  try
  {
    while (TRUE)
    {
      msg_code=GetMsg(&arg1,&arg2,1<<MSG_MS_L_DOWN|1<<MSG_MS_R_DOWN|
                      1<<MSG_MS_L_UP|1<<MSG_KEY_DOWN|1<<MSG_MS_MOVE);
      DCFill(dc);
      switch (msg_code)
      {
      case MSG_MS_L_DOWN:
        cursor_mass=tmpm1=tmpm2=NULL;
        switch (mode_bttn.state)
        {
        case MD_MASS:
          PlaceMass(arg1,arg2);
          break;
        case MD_CONCRETE:
        case MD_STEEL:
        case MD_WIRE:
          tmpm1=MassFind(ode,arg1,arg2);
          break;
        case MD_MOVE:
          if (run_bttn.state)
            cursor_mass=MassFind(ode,arg1,arg2);
          else if (tmpm1=MassFind(ode,arg1,arg2))
            MoveMass(tmpm1,arg1,arg2);
          break;
        case MD_DELETE:
          MassOrSpringFind(ode,&tmpm1,&tmps,arg1,arg2);
          if (tmpm1)
            DelMass(tmpm1);
          if (tmps)
            DelSpring(tmps);
          break;
        }
        break;
      case MSG_MS_L_UP:
        switch (mode_bttn.state)
        {
        case MD_CONCRETE:
        case MD_STEEL:
        case MD_WIRE:
          if (tmpm1 && (tmpm2=MassFind(ode,arg1,arg2)) && tmpm1!=tmpm2)
            PlaceSpring(tmpm1,tmpm2);
          break;
        case MD_MOVE:
          if (!run_bttn.state && tmpm1)
            MoveMass(tmpm1,arg1,arg2);
          break;
        }
        cursor_mass=tmpm1=tmpm2=NULL;
        break;
      case MSG_MS_MOVE:
        switch (mode_bttn.state)
        {
        case MD_MOVE:
          if (!run_bttn.state && tmpm1)
            MoveMass(tmpm1,arg1,arg2);
          break;
        case MD_CONCRETE:
        case MD_STEEL:
        case MD_WIRE:
          if (tmpm1)
          {
            DrawSpring(dc,tmpm1,arg1,arg2);
          }
          break;
        }
        break;
      case MSG_MS_R_DOWN:
        mode_bttn.state++;
        if (mode_bttn.state>=MD_MODES_NUM)
          mode_bttn.state=0;
        cursor_mass=tmpm1=tmpm2=NULL;
        break;
      case MSG_KEY_DOWN:
        switch (arg1)
        {
        case '\n':
          if (!SpanTime || !a.saved_ode)
          {
            Free(a.saved_ode);
            a.saved_ode=SpanSave(ode);
          }
          Init(&a);
          break;
        case CH_CTRLN:
          Free(a.saved_ode);
          a.saved_ode=NULL;
          Init(&a);
          break;
        case CH_CTRLO:
          if (src=SpanRead)
          {
            Free(a.saved_ode);
            a.saved_ode=src;
            Init(&a);
          }
          break;
        case CH_CTRLA:
          if (!SpanTime || !a.saved_ode)
          {
            Free(a.saved_ode);
            a.saved_ode=SpanSave(ode);
          }
          Init(&a);
          SpanWrite(ode);
          break;
        case CH_SPACE:
          run_bttn.state=!run_bttn.state;
          break;
        case 'c':
          mode_bttn.state=MD_CONCRETE;
          break;
        case 's':
          mode_bttn.state=MD_STEEL;
          break;
        case 'w':
          mode_bttn.state=MD_WIRE;
          break;
        case 'm':
          mode_bttn.state=MD_MASS;
          break;
        case 'v':
          mode_bttn.state=MD_MOVE;
          break;
        case 'd':
          mode_bttn.state=MD_DELETE;
          break;
        case CH_ESC:
          if (!SpanTime || !a.saved_ode)
          {
            Free(a.saved_ode);
            a.saved_ode=SpanSave(ode);
          }
          Init(&a);
          SpanWrite(ode);
        case CH_SHIFT_ESC:
          goto span_done;
        }
        break;
      }
    }
span_done: //Don't goto out of try
    GetMsg(,,1<<MSG_KEY_UP);
  }
  catch
    PutExcept;
  DocClear;
  SettingsPop;
  CtrlBttnDel(bt_run);
  CtrlBttnDel(bt_mode);
  SpanDel(ode);
  DCFill(dc);
  DCDel(dc);
  MenuPop;
}