Hi, > I've slightly modified the CG implementation to handle preconditionned > CG and unpreconditionned CG in the same routine (without adding any > extra allocation/copy), but I can't find a routine to test that I > haven't introduced a bug. Running > grep "cg_tag" tests/src/* > returns no result. > Do we have a test for the iterative solvers, or just examples?
unfortunately there is no native test in ViennaCL itself. Some other Vienna-libraries test their solvers against viennacl-dev, so this part of functionality is checked indirectly. The easiest way to verify proper operation 'by hand' is to check the residuals reported through 'solverbench'. Btw: How does your unified implementation behave performance-wise? I think the reason this was not unified is that a plain CG somewhere saves a memory transfer (norm() instead of inner_prod() or related). On the other hand, the pipelined CG version I'll push soon will require some extra tweaks for ViennaCL-types, so it's good to have a single 'generic' implementation for non-ViennaCL types. Best regards, Karli ------------------------------------------------------------------------------ "Accelerate Dev Cycles with Automated Cross-Browser Testing - For FREE Instantly run your Selenium tests across 300+ browser/OS combos. Get unparalleled scalability from the best Selenium testing platform available Simple to use. Nothing to install. Get started now for free." http://p.sf.net/sfu/SauceLabs _______________________________________________ ViennaCL-devel mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/viennacl-devel
