QR decomposition: Difference between revisions
m
→{{header|Wren}}: Minor tidy
m (Minor refactoring.) |
m (→{{header|Wren}}: Minor tidy) |
||
| (2 intermediate revisions by one other user not shown) | |||
Line 2,162:
vectorE.set_entry(i, 0, ( i == k ) ? 1 : 0);
}
vectorE = vectorE.scalarMultiply(magnitudeX).add(vectorX).unit();
versions_of_Q.emplace_back(householder_factor(vectorE));
z = versions_of_Q[k].multiply(z1);
Line 3,316 ⟶ 3,315:
for ( int k = 0; k < columnCount && k < rowCount - 1; k++ ) {
Matrix vectorE = new Matrix(rowCount,
z1 = z.minor(k);
Matrix vectorX = z1.column(k);
Line 3,327 ⟶ 3,326:
vectorE.setEntry(i, 0, ( i == k ) ? 1 : 0);
}
vectorE = vectorE.scalarMultiply(magnitudeX).add(vectorX).unit();
versionsOfQ.add(householderFactor(vectorE));
z = versionsOfQ.get(k).multiply(z1);
Line 3,349 ⟶ 3,347:
final int size = aVector.size();
for ( int i = 0; i < size; i++ ) {
for ( int j = 0; j < size; j++ ) {
}
}
for ( int i = 0; i < size; i++ ) {
result.setEntry(i, i,
}
return
}
Line 3,412 ⟶ 3,410:
}
public Matrix(int
this( new double[
}
Line 3,470 ⟶ 3,468:
}
public Matrix column(int
Matrix result = new Matrix(rowCount, 1);
for ( int i = 0; i < rowCount; i++ ) {
result.setEntry(i, 0, data[i][
}
return result;
Line 3,484 ⟶ 3,482:
Matrix result = new Matrix(rowCount, columnCount);
for ( int i = 0; i <
result.data[i][0] = data[i][0] * aValue;
}
Line 3,497 ⟶ 3,495:
final double magnitude = magnitude();
Matrix result = new Matrix(rowCount, columnCount);
for ( int i = 0; i <
result.data[i][0] = data[i][0] / magnitude;
}
Line 3,519 ⟶ 3,517:
throw new IllegalArgumentException("Incompatible matrix dimensions.");
}
return
}
Line 3,533 ⟶ 3,531:
}
public double getEntry(int aRow, int
return data[aRow][
}
public void setEntry(int aRow, int
data[aRow][
}
Line 5,444 ⟶ 5,442:
{{libheader|Wren-matrix}}
{{libheader|Wren-fmt}}
<syntaxhighlight lang="
import "./fmt" for Fmt
var minor = Fn.new { |x, d|
| |||