[INFRA-10] Convert Blitz++ code to Eigen Created: 17/Jul/14  Updated: 24/Jun/16  Resolved: 24/Jun/16

Status: Done
Project: Software Development Infrastructure
Component/s: None
Affects Version/s: None
Fix Version/s: None

Type: Story Priority: Major
Reporter: rhl Assignee: aritter
Resolution: Done Votes: 0
Labels: None
Remaining Estimate: Not Specified
Time Spent: Not Specified
Original Estimate: Not Specified

Issue Links:
Blocks
is blocked by INFRA-25 Convert Stella to use an explicit nam... Done
Relates
relates to INFRA-8 Convert basic image basic structure t... Done
Story Points: 8
Epic Link: Restructure STELLA

 Description   

This may or may not be a good idea, but I'd like to think about it.

  • Is Blitz++ still maintained?
  • Does it require lapack or blas?
  • Is it fast?
  • Does it play nicely with numpy?


 Comments   
Comment by rhl [ 23/Jul/14 ]

Given that STELLA uses Blitz++ as its primary image class we may need to move this up. One possibility would be to keep Blitz++ internally for now, and construct Blitz++ arrays from afw::Images. Note that afwImage supports a getArray() method to recover a strided array, and that Blitz++ can use such arrays in its ctors (http://ehc.ac/p/blitz/mailman/message/2406262/)

Comment by aritter [ 28/Feb/15 ]

The blitz-intensive functions FiberTrace::SlitFunc and FiberTrace::MkSlitFunc as well as all blitz functions from Math.h and Utils.h needed by them have been moved to MathBlitz.h and UtilsBlitz.h. Many Math.h functions as well as MpFit have been changed to use ndarrays instead of blitz arrays. The remaining source code is now free of blitz. Tests and Tasks have been changed to the new environment.

Comment by aritter [ 28/Feb/15 ]

The Issue is ready for review...

Comment by aritter [ 24/Jun/16 ]

Everything related to the Blitz++ library has been converted to Eigen or ndarrays.

Generated at Sat Feb 10 16:48:29 JST 2024 using Jira 8.3.4#803005-sha1:1f96e09b3c60279a408a2ae47be3c745f571388b.