U8 *master_bitmap;
I64 prime_range,my_mp_cnt,pending;

U0 PrimesJob(I64 i)
{
  I64 j,k,l=Sqrt(prime_range),
          lo=i*prime_range/my_mp_cnt,
          hi=(i+1)*prime_range/my_mp_cnt,
          lo2=lo+sys_cache_line_width,
          hi2=hi-sys_cache_line_width;
  if (lo2>hi2) lo2=hi2;
  for (j=2; j<=l; j++)
  {
    if (!Bt(master_bitmap,j))
    {
      if (j>=lo)
        k=j+j;
      else
      {
        k=lo-lo%j;
        if (k<lo)
          k+=j;
      }
      while (k<lo2)
      {
        LBts(master_bitmap,k);
        k+=j;
      }
      while (k<hi2)
      {
        Bts(master_bitmap,k);
        k+=j;
      }
      while (k<hi)
      {
        LBts(master_bitmap,k);
        k+=j;
      }
    }
  }
  lock pending--;
}

I64 Primes(I64 range,I64 _my_mp_cnt)
{
  I64 i,prime_cnt=0;
  F64 t0=tS,tf;
  prime_range=range;
  if (_my_mp_cnt<=mp_cnt)
    my_mp_cnt=_my_mp_cnt;
  else
    my_mp_cnt=mp_cnt;
  master_bitmap=CAlloc((prime_range+7)/8+1);
  Bts(master_bitmap,0);
  Bts(master_bitmap,1);
  pending=my_mp_cnt;
  for (i=0; i<my_mp_cnt; i++)
    JobQue(&PrimesJob,i,i);
  while (pending)
    Yield;
  tf=tS;
  for (i=0; i<prime_range; i++)
    if (!Bt(master_bitmap,i))
      prime_cnt++;
  for (i=MaxI64(prime_range-100,0); i<prime_range; i++)
    if (!Bt(master_bitmap,i))
      "%d ",i;
  "\n$RED$CPUs:%d PrimeRange:%,d PrimeCnt:%,d Time:%9.7,f$FG$\n"
  ,
  my_mp_cnt,prime_range,prime_cnt,tf-t0;
  Free(master_bitmap);
  return prime_cnt;
}

Primes(100,1);
Primes(100,mp_cnt);
Primes(1000000,1);
Primes(1000000,mp_cnt);