[GRASS-dev] calculating inverse matrix

i.pan.sharpen contains python code to do an inverse pca. It uses some of the new r.mapcalc capabilities to do parallel processing on multi-core processors. So it’s pretty fast.

Michael


C. Michael Barton
Director, Center for Social Dynamics & Complexity
Professor of Anthropology, School of Human Evolution & Social Change
Arizona State University

voice: 480-965-6262 (SHESC), 480-727-9746 (CSDC)
fax: 480-965-7671 (SHESC), 480-727-0709 (CSDC)

www: http://www.public.asu.edu/~cmbarton, http://csdc.asu.edu

On Mar 26, 2013, at 4:24 AM, <grass-dev-request@lists.osgeo.org>
wrote:

From: Martin Landa <landa.martin@gmail.com>

Subject: [GRASS-dev] calculating inverse matrix

Date: March 26, 2013 4:20:08 AM MST

To: GRASS developers list <grass-dev@lists.osgeo.org>

Hi all,

I found several methods in GRASS code which compute inverse matrix

  1. v.generalize

http://trac.osgeo.org/grass/browser/grass/trunk/vector/v.generalize/matrix.c#L134

  1. i.orto.photo

http://trac.osgeo.org/grass/browser/grass/trunk/imagery/i.ortho.photo/lib/m_inverse.c#L16

  1. GMath lib (depends on lapack/blas)

http://trac.osgeo.org/grass/browser/grass/trunk/lib/gmath/la.c#L566

It would be nice to consolidate it including data structures (modules
noted above define their own data structures for matrix). In other
worlds to incorporate v.generalize/i.orto.photo code into GMath lib.
Any pointers here?

Thanks, Martin


Martin Landa <landa.martin gmail.com> * http://geo.fsv.cvut.cz/~landa

Does that require compiling with OpenMP?
Doug

On Tue, Mar 26, 2013 at 1:27 PM, Michael Barton <Michael.Barton@asu.edu> wrote:

i.pan.sharpen contains python code to do an inverse pca. It uses some of the new r.mapcalc capabilities to do parallel processing on multi-core processors. So it’s pretty fast.

Michael


C. Michael Barton
Director, Center for Social Dynamics & Complexity
Professor of Anthropology, School of Human Evolution & Social Change
Arizona State University

voice: 480-965-6262 (SHESC), 480-727-9746 (CSDC)
fax: 480-965-7671 (SHESC), 480-727-0709 (CSDC)

www: http://www.public.asu.edu/~cmbarton, http://csdc.asu.edu

On Mar 26, 2013, at 4:24 AM, <grass-dev-request@lists.osgeo.org>
wrote:

From: Martin Landa <landa.martin@gmail.com>

Subject: [GRASS-dev] calculating inverse matrix

Date: March 26, 2013 4:20:08 AM MST

To: GRASS developers list <grass-dev@lists.osgeo.org>

Hi all,

I found several methods in GRASS code which compute inverse matrix

  1. v.generalize

http://trac.osgeo.org/grass/browser/grass/trunk/vector/v.generalize/matrix.c#L134

  1. i.orto.photo

http://trac.osgeo.org/grass/browser/grass/trunk/imagery/i.ortho.photo/lib/m_inverse.c#L16

  1. GMath lib (depends on lapack/blas)

http://trac.osgeo.org/grass/browser/grass/trunk/lib/gmath/la.c#L566

It would be nice to consolidate it including data structures (modules
noted above define their own data structures for matrix). In other
worlds to incorporate v.generalize/i.orto.photo code into GMath lib.
Any pointers here?

Thanks, Martin


Martin Landa <landa.martin gmail.com> * http://geo.fsv.cvut.cz/~landa


grass-dev mailing list
grass-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-dev

Doug Newcomb
USFWS
Raleigh, NC
919-856-4520 ext. 14 doug_newcomb@fws.gov

The opinions I express are my own and are not representative of the official policy of the U.S.Fish and Wildlife Service or Dept. of the Interior. Life is too short for undocumented, proprietary data formats.

No it does not. It uses new grass python libraries to initiate multiple processes on different cores.

Michael


C. Michael Barton
Director, Center for Social Dynamics & Complexity
Professor of Anthropology, School of Human Evolution & Social Change
Arizona State University
Tempe, AZ 85287-2402
USA

voice: 480-965-6262 (SHESC), 480-727-9746 (CSDC)
fax: 480-965-7671(SHESC), 480-727-0709 (CSDC)
www: http://csdc.asu.edu, http://shesc.asu.edu
http://www.public.asu.edu/~cmbarton

On Mar 26, 2013, at 10:41 AM, “Newcomb, Doug” <doug_newcomb@fws.gov>
wrote:

Does that require compiling with OpenMP?
Doug

On Tue, Mar 26, 2013 at 1:27 PM, Michael Barton <Michael.Barton@asu.edu> wrote:

i.pan.sharpen contains python code to do an inverse pca. It uses some of the new r.mapcalc capabilities to do parallel processing on multi-core processors. So it’s pretty fast.

Michael


C. Michael Barton
Director, Center for Social Dynamics & Complexity
Professor of Anthropology, School of Human Evolution & Social Change
Arizona State University

voice: 480-965-6262 (SHESC), 480-727-9746 (CSDC)
fax: 480-965-7671 (SHESC), 480-727-0709 (CSDC)

www: http://www.public.asu.edu/~cmbarton, http://csdc.asu.edu

On Mar 26, 2013, at 4:24 AM, <grass-dev-request@lists.osgeo.org>
wrote:

From: Martin Landa <landa.martin@gmail.com>

Subject: [GRASS-dev] calculating inverse matrix

Date: March 26, 2013 4:20:08 AM MST

To: GRASS developers list <grass-dev@lists.osgeo.org>

Hi all,

I found several methods in GRASS code which compute inverse matrix

  1. v.generalize

http://trac.osgeo.org/grass/browser/grass/trunk/vector/v.generalize/matrix.c#L134

  1. i.orto.photo

http://trac.osgeo.org/grass/browser/grass/trunk/imagery/i.ortho.photo/lib/m_inverse.c#L16

  1. GMath lib (depends on lapack/blas)

http://trac.osgeo.org/grass/browser/grass/trunk/lib/gmath/la.c#L566

It would be nice to consolidate it including data structures (modules
noted above define their own data structures for matrix). In other
worlds to incorporate v.generalize/i.orto.photo code into GMath lib.
Any pointers here?

Thanks, Martin


Martin Landa <landa.martin gmail.com> * http://geo.fsv.cvut.cz/~landa


grass-dev mailing list
grass-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-dev

Doug Newcomb
USFWS
Raleigh, NC
919-856-4520 ext. 14 doug_newcomb@fws.gov

The opinions I express are my own and are not representative of the official policy of the U.S.Fish and Wildlife Service or Dept. of the Interior. Life is too short for undocumented, proprietary data formats.