Skip to content

Trying to use std::pow with std::valarray for compute_gamma() #1

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
wants to merge 1 commit into
base: master
Choose a base branch
from
Open

Trying to use std::pow with std::valarray for compute_gamma() #1

Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
177 changes: 139 additions & 38 deletions isp-hls/src/xf_isp_tb.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -210,7 +210,7 @@ static void MultiPixelAXIvideo2Mat_yuv(OutVideoStrm_t& AXI_video_strm, cv::Mat&
}
}
void compute_gamma(float r_g, float g_g, float b_g, uchar gamma_lut[256 * 3]) {
float gamma_inv[256] = {
float gamma_inv[256 * 3] = {
0.000000, 0.003922, 0.007843, 0.011765, 0.015686, 0.019608, 0.023529, 0.027451, 0.031373, 0.035294, 0.039216,
0.043137, 0.047059, 0.050980, 0.054902, 0.058824, 0.062745, 0.066667, 0.070588, 0.074510, 0.078431, 0.082353,
0.086275, 0.090196, 0.094118, 0.098039, 0.101961, 0.105882, 0.109804, 0.113725, 0.117647, 0.121569, 0.125490,
Expand All @@ -234,43 +234,144 @@ void compute_gamma(float r_g, float g_g, float b_g, uchar gamma_lut[256 * 3]) {
0.862745, 0.866667, 0.870588, 0.874510, 0.878431, 0.882353, 0.886275, 0.890196, 0.894118, 0.898039, 0.901961,
0.905882, 0.909804, 0.913725, 0.917647, 0.921569, 0.925490, 0.929412, 0.933333, 0.937255, 0.941176, 0.945098,
0.949020, 0.952941, 0.956863, 0.960784, 0.964706, 0.968627, 0.972549, 0.976471, 0.980392, 0.984314, 0.988235,
0.992157, 0.996078, 1.000000};

unsigned char gam_r = 0, gam_g = 0, gam_b = 0;

for (int i = 0; i < 256; ++i) {
float r_inv = (float)1 / r_g;
float g_inv = (float)1 / g_g;
float b_inv = (float)1 / b_g;
float powval_r = (float)std::pow(gamma_inv[i], r_inv);
short tempgamma_r = (powval_r * 255.0);

if (tempgamma_r > 255) {
gam_r = 255;
} else {
gam_r = tempgamma_r;
}

float powval_g = (float)std::pow(gamma_inv[i], g_inv);
short tempgamma_g = (powval_g * 255.0);

if (tempgamma_g > 255) {
gam_g = 255;
} else {
gam_g = tempgamma_g;
}

float powval_b = (float)std::pow(gamma_inv[i], b_inv);
short tempgamma_b = (powval_b * 255.0);

if (tempgamma_b > 255) {
gam_b = 255;
} else {
gam_b = tempgamma_b;
}
gamma_lut[i] = gam_r;
gamma_lut[i + 256] = gam_g;
gamma_lut[i + 512] = gam_b;
0.992157, 0.996078, 1.000000,

0.000000, 0.003922, 0.007843, 0.011765, 0.015686, 0.019608, 0.023529, 0.027451, 0.031373, 0.035294, 0.039216,
0.043137, 0.047059, 0.050980, 0.054902, 0.058824, 0.062745, 0.066667, 0.070588, 0.074510, 0.078431, 0.082353,
0.086275, 0.090196, 0.094118, 0.098039, 0.101961, 0.105882, 0.109804, 0.113725, 0.117647, 0.121569, 0.125490,
0.129412, 0.133333, 0.137255, 0.141176, 0.145098, 0.149020, 0.152941, 0.156863, 0.160784, 0.164706, 0.168627,
0.172549, 0.176471, 0.180392, 0.184314, 0.188235, 0.192157, 0.196078, 0.200000, 0.203922, 0.207843, 0.211765,
0.215686, 0.219608, 0.223529, 0.227451, 0.231373, 0.235294, 0.239216, 0.243137, 0.247059, 0.250980, 0.254902,
0.258824, 0.262745, 0.266667, 0.270588, 0.274510, 0.278431, 0.282353, 0.286275, 0.290196, 0.294118, 0.298039,
0.301961, 0.305882, 0.309804, 0.313725, 0.317647, 0.321569, 0.325490, 0.329412, 0.333333, 0.337255, 0.341176,
0.345098, 0.349020, 0.352941, 0.356863, 0.360784, 0.364706, 0.368627, 0.372549, 0.376471, 0.380392, 0.384314,
0.388235, 0.392157, 0.396078, 0.400000, 0.403922, 0.407843, 0.411765, 0.415686, 0.419608, 0.423529, 0.427451,
0.431373, 0.435294, 0.439216, 0.443137, 0.447059, 0.450980, 0.454902, 0.458824, 0.462745, 0.466667, 0.470588,
0.474510, 0.478431, 0.482353, 0.486275, 0.490196, 0.494118, 0.498039, 0.501961, 0.505882, 0.509804, 0.513725,
0.517647, 0.521569, 0.525490, 0.529412, 0.533333, 0.537255, 0.541176, 0.545098, 0.549020, 0.552941, 0.556863,
0.560784, 0.564706, 0.568627, 0.572549, 0.576471, 0.580392, 0.584314, 0.588235, 0.592157, 0.596078, 0.600000,
0.603922, 0.607843, 0.611765, 0.615686, 0.619608, 0.623529, 0.627451, 0.631373, 0.635294, 0.639216, 0.643137,
0.647059, 0.650980, 0.654902, 0.658824, 0.662745, 0.666667, 0.670588, 0.674510, 0.678431, 0.682353, 0.686275,
0.690196, 0.694118, 0.698039, 0.701961, 0.705882, 0.709804, 0.713725, 0.717647, 0.721569, 0.725490, 0.729412,
0.733333, 0.737255, 0.741176, 0.745098, 0.749020, 0.752941, 0.756863, 0.760784, 0.764706, 0.768627, 0.772549,
0.776471, 0.780392, 0.784314, 0.788235, 0.792157, 0.796078, 0.800000, 0.803922, 0.807843, 0.811765, 0.815686,
0.819608, 0.823529, 0.827451, 0.831373, 0.835294, 0.839216, 0.843137, 0.847059, 0.850980, 0.854902, 0.858824,
0.862745, 0.866667, 0.870588, 0.874510, 0.878431, 0.882353, 0.886275, 0.890196, 0.894118, 0.898039, 0.901961,
0.905882, 0.909804, 0.913725, 0.917647, 0.921569, 0.925490, 0.929412, 0.933333, 0.937255, 0.941176, 0.945098,
0.949020, 0.952941, 0.956863, 0.960784, 0.964706, 0.968627, 0.972549, 0.976471, 0.980392, 0.984314, 0.988235,
0.992157, 0.996078, 1.000000,

0.000000, 0.003922, 0.007843, 0.011765, 0.015686, 0.019608, 0.023529, 0.027451, 0.031373, 0.035294, 0.039216,
0.043137, 0.047059, 0.050980, 0.054902, 0.058824, 0.062745, 0.066667, 0.070588, 0.074510, 0.078431, 0.082353,
0.086275, 0.090196, 0.094118, 0.098039, 0.101961, 0.105882, 0.109804, 0.113725, 0.117647, 0.121569, 0.125490,
0.129412, 0.133333, 0.137255, 0.141176, 0.145098, 0.149020, 0.152941, 0.156863, 0.160784, 0.164706, 0.168627,
0.172549, 0.176471, 0.180392, 0.184314, 0.188235, 0.192157, 0.196078, 0.200000, 0.203922, 0.207843, 0.211765,
0.215686, 0.219608, 0.223529, 0.227451, 0.231373, 0.235294, 0.239216, 0.243137, 0.247059, 0.250980, 0.254902,
0.258824, 0.262745, 0.266667, 0.270588, 0.274510, 0.278431, 0.282353, 0.286275, 0.290196, 0.294118, 0.298039,
0.301961, 0.305882, 0.309804, 0.313725, 0.317647, 0.321569, 0.325490, 0.329412, 0.333333, 0.337255, 0.341176,
0.345098, 0.349020, 0.352941, 0.356863, 0.360784, 0.364706, 0.368627, 0.372549, 0.376471, 0.380392, 0.384314,
0.388235, 0.392157, 0.396078, 0.400000, 0.403922, 0.407843, 0.411765, 0.415686, 0.419608, 0.423529, 0.427451,
0.431373, 0.435294, 0.439216, 0.443137, 0.447059, 0.450980, 0.454902, 0.458824, 0.462745, 0.466667, 0.470588,
0.474510, 0.478431, 0.482353, 0.486275, 0.490196, 0.494118, 0.498039, 0.501961, 0.505882, 0.509804, 0.513725,
0.517647, 0.521569, 0.525490, 0.529412, 0.533333, 0.537255, 0.541176, 0.545098, 0.549020, 0.552941, 0.556863,
0.560784, 0.564706, 0.568627, 0.572549, 0.576471, 0.580392, 0.584314, 0.588235, 0.592157, 0.596078, 0.600000,
0.603922, 0.607843, 0.611765, 0.615686, 0.619608, 0.623529, 0.627451, 0.631373, 0.635294, 0.639216, 0.643137,
0.647059, 0.650980, 0.654902, 0.658824, 0.662745, 0.666667, 0.670588, 0.674510, 0.678431, 0.682353, 0.686275,
0.690196, 0.694118, 0.698039, 0.701961, 0.705882, 0.709804, 0.713725, 0.717647, 0.721569, 0.725490, 0.729412,
0.733333, 0.737255, 0.741176, 0.745098, 0.749020, 0.752941, 0.756863, 0.760784, 0.764706, 0.768627, 0.772549,
0.776471, 0.780392, 0.784314, 0.788235, 0.792157, 0.796078, 0.800000, 0.803922, 0.807843, 0.811765, 0.815686,
0.819608, 0.823529, 0.827451, 0.831373, 0.835294, 0.839216, 0.843137, 0.847059, 0.850980, 0.854902, 0.858824,
0.862745, 0.866667, 0.870588, 0.874510, 0.878431, 0.882353, 0.886275, 0.890196, 0.894118, 0.898039, 0.901961,
0.905882, 0.909804, 0.913725, 0.917647, 0.921569, 0.925490, 0.929412, 0.933333, 0.937255, 0.941176, 0.945098,
0.949020, 0.952941, 0.956863, 0.960784, 0.964706, 0.968627, 0.972549, 0.976471, 0.980392, 0.984314, 0.988235,
0.992157, 0.996078, 1.000000
};

// std::valarray<float> g = { (float)1 / r_g, (float)1 / g_g, (float)1 / b_g };
float r_g_inv = (float)1 / r_g;
float g_g_inv = (float)1 / g_g;
float b_g_inv = (float)1 / b_g;

float g_inv[256 * 3] = {
r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv,
r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv,
r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv,
r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv,
r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv,
r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv,
r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv,
r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv,
r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv,
r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv,
r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv,
r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv,
r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv,
r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv,
r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv,
r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv,
r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv,
r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv,
r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv,
r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv,
r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv,
r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv,
r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv, r_g_inv,
r_g_inv, r_g_inv, r_g_inv,

g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv,
g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv,
g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv,
g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv,
g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv,
g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv,
g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv,
g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv,
g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv,
g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv,
g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv,
g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv,
g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv,
g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv,
g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv,
g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv,
g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv,
g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv,
g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv,
g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv,
g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv,
g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv,
g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv, g_g_inv,
g_g_inv, g_g_inv, g_g_inv,

b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv,
b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv,
b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv,
b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv,
b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv,
b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv,
b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv,
b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv,
b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv,
b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv,
b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv,
b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv,
b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv,
b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv,
b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv,
b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv,
b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv,
b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv,
b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv,
b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv,
b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv,
b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv,
b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv, b_g_inv,
b_g_inv, b_g_inv, b_g_inv
};

std::valarray<float> gamma_lut_f = std::pow(gamma_inv, g_inv);
for (int i = 0; i < 256 * 3; i++) {
gamma_lut[i] = std::min((uchar)255, (uchar)(gamma_lut_f[i] * 255));
}
}
int main(int argc, char** argv) {
Expand Down