Skip to content
GitLab
  • Menu
Projects Groups Snippets
  • /
  • Help
    • Help
    • Support
    • Community forum
    • Submit feedback
    • Contribute to GitLab
    • Switch to GitLab Next
  • Sign in / Register
  • eigen eigen
  • Project information
    • Project information
    • Activity
    • Labels
    • Members
  • Repository
    • Repository
    • Files
    • Commits
    • Branches
    • Tags
    • Contributors
    • Graph
    • Compare
    • Locked Files
  • Issues 662
    • Issues 662
    • List
    • Boards
    • Service Desk
    • Milestones
    • Iterations
    • Requirements
  • Custom issue tracker
    • Custom issue tracker
  • Merge requests 19
    • Merge requests 19
  • CI/CD
    • CI/CD
    • Pipelines
    • Jobs
    • Schedules
    • Test Cases
  • Deployments
    • Deployments
    • Environments
    • Releases
  • Packages & Registries
    • Packages & Registries
    • Package Registry
    • Container Registry
    • Infrastructure Registry
  • Monitor
    • Monitor
    • Incidents
  • Analytics
    • Analytics
    • Value stream
    • CI/CD
    • Code review
    • Insights
    • Issue
    • Repository
  • Snippets
    • Snippets
  • Activity
  • Graph
  • Create a new issue
  • Jobs
  • Commits
  • Issue Boards
Collapse sidebar
  • libeigenlibeigen
  • eigeneigen
  • Issues
  • #371
Closed
Open
Created Dec 04, 2019 by Eigen Bugzilla@eigenbzReporter

Automatically transforming a*b.transpose()*c into a*b.dot(c)

Submitted by Rhys Ulerich

Assigned to Nobody

Link to original bugzilla bug (#371)
Version: 3.0

Description

Awhile back (http://listengine.tuxfamily.org/lists.tuxfamily.org/eigen/2009/09/msg00280.html)there was some discussion of automatically transforming a*b.transpose()c into the equivalent (and faster) ab.dot(c).

A quick test 3.0.3 shows the latter expression being appreciably faster than the former at -O3 using the (admittedly lousy) benchmark

Eigen::Vector3d a, b, c, d;  

for (unsigned long i = 0; i < 1000000; ++i) {  
    grvy_timer_begin("Outer");  
    a.setRandom();  
    b.setRandom();  
    c.setRandom();  
    d = a*b.transpose()*c;  
    grvy_timer_end("Outer");  

    grvy_timer_begin("Dot");  
    a.setRandom();  
    b.setRandom();  
    c.setRandom();  
    d = a*b.dot(c);  
    grvy_timer_end("Dot");  
}  

Very, very low priority. Just wanted to capture it as an issue.

  • Rhys

Depends on

#99 (closed)

Edited Dec 05, 2019 by Eigen Bugzilla
Assignee
Assign to
Time tracking