bias_scale.cpp 3.2 KB
Newer Older
Z
zhangyang 已提交
1 2 3 4 5 6 7 8 9 10 11 12 13 14
/* 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
limitations under the License. */

15
#include "fpga/bias_scale.h"
Z
zhangyang 已提交
16
#include <memory.h>
17
#include "fpga/api.h"
Z
zhangyang 已提交
18 19

namespace paddle_mobile {
Z
zhangyang 已提交
20 21 22 23
namespace fpga {
namespace bias_scale {

void align_element(float **data_in, int num_per_div_before_alignment, int num) {
Z
zhangyang 已提交
24
  int copynum = 0;
Z
zhangyang 已提交
25 26 27 28 29
  float *ptr_unaligned = *data_in;
  int div_num =
      (num + num_per_div_before_alignment - 1) / num_per_div_before_alignment;
  int num_per_div_after_alignment =
      align_to_x(num_per_div_before_alignment, BS_NUM_ALIGNMENT);
Z
zhangyang 已提交
30
  if (num_per_div_before_alignment == num_per_div_after_alignment) {
31 32
    return;
  }
Z
zhangyang 已提交
33 34
  int num_element =
      2 * div_num * num_per_div_after_alignment;  // including bias & scale
35 36
  float *ptr_aligned =
      (float *)fpga_malloc(num_element * sizeof(float));  // NOLINT
Z
zhangyang 已提交
37 38 39 40

  memset(ptr_aligned, 0, num_element * sizeof(float));

  for (int i = 0; i < div_num; i++) {
Z
zhangyang 已提交
41 42 43 44 45 46 47 48 49 50 51 52 53 54
    if (i == div_num - 1) {
      copynum = (num_per_div_after_alignment * div_num > num)
                    ? (num % num_per_div_after_alignment)
                    : (num_per_div_before_alignment);
    } else {
      copynum = num_per_div_before_alignment;
    }

    memcpy(ptr_aligned + i * num_per_div_after_alignment,
           ptr_unaligned + num_per_div_before_alignment * i,
           copynum * sizeof(float));
    memcpy(ptr_aligned + (div_num + i) * num_per_div_after_alignment,
           ptr_unaligned + num_per_div_before_alignment * i + num,
           copynum * sizeof(float));
Z
zhangyang 已提交
55 56 57 58 59 60 61 62 63 64 65
  }

  fpga_free(ptr_unaligned);
  *data_in = ptr_aligned;
}

void interleave(float **data_in, int num_after_alignment) {
  // num_after_alignment: number of bias after alignment

  float *ptr_uninterleaved = *data_in;
  float *ptr_interleaved =
66
      (float *)fpga_malloc(2 * num_after_alignment * sizeof(float));  // NOLINT
Z
zhangyang 已提交
67 68 69 70 71
  int num = num_after_alignment / 4;
  for (int i = 0; i < num; i++) {
    memcpy(ptr_interleaved + 8 * i, ptr_uninterleaved + 4 * i,
           4 * sizeof(float));
    memcpy(ptr_interleaved + 8 * i + 4,
Z
zhangyang 已提交
72
           ptr_uninterleaved + num_after_alignment + 4 * i, 4 * sizeof(float));
Z
zhangyang 已提交
73 74 75 76 77 78
  }

  fpga_free(ptr_uninterleaved);
  *data_in = ptr_interleaved;
}

Z
zhangyang 已提交
79 80 81 82 83 84 85
void format_bias_scale_array(float **bias_scale_array,
                             int element_num_per_division, int num) {
  align_element(bias_scale_array, element_num_per_division, num);
  int div_num = (num + element_num_per_division - 1) / element_num_per_division;
  int element_num_after_division =
      align_to_x(element_num_per_division, BS_NUM_ALIGNMENT);
  interleave(bias_scale_array, div_num * element_num_after_division);
86
  fpga_flush(*bias_scale_array, 2 * element_num_after_division * sizeof(float));
Z
zhangyang 已提交
87 88
}

Z
zhangyang 已提交
89 90 91
}  // namespace bias_scale
}  // namespace fpga
}  // namespace paddle_mobile