Displaying 20 results from an estimated 100 matches similar to: "Calling Fortran from C++"
2005 Dec 31
3
Problems with updating R-packages
Dear Helpers,
I am new in Linux and R, trying to update packages I need in my work,
but repeatedly I failed for reasons I do not understand.
First, I use:
> update.packages()
Then, I choose nearest mirror and accept update of 'Design' package:
> Design :
> Version 2.0-9 installed in /usr/lib/R/site-library
> Version 2.0-12 available at http://cran.r-mirror.de
> Update
2009 Feb 06
3
Sweave and backslashes
Hello Everyone,
I am an avid Sweave user and I am trying to pretty print floating point
numbers for latex output. For example in my document, I would like:
4.2\cdot 10^-{8}
Instead of:
4.2e-08
The Hmisc package has a nice function for doing this- but Hmisc has a ton of
dependencies and has proved very unportable to the various machines I work
on. So, I set out to write my own function that
2008 Oct 07
11
Creating GUIs for R
Hi all,
I have looked around for help on creating GUIs for R, but haven't found
anything. I would be interested in any advice or webpages that have
information on the best language, tutorials etc. for creating simple GUIs.
Mainly I want to do this as a heuristic exercise.
Thanks for any help.
Wade Wall
[[alternative HTML version deleted]]
2011 Nov 10
2
performance of adaptIntegrate vs. integrate
Dear list,
[cross-posting from Stack Overflow where this question has remained
unanswered for two weeks]
I'd like to perform a numerical integration in one dimension,
I = int_a^b f(x) dx
where the integrand f: x in IR -> f(x) in IR^p is vector-valued.
integrate() only allows scalar integrands, thus I would need to call
it many (p=200 typically) times, which sounds suboptimal. The
2008 Aug 12
1
X11 plotting error
Hey everyone,
I'm part of a team working at NOAA and we have developed a swan script that
analyzes data and generates plots for the forecasters. We are ready to put
the script on the cron and call it a summer but a technical issue has popped
up which has us puzzled.
When we run the script as an unprivileged user, it bails when png() gets
called to open a device for the output. The error is
2009 Nov 25
2
XML package example code?
I'm interested in parsing an html page. I should use XML, right? Could
you somebody show me some example code? Is there a tutorial for this
package?
2008 Aug 03
1
Convert date to decimal days
Hello all,
I have a quick question about formatting date strings..
I am currently debugging some Matlab code someone else wrote and since it is
so bad that I have to go through it line by line I figured that I would just
rewrite the thing in R.
The code produces plots of wave spectra with decimal days since the Epoch as
the x axis and wave period as the Y axis. I am able to convert the date
2008 Oct 14
3
gWidgets install
Hi,
I wanted to install the gWidgets package and ran the following code:
------------------------------------------
> install.packages("gWidgets")
trying URL 'http://www.ibiblio.org/pub/languages/R/CRAN/bin/windows/contrib/2.7/gWidgets_0.0-30.zip'
Content type 'application/zip' length 1354268 bytes (1.3 Mb)
opened URL
downloaded 1.3 Mb
package 'gWidgets'
2009 Jul 19
4
space in column name
I read a table from Microsoft Access using RODBC. Some of the variables had
a name with a space in it.
R has no problem with it but I do.
I cannot find out how to specify the space
names(alltime)
[1] "ID" "LVL7" "Ref Pv No" "Ref Pv Name" "DOS"
"Pt Last Name" "Pt First Name" "MRN"
2009 Sep 23
4
save txt file
Hi, I have 2 questions:
Question 1:
I define 2 variables: "a", "b":
> a<-rbinom(4,10,0.8)
output:
[1] 9 7 8 8
> b<-rbinom(2,6,0.7)
output:
[1] 4 5
if I write:
> write.table(a, file = "filename", etc. etc. .... )
it save only the values of variable "a".
There is a way to save in a .txt file the values "a" and "b" as
2009 Sep 30
5
Rounding error in seq(...)
Hi,
Today I was flabbergasted to see something that looks like a rounding
error in the very basic seq function in R.
> a = seq(0.1,0.9,by=0.1)
> a
[1] 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9
> a[1] == 0.1
[1] TRUE
> a[2] == 0.2
[1] TRUE
> a[3] == 0.3
[1] FALSE
It turns out that the alternative
> a = (1:9)/10
works just fine. Are there any good guides out there on how to deal
2010 Jan 03
3
F77_CALL, F77_NAME definition
I give up. Maybe it is my search (Windows) but I cannot seem to find the definition of the F77_CALL or F77_NAME macros. Either there are too many matches or the search just doesn't find it. For example where is the source for:
F77_CALL(dpotri)
?
Thank you.
Kevin
2008 Oct 11
6
Sweave-LaTEX question
Hi:
I am working on a publication and I have heard about LaTEX but I haven't actually tried to learn about it until today. I've found a few examples but I can't actually make them work properly. I have a couple of questions:
Does LATEX have to be installed on your computer? How does the xtable package and Sweave work together? How can I make the code below create a table as pdf file?
2014 Dec 20
2
Unexplained difference between results of dppsv and dpotri LAPACK routines
Dear R contributors,
Considering the following sample C code, that illustrates two possible
uses of a Cholesky decomp for inverting a matrix, equally valid at
least in theory:
SEXP test() {
int d = 2;
int info = 0;
double mat[4] = {2.5, 0.4, 0.4, 1.7};
double id[4] = {1.0, 0.0, 0.0, 1.0};
double lmat[3];
F77_CALL(dpotrf)("L", &d, mat, &d, &info);
lmat[0] = mat[0];
lmat[1]
2020 Jul 15
2
Openblas?
On 2020-07-15 14:36, Dirk Eddelbuettel wrote:
>
> G?ran,
>
> This is not an easy email to reply to because it _contains nothing
> reproducible_.
Thanks Dirk,
Sorry about that, but my real question was (see below): "Is the problem
that openblas uses C versions of blas?" That is, do I need to change
F77_CALL(name)(...);
to
cblas_name(...);
everywhere? And if so, is
2017 Feb 09
3
Ancient C /Fortran code linpack error
In my package 'glmmML' I'm using old C code and linpack in the
optimizing procedure. Specifically, one part of the code looks like this:
F77_CALL(dpoco)(*hessian, &bdim, &bdim, &rcond, work, info);
if (*info == 0){
F77_CALL(dpodi)(*hessian, &bdim, &bdim, det, &job);
........
This usually works OK, but with an ill-conditioned data
2008 Mar 27
1
Cannot update packages on F8
Dear All,
I have just updated R to the version 2.6.2 on F8 (with the official F8
rpm). However, when running as root the following command:
update.packages(checkBuilt=T)
I get a bunch of errors like the ones below. Any ideas?
Thanks in advance,
Paul
-----------------------------------------------
* Installing *source* package 'nlme' ...
** libs
gcc -m32 -std=gnu99 -I/usr/include/R
2012 Mar 06
2
Calling FORTRAN function from R issue?
Hello,
I am trying to call the BLAS Level1 function zdotc from R via
a .C call like this:
#include "R.h"
#include "R_ext/BLAS.h"
void testzdotc() {
Rcomplex zx[3], zy[3], ret_val;
zx[0].r = 1.0; zx[0].i = 0.0;
zx[1].r = 2.0; zx[0].i = 0.0;
zx[2].r = 3.0; zx[0].i = 0.0;
zy[0].r = 1.0; zy[0].i = 0.0;
zy[1].r = 2.0; zy[0].i = 0.0;
zy[2].r = 3.0;
2020 Jul 15
2
Openblas?
Hello,
I thought that I should try openblas when building a CRAN package
containing lots of old (twentieth century) C-code with frequent calls to
blas and lapack routines. I have the following options on my Ubuntu
20.04 machine:
Selection Path Priority Status
------------------------------------------------------------
* 0
2009 Jun 12
1
Can't get F77_CALL(dgemm) to work [SEC=UNCLASSIFIED]
Hi
I am new to writing C code and am trying to write an R extension in C. I
have hit a wall with F77_CALL(dgemm) in that it produces wrong results.
The code below is a simplified example that multiplies the matrices Ab and
Bm to give Cm. The results below show clearly that Cm is wrong.
Am=
1 2 3 4
5 6 7 8
9 10 11 12
13 14 15 16
17 18 19 20
Bm=
1 1 1
1 1