DNA C programing

Code:

---

#include <stdio.h>
#include <stdlib.h>
#include <time.h>

int main()
{

        time_t time_start, time_end, time_period;
        time_start = time(0);
        //time_start = time(NULL);

        FILE *fPtr;

        if ((fPtr = fopen("basePair_125kb.bin", "wb+")) == NULL)
        {
                printf("File could not be opened.\n");
                return -1;
        }
       
        srand(time(NULL));

        unsigned long basePair_0, basePair_1, basePair_2, basePair_32bit;
       
        int loop;
       
        for (loop = 0; loop < 7812; loop++)
        {
                basePair_0 = rand();
                basePair_1 = rand();
                basePair_2 = rand();
                basePair_32bit = (basePair_2 << 30) + (basePair_1 << 15) + basePair_0;
               
                fprintf(fPtr, "%.8X", basePair_32bit);

                if (loop % 1000 == 0)
                {
                        printf("loop count = %d\n", loop);
                }
        }

        basePair_0 = rand();
        basePair_1 = rand();
        basePair_2 = rand();

        basePair_32bit = (basePair_2 << 30) + (basePair_1 << 15) + basePair_0;
        basePair_32bit = basePair_32bit >> 16;

        printf("Random 15-bit basePair_0: %u\n", basePair_0);
        printf("Random 15-bit basePair_1: %u\n", basePair_1);
        printf("Random 2-bit basePair_2: %u\n", basePair_2);
        printf("Random 32-bit basePair_32bit: %u\n", basePair_32bit);
        printf("RAND_MAX = %u\n", RAND_MAX);
        printf("sizeof(long) = %d\n", sizeof(long));

        fprintf(fPtr, "%.4X", basePair_32bit);
       
        fclose(fPtr);

        time_end = time(0);
        //time_end = time(NULL);
        time_period = time_end - time_start;
       
        printf("Time to run the random DNA generation = %d\n", time_period);

        return 0;
}

---

hi guys.... afta consult our lecturer, we successfully done few parts of our assignment... somehow, we do not know how to convert it from hexadecimal form to binary form and finally to character form (ATGC)... anyone can suggest us how to convert?

http://cboard.cprogramming.com/c-pro...aming-dna.html
And what was wrong with continuing the previous thread?

I have no idea what you think you are doing, as you seem to assume that rand() generates a 15-bit number, which I'm sure some implementations do, but it's far from guaranteed.

--
Mats

in order to get the number into AGCT form take the number and convert it into base 4, then use that as an array index for a character array holding "AGCT"