Skip to content

Commit

Permalink
Remove extra function parameter for posv (#24)
Browse files Browse the repository at this point in the history
  • Loading branch information
jmh530 committed Dec 8, 2023
1 parent f072a0a commit 1f5eb8c
Showing 1 changed file with 6 additions and 11 deletions.
17 changes: 6 additions & 11 deletions source/mir/lapack.d
Expand Up @@ -2275,21 +2275,19 @@ template posv(T)
size_t posv(
char uplo,
Slice!(T*, 2, Canonical) a,
Slice!(T*, 2, Canonical) b,
Slice!(T*, 2, Canonical) x
Slice!(T*, 2, Canonical) b
)
in {
assert(uplo == 'U' || uplo == 'L');
auto n = a.length!0;
auto nrhs = x.length;
auto nrhs = b.length!0;
assert(a.length!1 == n);
assert(x.length!1 == n);
assert(b.length!1 == n);
assert(b.length!0 == nrhs);
}
do {
lapackint n = cast(lapackint) a.length!0;
lapackint nrhs = cast(lapackint) x.length;
lapackint nrhs = cast(lapackint) b.length!0;
lapackint lda = cast(lapackint) a._stride.max(1);
lapackint ldb = cast(lapackint) b._stride.max(1);
lapackint info;
Expand All @@ -2302,19 +2300,16 @@ template posv(T)
size_t posv(
char uplo,
Slice!(T*, 2, Canonical) a,
Slice!(T*) b,
Slice!(T*) x
Slice!(T*) b
)
in {
assert(uplo == 'U' || uplo == 'L');
auto n = a.length!0;
assert(a.length!1 == n);
assert(x.length == n);
assert(a.length!0 == a.length!1);
assert(b.length == n);
}
do {
import mir.ndslice.topology: canonical;
return posv(uplo, a, b.sliced(1, b.length).canonical, x.sliced(1, x.length).canonical);
return posv(uplo, a, b.sliced(1, b.length).canonical);
}
}

Expand Down

0 comments on commit 1f5eb8c

Please sign in to comment.