image.cpp 5.0 KB
Newer Older
Z
zhangyang 已提交
1 2 3 4 5 6 7 8 9 10 11 12
/* Copyright (c) 2018 PaddlePaddle Authors. All Rights Reserved.

Licensed under the Apache License, Version 2.0 (the "License");
you may not use this file except in compliance with the License.
You may obtain a copy of the License at

    http://www.apache.org/licenses/LICENSE-2.0

Unless required by applicable law or agreed to in writing, software
distributed under the License is distributed on an "AS IS" BASIS,
WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
See the License for the specific language governing permissions and
Z
zhangyang 已提交
13
limitations under the License. */
Z
zhangyang 已提交
14

Z
zhangyang 已提交
15
#include "fpga/V1/image.h"
Z
zhangyang 已提交
16
#include <memory.h>
17
#include <algorithm>
Z
zhangyang 已提交
18
#include "fpga/common/fpga_common.h"
Z
zhangyang 已提交
19 20 21 22 23 24 25 26

namespace paddle_mobile {
namespace fpga {
namespace image {

void convert_to_hwc(float **data_in, int channel, int height, int width) {
  float *tmp = *data_in;
  float *data_tmp =
27
      (float *)fpga_malloc(channel * height * width * sizeof(float));  // NOLINT
Z
zhangyang 已提交
28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45
  int64_t amount_per_row = width * channel;
  for (int c = 0; c < channel; c++) {
    for (int h = 0; h < height; h++) {
      int64_t offset_height = h * amount_per_row;
      for (int w = 0; w < width; w++) {
        *(data_tmp + offset_height + w * channel + c) = *((*data_in)++);
      }
    }
  }
  *data_in = data_tmp;
  fpga_free(tmp);
}

void align_element_conv(float **data_in, int height, int cw) {
  int h = 0;
  int align_cw = align_to_x(cw, IMAGE_ALIGNMENT);
  if (align_cw != cw) {
    float *tmp = *data_in;
46 47
    float *data_tmp =
        (float *)fpga_malloc(height * align_cw * sizeof(float));  // NOLINT
Z
zhangyang 已提交
48 49 50 51

    memset(data_tmp, 0, height * align_cw * sizeof(float));

    for (h = 0; h < height; h++) {
52 53
      memcpy((void *)(data_tmp + h * align_cw),  // NOLINT
             (void *)(*data_in + h * cw),        // NOLINT
Z
zhangyang 已提交
54 55 56 57 58 59 60 61 62 63 64
             cw * sizeof(float));
    }

    *data_in = data_tmp;
    fpga_free(tmp);
  }
}

void format_image(float **data_in, int channel, int height, int width) {
  convert_to_hwc(data_in, channel, height, width);
  align_element_conv(data_in, height, channel * width);
65 66
  fpga_flush(*data_in, align_to_x(channel * width, IMAGE_ALIGNMENT) * height *
                           sizeof(float));
Z
zhangyang 已提交
67 68
}

Z
zhangyang 已提交
69 70
void concat_images(int16_t **images_in, float **scales_in, void *image_out,
                   float *scale_out, int image_num, uint32_t *channel_num,
71 72 73 74 75 76 77 78 79
                   int height, int width) {
  int i = 0;
  int j = 0;
  int k = 0;
  int each_out_line_channel = 0;
  int align_each_out_area_cw = 0;
  int align_each_in_area_cw = 0;
  int align_each_out_area_cw_differ = 0;
  int tmp_channel = 0;
80 81
  scale_out[0] = 0.0;
  scale_out[1] = 0.0;
82 83
  for (i = 0; i < image_num; i++) {
    each_out_line_channel += channel_num[i];
84 85 86 87 88
    scale_out[0] = std::max(*scale_out, scales_in[i][0]);
    fpga_invalidate(images_in[i],
                    height *
                        align_to_x(channel_num[i] * width, IMAGE_ALIGNMENT) *
                        sizeof(int16_t));
89
  }
90
  scale_out[1] = 1 / scale_out[0];
91 92 93 94 95 96 97 98 99 100
  align_each_out_area_cw =
      align_to_x(each_out_line_channel * width, IMAGE_ALIGNMENT);
  align_each_out_area_cw_differ =
      align_each_out_area_cw - each_out_line_channel * width;

  for (k = 0; k < height; k++) {
    for (j = 0; j < width; j++) {
      for (i = 0; i < image_num; i++) {
        align_each_in_area_cw =
            align_to_x(channel_num[i] * width, IMAGE_ALIGNMENT);
101
        memcpy((int16_t *)image_out + tmp_channel +  // NOLINT
102 103 104 105 106 107 108 109
                   k * align_each_out_area_cw_differ,
               images_in[i] + j * channel_num[i] + k * align_each_in_area_cw,
               channel_num[i] * sizeof(int16_t));

        tmp_channel += channel_num[i];
      }
    }
  }
110 111

  fpga_flush(image_out, height * align_each_out_area_cw * sizeof(int16_t));
112
}
Z
zhangyang 已提交
113

114 115 116
void split_image(int16_t *image_in, const float *scale_in, void **images_out,
                 float **scales_out, int image_num,
                 const uint32_t *channel_nums, int height, int width) {
117 118 119 120 121 122
  int total_channel = 0;
  for (int i = 0; i < image_num; i++) {
    scales_out[i][0] = scale_in[0];
    scales_out[i][1] = scale_in[1];
    total_channel += channel_nums[i];
  }
123 124
  int element_num = height * align_to_x(width * total_channel, IMAGE_ALIGNMENT);
  fpga_invalidate(image_in, element_num * sizeof(int16_t));
125

126
  int src_offset = 0, des_offset = 0;
127
  for (int h = 0; h < height; h++) {
128 129 130 131 132 133 134 135 136 137
    for (int w = 0; w < width; w++) {
      src_offset = h * align_to_x(total_channel * width, IMAGE_ALIGNMENT) +
                   w * total_channel;
      for (int i = 0; i < image_num; i++) {
        des_offset = h * align_to_x(channel_nums[i] * width, IMAGE_ALIGNMENT) +
                     w * channel_nums[i];
        memcpy((int16_t *)images_out[i] + des_offset, image_in + src_offset,
               channel_nums[i] * sizeof(int16_t));
        src_offset += channel_nums[i];
      }
138 139
    }
  }
140 141 142 143 144

  for (int i = 0; i < image_num; i++) {
    element_num = height * align_to_x(width * channel_nums[i], IMAGE_ALIGNMENT);
    fpga_flush(images_out[i], element_num * sizeof(int16_t));
  }
145 146
}

Z
zhangyang 已提交
147 148 149
}  // namespace image
}  // namespace fpga
}  // namespace paddle_mobile