summaryrefslogtreecommitdiff
path: root/benchmarks/CUDA/MUM/mummergpu_gold.cpp
diff options
context:
space:
mode:
Diffstat (limited to 'benchmarks/CUDA/MUM/mummergpu_gold.cpp')
-rw-r--r--benchmarks/CUDA/MUM/mummergpu_gold.cpp440
1 files changed, 0 insertions, 440 deletions
diff --git a/benchmarks/CUDA/MUM/mummergpu_gold.cpp b/benchmarks/CUDA/MUM/mummergpu_gold.cpp
deleted file mode 100644
index e5507d3..0000000
--- a/benchmarks/CUDA/MUM/mummergpu_gold.cpp
+++ /dev/null
@@ -1,440 +0,0 @@
-#include <stdio.h>
-#include <string.h>
-
-#include <fstream>
-#define ulong4 uint32_t
-#include "mummergpu.h"
-
-static const int maxdim = 4096;
-
-#define __VERBOSE___no
-
-#ifdef __VERBOSE__
-#define XPRINTF(...) printf(__VA_ARGS__)
-#else
-#define XPRINTF(...) do{}while(0)
-#endif
-
-inline void set_result(const TextureAddress& cur,
- MatchCoord* result,
- int edge_match_length,
- int qry_match_len,
- int min_match_len,
- int rc)
-{
- if (qry_match_len > min_match_len)
- {
- int blocky = cur.y & 0x1F;
- int bigy = cur.y >> 5;
- int bigx = (cur.x << 5) + blocky;
- int nodeid = bigx + (bigy << 17);
-
- edge_match_length |= rc;
- result->node = nodeid;
- result->edge_match_length = edge_match_length;
- }
-}
-
-int kernel_gold(int qryid,
- MatchResults* results,
- char* refstr,
- char* queries,
- int* queryAddrs,
- int* queryLengths,
- PixelOfNode* nodeTexture,
- PixelOfChildren* childrenTexture,
- int numQueries,
- int min_match_len)
-{
- if (qryid >= numQueries) { return 0; }
-
- XPRINTF("> qryid: %d\n", qryid);
-
- int qlen = queryLengths[qryid];
-
- // start at root for first query character
- TextureAddress cur;
- cur.data = 0;
- int mustmatch = 0;
- int qry_match_len = 0;
-
- int qryAddr=queryAddrs[qryid];
- MatchCoord * result = results->h_match_coords + qryAddr - (qryid * ( min_match_len + 1));
- queries += qryAddr;
-
- int last = qlen - min_match_len;
-
- for (int qrystart = 0;
- qrystart <= last;
- qrystart++, result++, queries++)
- {
- XPRINTF("qry: %s\n", queries + 1);
-
- PixelOfNode node;
- TextureAddress prev;
-
- if ((cur.data == 0) || (qry_match_len < 1))
- {
- // start at root of tree
- cur.x = 0; cur.y = 1;
- qry_match_len = 1;
- mustmatch = 0;
- }
-
- char c = queries[qry_match_len];
-
- XPRINTF("In node (%d,%d): starting with %c [%d] => \n", cur.x, cur.y, c, qry_match_len);
-
- int refpos = 0;
- while ((c != '\0'))
- {
- XPRINTF("Next edge to follow: %c (%d)\n", c, qry_match_len);
-
- PixelOfChildren children;
- children = *(childrenTexture + (cur.x) + ((cur.y) * maxdim));
- prev = cur;
-
- switch(c)
- {
- case 'A': cur=children.children[0]; break;
- case 'C': cur=children.children[1]; break;
- case 'G': cur=children.children[2]; break;
- case 'T': cur=children.children[3]; break;
- default: cur.data = 0; break;
- };
-
- XPRINTF(" In node: (%d,%d)\n", cur.x, cur.y);
-
- // No edge to follow out of the node
- if (cur.data == 0)
- {
- XPRINTF(" no edge\n");
- set_result(prev, result, 0, qry_match_len, min_match_len, FORWARD);
-
- qry_match_len -= 1;
- mustmatch = 0;
-
- goto NEXT_SUBSTRING;
- }
-
- {
- unsigned short xval = cur.data & 0xFFFF;
- unsigned short yval = (cur.data & 0xFFFF0000) >> 16;
- node = *(nodeTexture + (xval) + (yval * maxdim));
- }
-
- XPRINTF(" Edge coordinates: %d - %d\n", node.start, node.end);
-
- if (mustmatch)
- {
- int edgelen = node.end - node.start+1;
- if (mustmatch >= edgelen)
- {
- XPRINTF(" mustmatch(%d) >= edgelen(%d), skipping edge\n", mustmatch, edgelen);
-
- refpos = node.end+1;
- qry_match_len += edgelen;
- mustmatch -= edgelen;
- }
- else
- {
- XPRINTF(" mustmatch(%d) < edgelen(%d), skipping to:%d\n",
- mustmatch, edgelen, node.start+mustmatch);
-
- qry_match_len += mustmatch;
- refpos = node.start + mustmatch;
- mustmatch = 0;
- }
- }
- else
- {
- // Try to walk the edge, the first char definitely matches
- qry_match_len++;
- refpos = node.start+1;
- }
-
- c = queries[qry_match_len];
-
- while (refpos <= node.end && c != '\0')
- {
- //char r = getRef(refpos);
- char r = refstr[refpos];
- XPRINTF(" Edge cmp ref: %d %c, qry: %d %c\n", refpos, r, qry_match_len, c);
-
- if (r != c)
- {
- // mismatch on edge
- XPRINTF("mismatch on edge: %d, edge_pos: %d\n", qry_match_len,refpos - (node.start));
- goto RECORD_RESULT;
- }
-
- qry_match_len++;
- refpos++;
- c = queries[qry_match_len];
- }
- }
-
- XPRINTF("end of string\n");
-
- RECORD_RESULT:
-
- set_result(cur, result, refpos - node.start, qry_match_len, min_match_len, FORWARD);
- mustmatch = refpos - node.start;
- qry_match_len -= mustmatch + 1;
-
- NEXT_SUBSTRING:
-
- //node.data = texfetch(nodetex, prev.x, prev.y);
- node = *(nodeTexture + (prev.x) + ((prev.y) * maxdim));
- cur = node.suffix;
-
- XPRINTF(" following suffix link. mustmatch:%d qry_match_len:%d sl:(%d,%d)\n",
- mustmatch, qry_match_len, cur.x, cur.y);
-
- do {} while(0);
- }
-
- return 0 ;
-}
-
-
-
-inline char rc(char c)
-{
- switch(c)
- {
- case 'A': return 'T';
- case 'C': return 'G';
- case 'G': return 'C';
- case 'T': return 'A';
- case 'q': return '\0';
- default: return c;
- };
-}
-
-///////////////////////////////////////
-//// Compute reverse substring matches
-///////////////////////////////////////
-
-
-void
-rc_kernel_gold(int qryid,
- MatchResults* results,
- char* refstr,
- char* queries,
- const int* queryAddrs,
- const int* queryLengths,
- PixelOfNode* nodeTexture,
- PixelOfChildren* childrenTexture,
- const int numQueries,
- const int min_match_len)
-{
-
- if (qryid >= numQueries) { return; }
-
- int qlen = queryLengths[qryid];
-
- XPRINTF("> rc qryid: %d\n", qryid);
-
- queries++; // skip the 'q' character
-
-
- // start at root for first query character
- TextureAddress cur;
-
- int mustmatch = 0;
- int qry_match_len = 0;
-
- int qryAddr=queryAddrs[qryid];
- MatchCoord * result = results->h_match_coords + qryAddr - (qryid * ( min_match_len + 1));
- queries += qryAddr;
-
- for (int qrystart = qlen;
- qrystart >= min_match_len ;
- qrystart--, result++)
- {
- #ifdef VERBOSE
- queries[qrystart] = '\0';
- XPRINTF("qry: ", queries);
- for (int j = qrystart-1; j >= 0; j--)
- { XPRINTF("%c", rc(queries[j])); }
- XPRINTF("\n");
- #endif
-
- PixelOfNode node;
- TextureAddress prev;
-
- if (((cur.data == 0)) || (qry_match_len < 1))
- {
- // start at root of tree
- cur.x = 0; cur.y = 1;
- qry_match_len = 1;
- mustmatch = 0;
- }
-
- char c = rc(queries[qrystart-qry_match_len]);
-
- XPRINTF("In node (%d,%d): starting with %c [%d] => \n", cur.x, cur.y, c, qry_match_len);
-
- int refpos = 0;
- while ((c != '\0'))
- {
- XPRINTF("Next edge to follow: %c (%d)\n", c, qry_match_len);
-
- PixelOfChildren children;
- children = *(childrenTexture + (cur.x) + ((cur.y) * maxdim));
- prev = cur;
-
- switch(c)
- {
- case 'A': cur=children.children[0]; break;
- case 'C': cur=children.children[1]; break;
- case 'G': cur=children.children[2]; break;
- case 'T': cur=children.children[3]; break;
- default: cur.data = 0; break;
- };
-
- XPRINTF(" In node: (%d,%d)\n", cur.x, cur.y);
-
- // No edge to follow out of the node
- if (cur.data == 0)
- {
- XPRINTF(" no edge\n");
- set_result(prev, result, 0, qry_match_len, min_match_len,
- REVERSE);
-
- qry_match_len -= 1;
- mustmatch = 0;
-
- goto NEXT_SUBSTRING;
- }
-
- {
- unsigned short xval = cur.data & 0xFFFF;
- unsigned short yval = (cur.data & 0xFFFF0000) >> 16;
- node = *(nodeTexture + (xval) + (yval * maxdim));
- }
-
- XPRINTF(" Edge coordinates: %d - %d\n", node.start, node.end);
-
- if (mustmatch)
- {
- int edgelen = node.end - node.start+1;
- if (mustmatch >= edgelen)
- {
- XPRINTF(" mustmatch(%d) >= edgelen(%d), skipping edge\n", mustmatch, edgelen);
-
- refpos = node.end+1;
- qry_match_len += edgelen;
- mustmatch -= edgelen;
- }
- else
- {
- XPRINTF(" mustmatch(%d) < edgelen(%d), skipping to:%d\n",
- mustmatch, edgelen, node.start+mustmatch);
-
- qry_match_len += mustmatch;
- refpos = node.start + mustmatch;
- mustmatch = 0;
- }
- }
- else
- {
- // Try to walk the edge, the first char definitely matches
- qry_match_len++;
- refpos = node.start+1;
- }
-
- c = rc(queries[qrystart-qry_match_len]);
-
- while (refpos <= node.end && c != '\0')
- {
- //char r = getRef(refstr, refpos);
- char r = refstr[refpos];
- XPRINTF(" Edge cmp ref: %d %c, qry: %d %c\n", refpos, r, qry_match_len, c);
-
- if (r != c)
- {
- // mismatch on edge
- XPRINTF("mismatch on edge: %d, edge_pos: %d\n", qry_match_len,refpos - (node.start));
- goto RECORD_RESULT;
- }
-
- qry_match_len++;
- refpos++;
- c = rc(queries[qrystart-qry_match_len]);
- }
- }
-
- XPRINTF("end of string\n");
-
- RECORD_RESULT:
-
- set_result(cur, result, refpos - node.start, qry_match_len,
- min_match_len, REVERSE);
-
- mustmatch = refpos - node.start;
- qry_match_len -= mustmatch + 1;
-
- NEXT_SUBSTRING:
-
- node = *(nodeTexture + (prev.x) + ((prev.y) * maxdim));
- cur = node.suffix;
-
- XPRINTF(" following suffix link. mustmatch:%d qry_match_len:%d sl:(%d,%d)\n",
- mustmatch, qry_match_len, cur.x, cur.y);
-
- do {} while(0);
- }
-
- return;
-}
-
-
-extern "C"
-void computeGold(MatchResults* results,
- char* refstr,
- char* queries,
- int* queryAddrs,
- int* queryLengths,
- PixelOfNode* nodeTexture,
- PixelOfChildren* childrenTexture,
- int numQueries,
- int match_length,
- int rc)
-{
-
- if (rc == REVERSE)
- {
- for (int i = 0; i < numQueries; ++i)
- {
- rc_kernel_gold(i,
- results,
- refstr,
- queries,
- queryAddrs,
- queryLengths,
- nodeTexture,
- childrenTexture,
- numQueries,
- match_length);
- }
- }
- else
- {
- for (int i = 0; i < numQueries; ++i)
- {
- kernel_gold(i,
- results,
- refstr,
- queries,
- queryAddrs,
- queryLengths,
- nodeTexture,
- childrenTexture,
- numQueries,
- match_length);
- }
- }
-}
-
generated by cgit v1.3 (git 2.53.0) at 2026-07-10 23:26:45 +0000