google-research

1
// Copyright 2024 The Google Research Authors.
2
//
3
// Licensed under the Apache License, Version 2.0 (the "License");
4
// you may not use this file except in compliance with the License.
5
// You may obtain a copy of the License at
6
//
7
//     http://www.apache.org/licenses/LICENSE-2.0
8
//
9
// Unless required by applicable law or agreed to in writing, software
10
// distributed under the License is distributed on an "AS IS" BASIS,
11
// WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
12
// See the License for the specific language governing permissions and
13
// limitations under the License.
14

15
#include "sparse/ops/cc/sddmm_launcher.h"
16

17
namespace sgk {
18

19
// Simple CPU kernel launcher.
20
void LaunchSddmm(const Eigen::ThreadPoolDevice &d, int m, int k, int n,
21
                 int nonzeros, const int *row_indices, const int *row_offsets,
22
                 const int *column_indices, const float *lhs_matrix,
23
                 const float *rhs_matrix, float *output_values) {
24
  for (int i = 0; i < m; ++i) {
25
    for (int j = row_offsets[i]; j < row_offsets[i + 1]; ++j) {
26
      int idx_n = column_indices[j];
27
      float accumulator = 0.0f;
28
      for (int l = 0; l < k; ++l) {
29
        accumulator += lhs_matrix[i * k + l] * rhs_matrix[idx_n * k + l];
30
      }
31
      output_values[j] = accumulator;
32
    }
33
  }
34
}
35

36
}  // namespace sgk
37