Let's start our exploration of the new mapper functions with the 4-output mapper functions of 2, 3, and 4 arguments. In GeneXproTools these functions are represented by Map4A, Map4B, and Map4C, where the big letter at the end represents the number of arguments (A corresponds to 2, B to 3, and C to 4).

The Map4B and Map4C functions were the very first that I designed and tested (remember, I was trying to create functions with 3-4 discrete outputs that matched the argmin/argmax functions in performance), and both of them perform better than the argmax function of 4 arguments (98% hits for B, 99% for C, and 96% for the argmax).

I've already talked about the Map4B function in the post "Function Design: New 3-6 Output Functions" and as you can see from the code below, the Map4C function explores a similar design structure, with the difference that the middle point is in this case also determined independently by one of the arguments of the function.

The Map4A function is a late comer and is slightly inferior to the Map4B and C, but nonetheless performs slightly better than the argmax function of 4 arguments, even though it's only a function of 2 arguments (97% hits vs 96% for the argmax)!

Like I said, the Map4A function (and all mapper functions of the A series; I'll describe them in subsequent posts) was one of the latest mapper functions that I created. As they use only one argument for the map, I had to define a slack to generate the necessary number of intervals (here, for the Map4A function, we need to define 4 intervals for the four different outputs). And my first choice for the slack was 1000 but unfortunately it didn't work very well. I was very disappointed because I so wanted to create mapper functions of 2 arguments (functions of 2 arguments enjoy a special status in GeneXproTools because they require less bulky tree structures)! Just out of desperation I decided to try other values for the slack: 100, 10, and 1. And surprise: It made a huge difference, with 10 the best value for the Map4A function (a slack of 1 was the second best). (Who said all intervals are equal?)

For completeness I'm including here the C++ code for all the mappers of this series, including the Map4B function introduced in the post "Function Design: New 3-6 Output Functions":

// Map4A(x0,x1): 4-Output Mapper Function

const
double SLACK = 10.0;

double output = 0.0;

if (x[1] < (x[0] - SLACK))

output = 0.0;

else
if (x[1] >= (x[0] - SLACK) && x[1] < x[0])

output = 1.0;

else
if (x[1] >= x[0] && x[1] < (x[0] + SLACK))

output = 2.0;

else
if (x[1] >= (x[0] + SLACK))

output = 3.0;

return output;

// Map4B(x0,x1,x2): 4-Output Mapper Function

// evaluate min(x0,x1), max(x0,x1) and midrange

double min = x[0];

double max = x[1];

if (min > x[1])

{

min = x[1];

max = x[0];

}

double midrange = (max + min)/2.0;

double output = 0.0;

if (x[2] < min)

output = 0.0;

else
if (x[2] >= min && x[2] < midrange)

output = 1.0;

else
if (x[2] >= midrange && x[2] < max)

output = 2.0;

else
if (x[2] >= max)

output = 3.0;

return output;

// Map4C(x0,x1,x2,x3): 4-Output Mapper Function

// evaluate min(x,y,z), max(x,y,z) and midleValue(x,y,z)

//

// evaluate min(x,y,z)

double min = x[0];

int argmin = 0;

if (min > x[1])

{

min = x[1];

argmin = 1;

}

if (min > x[2])

{

min = x[2];

argmin = 2;

}

// evaluate max(x,y,z)

double max = x[0];

int argmax = 0;

if (max < x[1])

{

max = x[1];

argmax = 1;

}

if (max < x[2])

{

max = x[2];

argmax = 2;

}

// evaluate midleValue(x,y,z)

double midleValue = x[2];

if (0 != argmin && 0 != argmax)

midleValue = x[0];

if (1 != argmin && 1 != argmax)

midleValue = x[1];

double output = 0.0;

if (x[3] < min)

output = 0.0;

else
if (x[3] >= min && x[3] < midleValue)

output = 1.0;

else
if (x[3] >= midleValue && x[3] < max)

output = 2.0;

else
if (x[3] >= max)

output = 3.0;

return output;

In the next post I'll describe the 5-output mapper functions of 2, 3, and 4 arguments.

## Comments

There are currently no comments on this article.

## Comment

your_ip_is_blacklisted_by sbl.spamhaus.org