Patch Detail
Show a patch.
GET /api/1.1/patches/20261/?format=api
{ "id": 20261, "url": "https://patchwork.libcamera.org/api/1.1/patches/20261/?format=api", "web_url": "https://patchwork.libcamera.org/patch/20261/", "project": { "id": 1, "url": "https://patchwork.libcamera.org/api/1.1/projects/1/?format=api", "name": "libcamera", "link_name": "libcamera", "list_id": "libcamera_core", "list_email": "libcamera-devel@lists.libcamera.org", "web_url": "", "scm_url": "", "webscm_url": "" }, "msgid": "<20240611140207.520083-2-paul.elder@ideasonboard.com>", "date": "2024-06-11T14:02:05", "name": "[v7,1/3] ipa: libipa: Add Matrix class", "commit_ref": null, "pull_url": null, "state": "superseded", "archived": false, "hash": "b6190295c015dd95c2a4491195217fe29fa0da52", "submitter": { "id": 17, "url": "https://patchwork.libcamera.org/api/1.1/people/17/?format=api", "name": "Paul Elder", "email": "paul.elder@ideasonboard.com" }, "delegate": null, "mbox": "https://patchwork.libcamera.org/patch/20261/mbox/", "series": [ { "id": 4383, "url": "https://patchwork.libcamera.org/api/1.1/series/4383/?format=api", "web_url": "https://patchwork.libcamera.org/project/libcamera/list/?series=4383", "date": "2024-06-11T14:02:04", "name": "ipa: rkisp1: Add crosstalk algorithm", "version": 7, "mbox": "https://patchwork.libcamera.org/series/4383/mbox/" } ], "comments": "https://patchwork.libcamera.org/api/patches/20261/comments/", "check": "pending", "checks": "https://patchwork.libcamera.org/api/patches/20261/checks/", "tags": {}, "headers": { "Return-Path": "<libcamera-devel-bounces@lists.libcamera.org>", "X-Original-To": "parsemail@patchwork.libcamera.org", "Delivered-To": "parsemail@patchwork.libcamera.org", "Received": [ "from lancelot.ideasonboard.com (lancelot.ideasonboard.com\n\t[92.243.16.209])\n\tby patchwork.libcamera.org (Postfix) with ESMTPS id 34B54C31E9\n\tfor <parsemail@patchwork.libcamera.org>;\n\tTue, 11 Jun 2024 14:02:34 +0000 (UTC)", "from lancelot.ideasonboard.com (localhost [IPv6:::1])\n\tby lancelot.ideasonboard.com (Postfix) with ESMTP id DA6AB6546C;\n\tTue, 11 Jun 2024 16:02:33 +0200 (CEST)", "from perceval.ideasonboard.com (perceval.ideasonboard.com\n\t[213.167.242.64])\n\tby lancelot.ideasonboard.com (Postfix) with ESMTPS id 4261261A26\n\tfor <libcamera-devel@lists.libcamera.org>;\n\tTue, 11 Jun 2024 16:02:24 +0200 (CEST)", "from neptunite.hamster-moth.ts.net\n\t(h175-177-049-156.catv02.itscom.jp [175.177.49.156])\n\tby perceval.ideasonboard.com (Postfix) with ESMTPSA id 39888E0C;\n\tTue, 11 Jun 2024 16:02:07 +0200 (CEST)" ], "Authentication-Results": "lancelot.ideasonboard.com; dkim=pass (1024-bit key;\n\tunprotected) header.d=ideasonboard.com header.i=@ideasonboard.com\n\theader.b=\"SVRF4ACs\"; dkim-atps=neutral", "DKIM-Signature": "v=1; a=rsa-sha256; c=relaxed/simple; d=ideasonboard.com;\n\ts=mail; t=1718114531;\n\tbh=6SkWhaVo8yh2kxVL49xg198Jfjk9xEEAb6+rbG6pks8=;\n\th=From:To:Cc:Subject:Date:In-Reply-To:References:From;\n\tb=SVRF4ACsbbH4JyOv3uqn4iG8P6Q3qwGCd9u+uE+xOMGp+bWtqjNwQFej3Z9f0pd2G\n\trv972JsFbOFj0NKy9ST388vEA8gD3hac695P+ydpRh+6+g8FOB8A86BEkXnFz3xehe\n\tySTJKjOAMB/+vjBL9BGsX/NWcK4M4KO+gfjAeHwA=", "From": "Paul Elder <paul.elder@ideasonboard.com>", "To": "libcamera-devel@lists.libcamera.org", "Cc": "Paul Elder <paul.elder@ideasonboard.com>,\n\tStefan Klug <stefan.klug@ideasonboard.com>,\n\tKieran Bingham <kieran.bingham@ideasonboard.com>", "Subject": "[PATCH v7 1/3] ipa: libipa: Add Matrix class", "Date": "Tue, 11 Jun 2024 23:02:05 +0900", "Message-Id": "<20240611140207.520083-2-paul.elder@ideasonboard.com>", "X-Mailer": "git-send-email 2.39.2", "In-Reply-To": "<20240611140207.520083-1-paul.elder@ideasonboard.com>", "References": "<20240611140207.520083-1-paul.elder@ideasonboard.com>", "MIME-Version": "1.0", "Content-Transfer-Encoding": "8bit", "X-BeenThere": "libcamera-devel@lists.libcamera.org", "X-Mailman-Version": "2.1.29", "Precedence": "list", "List-Id": "<libcamera-devel.lists.libcamera.org>", "List-Unsubscribe": "<https://lists.libcamera.org/options/libcamera-devel>,\n\t<mailto:libcamera-devel-request@lists.libcamera.org?subject=unsubscribe>", "List-Archive": "<https://lists.libcamera.org/pipermail/libcamera-devel/>", "List-Post": "<mailto:libcamera-devel@lists.libcamera.org>", "List-Help": "<mailto:libcamera-devel-request@lists.libcamera.org?subject=help>", "List-Subscribe": "<https://lists.libcamera.org/listinfo/libcamera-devel>,\n\t<mailto:libcamera-devel-request@lists.libcamera.org?subject=subscribe>", "Errors-To": "libcamera-devel-bounces@lists.libcamera.org", "Sender": "\"libcamera-devel\" <libcamera-devel-bounces@lists.libcamera.org>" }, "content": "Add a class to represent a Matrix object and operations for adding\nmatrices, multipling a matrix by a scalar, and multiplying two matrices.\n\nSigned-off-by: Paul Elder <paul.elder@ideasonboard.com>\nReviewed-by: Stefan Klug <stefan.klug@ideasonboard.com>\nReviewed-by: Kieran Bingham <kieran.bingham@ideasonboard.com>\n\n---\nChanges in v7:\n- fix copyright and license\n\nChanges in v6:\n- fix doxygen\n\nChanges in v5:\n- add documentation\n\nChanges in v4:\n- remove stray semicolons\n- add operator<<\n- clean up/optimize constructor\n- replace get() and set() with operator[] (and a second [] can be used\n as operator[] returns a Span)\n\nChanges in v3:\n- fix template parameters of operator* to allow different types for the\n scalar multiplier and the matrix's number type\n- clear data in constructors\n- fix assert in constructor\n\nChanges v2:\n- make rows and columns into template arguments\n- initialize to identity matrix on construction\n- add getter and setter\n- change from struct to class\n- fix matrix multiplication\n- clean up unused includes\n- avoid dereferencing an absent std::optional\n---\n src/ipa/libipa/matrix.cpp | 123 ++++++++++++++++++++++++++\n src/ipa/libipa/matrix.h | 172 +++++++++++++++++++++++++++++++++++++\n src/ipa/libipa/meson.build | 2 +\n 3 files changed, 297 insertions(+)\n create mode 100644 src/ipa/libipa/matrix.cpp\n create mode 100644 src/ipa/libipa/matrix.h", "diff": "diff --git a/src/ipa/libipa/matrix.cpp b/src/ipa/libipa/matrix.cpp\nnew file mode 100644\nindex 000000000000..02091a0b4ce1\n--- /dev/null\n+++ b/src/ipa/libipa/matrix.cpp\n@@ -0,0 +1,123 @@\n+/* SPDX-License-Identifier: LGPL-2.1-or-later */\n+/*\n+ * Copyright (C) 2024, Paul Elder <paul.elder@ideasonboard.com>\n+ *\n+ * Matrix and related operations\n+ */\n+\n+#include \"matrix.h\"\n+\n+#include <libcamera/base/log.h>\n+\n+/**\n+ * \\file matrix.h\n+ * \\brief Matrix class\n+ */\n+\n+namespace libcamera {\n+\n+LOG_DEFINE_CATEGORY(Matrix)\n+\n+namespace ipa {\n+\n+/**\n+ * \\class Matrix\n+ * \\brief Matrix class\n+ * \\tparam T Type of numerical values to be stored in the matrix\n+ * \\tparam R Number of rows in the matrix\n+ * \\tparam C Number of columns in the matrix\n+ */\n+\n+/**\n+ * \\fn Matrix::Matrix()\n+ * \\brief Construct an identity matrix\n+ */\n+\n+/**\n+ * \\fn Matrix::Matrix(const std::vector<T> &data)\n+ * \\brief Construct matrix from supplied data\n+ * \\param data Data from which to construct a matrix\n+ *\n+ * \\a data is a one-dimensional vector and will be turned into a matrix in\n+ * row-major order. The size of \\a data must be equal to the product of the\n+ * number of rows and columns of the matrix (RxC).\n+ */\n+\n+/**\n+ * \\fn Matrix::readYaml\n+ * \\brief Populate the matrix with yaml data\n+ * \\param yaml Yaml data to populate the matrix with\n+ *\n+ * Any existing data in the matrix will be overwritten. The size of the data\n+ * read from \\a yaml must be equal to the product of the number of rows and\n+ * columns of the matrix (RxC).\n+ *\n+ * The yaml data is expected to be a list with elements of type T.\n+ *\n+ * \\return 0 on success, negative error code otherwise\n+ */\n+\n+/**\n+ * \\fn Matrix::toString\n+ * \\brief Assemble and return a string describing the matrix\n+ * \\return A string describing the matrix\n+ */\n+\n+/**\n+ * \\fn Span<const T, C> Matrix::operator[](size_t i) const\n+ * \\brief Index to a row in the matrix\n+ * \\param i Index of row to retrieve\n+ *\n+ * This operator[] returns a Span, which can then be indexed into again with\n+ * another operator[], allowing a convenient m[i][j] to access elements of the\n+ * matrix. Note that the lifetime of the Span returned by this first-level\n+ * operator[] is bound to that of the Matrix itself, so it is not recommended\n+ * to save the Span that is the result of this operator[].\n+ *\n+ * \\return Row \\a i from the matrix, as a Span\n+ */\n+\n+/**\n+ * \\fn Matrix::operator[](size_t i)\n+ * \\copydoc Matrix::operator[](size_t i) const\n+ */\n+\n+/**\n+ * \\fn Matrix::Matrix<U, R, C> operator*(T d, const Matrix<U, R, C> &m)\n+ * \\brief Scalar product\n+ * \\tparam T Type of the numerical scalar value\n+ * \\tparam U Type of numerical values in the matrix\n+ * \\tparam R Number of rows in the matrix\n+ * \\tparam C Number of columns in the matrix\n+ * \\param d Scalar\n+ * \\param m Matrix\n+ * \\return Product of scalar \\a d and matrix \\a m\n+ */\n+\n+/**\n+ * \\fn Matrix<T, R1, C2> operator*(const Matrix<T, R1, C1> &m1, const Matrix<T, R2, C2> &m2)\n+ * \\brief Matrix multiplication\n+ * \\tparam T Type of numerical values in the matrices\n+ * \\tparam R1 Number of rows in the first matrix\n+ * \\tparam C1 Number of columns in the first matrix\n+ * \\tparam R2 Number of rows in the second matrix\n+ * \\tparam C2 Number of columns in the second matrix\n+ * \\param m1 Multiplicand matrix\n+ * \\param m2 Multiplier matrix\n+ * \\return Matrix product of matrices \\a m1 and \\a m2\n+ */\n+\n+/**\n+ * \\fn Matrix<T, R, C> operator+(const Matrix<T, R, C> &m1, const Matrix<T, R, C> &m2)\n+ * \\brief Matrix addition\n+ * \\tparam T Type of numerical values in the matrices\n+ * \\tparam R Number of rows in the matrices\n+ * \\tparam C Number of columns in the matrices\n+ * \\param m1 Summand matrix\n+ * \\param m2 Summand matrix\n+ * \\return Matrix sum of matrices \\a m1 and \\a m2\n+ */\n+\n+} /* namespace ipa */\n+\n+} /* namespace libcamera */\ndiff --git a/src/ipa/libipa/matrix.h b/src/ipa/libipa/matrix.h\nnew file mode 100644\nindex 000000000000..90eaea03bd14\n--- /dev/null\n+++ b/src/ipa/libipa/matrix.h\n@@ -0,0 +1,172 @@\n+/* SPDX-License-Identifier: LGPL-2.1-or-later */\n+/*\n+ * Copyright (C) 2024, Paul Elder <paul.elder@ideasonboard.com>\n+ *\n+ * Matrix and related operations\n+ */\n+#pragma once\n+\n+#include <algorithm>\n+#include <cmath>\n+#include <sstream>\n+#include <vector>\n+\n+#include <libcamera/base/log.h>\n+#include <libcamera/base/span.h>\n+\n+#include \"libcamera/internal/yaml_parser.h\"\n+\n+namespace libcamera {\n+\n+LOG_DECLARE_CATEGORY(Matrix)\n+\n+namespace ipa {\n+\n+#ifndef __DOXYGEN__\n+template<typename T, unsigned int R, unsigned int C,\n+\t std::enable_if_t<std::is_arithmetic_v<T>> * = nullptr>\n+#else\n+template<typename T, unsigned int R, unsigned int C>\n+#endif /* __DOXYGEN__ */\n+class Matrix\n+{\n+public:\n+\tMatrix()\n+\t\t: data_(R * C, static_cast<T>(false))\n+\t{\n+\t\tfor (size_t i = 0; i < std::min(R, C); i++)\n+\t\t\t(*this)[i][i] = static_cast<T>(true);\n+\t}\n+\n+\tMatrix(const std::vector<T> &data)\n+\t{\n+\t\tASSERT(data.size() == R * C);\n+\n+\t\tdata_.clear();\n+\t\tfor (const T &x : data)\n+\t\t\tdata_.push_back(x);\n+\t}\n+\n+\t~Matrix() = default;\n+\n+\tint readYaml(const libcamera::YamlObject &yaml)\n+\t{\n+\t\tif (yaml.size() != R * C) {\n+\t\t\tLOG(Matrix, Error)\n+\t\t\t\t<< \"Wrong number of values in matrix: expected \"\n+\t\t\t\t<< R * C << \", got \" << yaml.size();\n+\t\t\treturn -EINVAL;\n+\t\t}\n+\n+\t\tunsigned int i = 0;\n+\t\tfor (const auto &x : yaml.asList()) {\n+\t\t\tauto value = x.get<T>();\n+\t\t\tif (!value) {\n+\t\t\t\tLOG(Matrix, Error) << \"Failed to read matrix value\";\n+\t\t\t\treturn -EINVAL;\n+\t\t\t}\n+\n+\t\t\tdata_[i++] = *value;\n+\t\t}\n+\n+\t\treturn 0;\n+\t}\n+\n+\tconst std::string toString() const\n+\t{\n+\t\tstd::stringstream out;\n+\n+\t\tout << \"Matrix { \";\n+\t\tfor (unsigned int i = 0; i < R; i++) {\n+\t\t\tout << \"[ \";\n+\t\t\tfor (unsigned int j = 0; j < C; j++) {\n+\t\t\t\tout << (*this)[i][j];\n+\t\t\t\tout << ((j + 1 < C) ? \", \" : \" \");\n+\t\t\t}\n+\t\t\tout << ((i + 1 < R) ? \"], \" : \"]\");\n+\t\t}\n+\t\tout << \" }\";\n+\n+\t\treturn out.str();\n+\t}\n+\n+\tSpan<const T, C> operator[](size_t i) const\n+\t{\n+\t\treturn Span<const T, C>{ &data_.data()[i * C], C };\n+\t}\n+\n+\tSpan<T, C> operator[](size_t i)\n+\t{\n+\t\treturn Span<T, C>{ &data_.data()[i * C], C };\n+\t}\n+\n+private:\n+\tstd::vector<T> data_;\n+};\n+\n+#ifndef __DOXYGEN__\n+template<typename T, typename U, unsigned int R, unsigned int C,\n+\t std::enable_if_t<std::is_arithmetic_v<T> && std::is_arithmetic_v<U>> * = nullptr>\n+#endif /* __DOXYGEN__ */\n+Matrix<U, R, C> operator*(T d, const Matrix<U, R, C> &m)\n+{\n+\tMatrix<U, R, C> result;\n+\n+\tfor (unsigned int i = 0; i < R; i++)\n+\t\tfor (unsigned int j = 0; j < C; j++)\n+\t\t\tresult[i][j] = d * m[i][j];\n+\n+\treturn result;\n+}\n+\n+#ifndef __DOXYGEN__\n+template<typename T,\n+\t unsigned int R1, unsigned int C1,\n+\t unsigned int R2, unsigned int C2,\n+\t std::enable_if_t<std::is_arithmetic_v<T> && C1 == R2> * = nullptr>\n+#endif /* __DOXYGEN__ */\n+Matrix<T, R1, C2> operator*(const Matrix<T, R1, C1> &m1, const Matrix<T, R2, C2> &m2)\n+{\n+\tMatrix<T, R1, C2> result;\n+\n+\tfor (unsigned int i = 0; i < R1; i++) {\n+\t\tfor (unsigned int j = 0; j < C2; j++) {\n+\t\t\tT sum = 0;\n+\n+\t\t\tfor (unsigned int k = 0; k < C1; k++)\n+\t\t\t\tsum += m1[i][k] * m2[k][j];\n+\n+\t\t\tresult[i][j] = sum;\n+\t\t}\n+\t}\n+\n+\treturn result;\n+}\n+\n+#ifndef __DOXYGEN__\n+template<typename T, unsigned int R, unsigned int C,\n+\t std::enable_if_t<std::is_arithmetic_v<T>> * = nullptr>\n+#endif /* __DOXYGEN__ */\n+Matrix<T, R, C> operator+(const Matrix<T, R, C> &m1, const Matrix<T, R, C> &m2)\n+{\n+\tMatrix<T, R, C> result;\n+\n+\tfor (unsigned int i = 0; i < R; i++)\n+\t\tfor (unsigned int j = 0; j < C; j++)\n+\t\t\tresult[i][j] = m1[i][j] + m2[i][j];\n+\n+\treturn result;\n+}\n+\n+} /* namespace ipa */\n+\n+#ifndef __DOXYGEN__\n+template<typename T, unsigned int R, unsigned int C>\n+std::ostream &operator<<(std::ostream &out, const ipa::Matrix<T, R, C> &m)\n+{\n+\tout << m.toString();\n+\treturn out;\n+}\n+#endif /* __DOXYGEN__ */\n+\n+} /* namespace libcamera */\ndiff --git a/src/ipa/libipa/meson.build b/src/ipa/libipa/meson.build\nindex b663afc7d9fe..067d0d273e0a 100644\n--- a/src/ipa/libipa/meson.build\n+++ b/src/ipa/libipa/meson.build\n@@ -7,6 +7,7 @@ libipa_headers = files([\n 'exposure_mode_helper.h',\n 'fc_queue.h',\n 'histogram.h',\n+ 'matrix.h',\n 'module.h',\n ])\n \n@@ -17,6 +18,7 @@ libipa_sources = files([\n 'exposure_mode_helper.cpp',\n 'fc_queue.cpp',\n 'histogram.cpp',\n+ 'matrix.cpp',\n 'module.cpp',\n ])\n \n", "prefixes": [ "v7", "1/3" ] }