sunwen/UR
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

single bug fix

Browse Source
This commit is contained in:
kradchen
2023-06-02 13:44:14 +08:00
parent c1eeadf72a
commit 9908202e9f
4 changed files with 21 additions and 10 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
src/transmissionReconstruction/reconstruction/buildMatrix/buildMatrix.cpp
View File

@@ -220,9 +220,7 @@ namespace Recon
uint *temp = new uint[path.getDataSize()]{0};
//uint32(path)
std::copy(path.getData(),path.getData()+path.getDataSize(),temp);
std::copy(temp,temp+path.getDataSize(),path.getData());
delete [] temp;
Aurora::round(path);
result.weighting = traceLineResult.ds;
result.path = correctPathToImgDimensions(path, dims);
result.pathLen = result.path.getDimSize(0);

9
src/transmissionReconstruction/reconstruction/solvingEquationSystem/TVAL/TVAL.cpp
View File

@@ -170,6 +170,10 @@ namespace Recon
double * csrValues;
sparse_index_base_t index;
mkl_sparse_d_export_csr(csrA, &index, &n_rows, &n_cols, &rows_start, &rows_end, &col_indx, &csrValues);
mkl_sparse_destroy(A);
delete [] xIdxs;
delete [] yIdxs;
int *row_idx = new int[n_rows+1];
std::copy(rows_start,rows_start+n_rows,row_idx);
row_idx[n_rows] = rows_end[n_rows-1];
@@ -178,14 +182,9 @@ namespace Recon
size_t bDims[3]={(size_t)b.getDimSize(0),(size_t)b.getDimSize(1),(size_t)b.getDimSize(2)};
size_t rdims[3] = {(size_t)dims[0], (size_t)dims[1], (size_t)dims[2]};
bool pagelocked = false;
auto result = TVALGPU(row_idx, col_indx, csrValues, M.getM(), M.getN(), nz, bData, bDims, rdims, opt, device, false);
mkl_sparse_destroy(A);
mkl_sparse_destroy(csrA);
delete [] row_idx;
delete [] xIdxs;
delete [] yIdxs;
delete [] bData;
return Aurora::Matrix::fromRawData(result.data, result.dims[0],result.dims[1],result.dims[2]);
// return Aurora::Matrix();