Skip to content
Projects
Groups
Snippets
Help
Loading...
Help
Support
Keyboard shortcuts
?
Submit feedback
Contribute to GitLab
Sign in / Register
Toggle navigation
L
LGPU FPP
Project overview
Project overview
Details
Activity
Releases
Repository
Repository
Files
Commits
Branches
Tags
Contributors
Graph
Compare
Issues
0
Issues
0
List
Boards
Labels
Milestones
Merge Requests
0
Merge Requests
0
CI / CD
CI / CD
Pipelines
Jobs
Schedules
Analytics
Analytics
CI / CD
Repository
Value Stream
Wiki
Wiki
Snippets
Snippets
Members
Members
Collapse sidebar
Close sidebar
Activity
Graph
Create a new issue
Jobs
Commits
Issue Boards
Open sidebar
Fernando Pérez Panadero
LGPU FPP
Commits
6503bdf3
Commit
6503bdf3
authored
1 year ago
by
Fernando Pérez Panadero
Browse files
Options
Browse Files
Download
Email Patches
Plain Diff
Change of csw done in CG
parent
8c20c084
Changes
1
Hide whitespace changes
Inline
Side-by-side
Showing
1 changed file
with
2 additions
and
2 deletions
+2
-2
src/Solvers/CG.jl
src/Solvers/CG.jl
+2
-2
No files found.
src/Solvers/CG.jl
View file @
6503bdf3
...
...
@@ -34,7 +34,7 @@ return sum(sumf)
end
function
CG
!
(
si
,
U
,
A
,
dpar
::
DiracParam
,
lp
::
SpaceParm
,
dws
::
DiracWorkspace
{
T
},
ymws
::
YMworkspace
,
maxiter
::
Int64
=
10
,
tol
=
1.0
)
where
{
T
}
function
CG
!
(
si
,
U
,
A
,
dpar
::
DiracParam
,
lp
::
SpaceParm
,
dws
::
DiracWorkspace
{
T
},
maxiter
::
Int64
=
10
,
tol
=
1.0
)
where
{
T
}
dws
.
sr
.=
si
dws
.
sp
.=
si
...
...
@@ -49,7 +49,7 @@ function CG!(si, U, A, dpar::DiracParam, lp::SpaceParm, dws::DiracWorkspace{T},
niter
=
0
for
i
in
1
:
maxiter
A
(
dws
.
sAp
,
U
,
dws
.
sp
,
dpar
,
dws
.
st
,
lp
,
ymws
)
A
(
dws
.
sAp
,
U
,
dws
.
sp
,
dpar
,
dws
.
st
,
lp
)
prod
=
field_dot
(
dws
.
sp
,
dws
.
sAp
,
sumf
,
lp
)
#prod = CUDA.mapreduce(x -> dot(x[1],x[2]), +, zip(dws.sp, dws.sAp))
...
...
This diff is collapsed.
Click to expand it.
Write
Preview
Markdown
is supported
0%
Try again
or
attach a new file
Attach a file
Cancel
You are about to add
0
people
to the discussion. Proceed with caution.
Finish editing this message first!
Cancel
Please
register
or
sign in
to comment