Can chowiluviennacl do ilu0? I need to solve a tri-diagonal system directly. If I apply the PCILU, I will obtain the exact solution with preonly + pcilu. However, the preonly + chowiluviennacl will not provide the exact solution. Any option keys to set the CHOWILUVIENNACL filling level or dropping off tolerance like the standard ilu?
Thank you. Best, Xiangdong On Tue, Jan 14, 2020 at 10:05 PM Matthew Knepley <[email protected]> wrote: > On Tue, Jan 14, 2020 at 9:56 PM Xiangdong <[email protected]> wrote: > >> Dear Developers, >> >> I have a quick question about the chowiluviennacl. When I tried to use >> it, I found that it only works for np=1, not np>1. However, in the >> description of chowiluviennacl.cxx, it says "the ViennaCL Chow-Patel >> parallel ILU preconditioner". >> > > By parallel, this means shared memory parallelism on the GPU. > > >> I am wondering whether I am using it correctly. Does chowiluviennacl work >> for np>1? >> > > I do not believe so. I do not see why it could not be extended, but that > would mean writing some more code. > > Thanks, > > Matt > > >> In addition, are there option keys for the chowiluviennacl one can try? >> Thank you. >> >> Best, >> Xiangdong >> > > > -- > What most experimenters take for granted before they begin their > experiments is infinitely more interesting than any results to which their > experiments lead. > -- Norbert Wiener > > https://www.cse.buffalo.edu/~knepley/ > <http://www.cse.buffalo.edu/~knepley/> >
