Pop_Slime/crtr
1
0
Fork 0
Code Issues 6 Pull Requests Projects 1 Releases Activity

Optimize GC for group frame.

Browse Source
This commit is contained in:
PopSlime
2022-12-20 16:05:52 +08:00
parent 571320630b
commit 06d8012675
3 changed files with 30 additions and 24 deletions
Download Patch File Download Diff File Expand all files Collapse all files

30
Assets/Cryville/Common/Math/SquareMatrix.cs
View File

@@ -1,9 +1,13 @@
namespace Cryville.Common.Math {
using System;
namespace Cryville.Common.Math {
/// <summary>
/// Represents a square matrix.
/// </summary>
public class SquareMatrix {
readonly float[,] content;
readonly float[,] buffer;
readonly int[] refl;
/// <summary>
/// The size of the matrix.
/// </summary>
@@ -17,6 +21,8 @@
/// <param name="size">The size of the matrix.</param>
public SquareMatrix(int size) {
content = new float[size, size];
buffer = new float[size, size];
refl = new int[size];
Size = size;
}
/// <summary>
@@ -38,38 +44,36 @@
/// <returns>The column vector eliminated.</returns>
public ColumnVector<T> Eliminate<T>(ColumnVector<T> v, IVectorOperator<T> o) {
int s = Size;
float[,] d = (float[,])content.Clone();
int[] refl = new int[s];
for (int i = 0; i < s; i++)
refl[i] = i;
Array.Copy(content, buffer, Size * Size);
for (int i = 0; i < s; i++) refl[i] = i;
for (int r = 0; r < s; r++) {
for (int r0 = r; r0 < s; r0++)
if (d[refl[r0], r] != 0) {
if (buffer[refl[r0], r] != 0) {
refl[r] = r0;
refl[r0] = r;
break;
}
int or = refl[r];
float sf0 = d[or, r];
float sf0 = buffer[or, r];
for (int c0 = r; c0 < s; c0++)
d[or, c0] /= sf0;
buffer[or, c0] /= sf0;
v[or] = o.ScalarMultiply(1 / sf0, v[or]);
for (int r1 = r + 1; r1 < s; r1++) {
int or1 = refl[r1];
float sf1 = d[or1, r];
float sf1 = buffer[or1, r];
for (int c1 = r; c1 < s; c1++)
d[or1, c1] -= d[or, c1] * sf1;
buffer[or1, c1] -= buffer[or, c1] * sf1;
v[or1] = o.Add(v[or1], o.ScalarMultiply(-sf1, v[or]));
}
}
T[] res = new T[s];
ColumnVector<T> res = new ColumnVector<T>(s);
for (int r2 = s - 1; r2 >= 0; r2--) {
var v2 = v[refl[r2]];
for (int c2 = r2 + 1; c2 < s; c2++)
v2 = o.Add(v2, o.ScalarMultiply(-d[refl[r2], c2], res[refl[c2]]));
v2 = o.Add(v2, o.ScalarMultiply(-buffer[refl[r2], c2], res[refl[c2]]));
res[refl[r2]] = v2;
}
return new ColumnVector<T>(res);
return res;
}
/// <summary>
/// Creates a square matrix and fills it with polynomial coefficients.